SYSTÈMES DYNAMIQUES DIFFÉRENTIABLES

SYSTÈMES DYNAMIQUES DIFFÉRENTIABLES

Sans doute née avec le mémoire que Poincaré écrivit en 1881 «sur les courbes définies par des équations différentielles», où l’étude quantitative (analytique) locale des équations différentielles dans le champ complexe est remplacée par leur étude qualitative (géométrique) globale dans le champ réel, la théorie des systèmes dynamiques a pour ambition une description aussi précise que possible du comportement asymptotique (en temps infini ) des solutions des (de la plupart des, de certaines) équations différentielles; la mécanique céleste et les accélérateurs de particules sont des exemples où l’intérêt de tels problèmes est manifeste, mais aussi les oscillations électriques et la turbulence hydrodynamique encore si mal comprise. Il est remarquable que beaucoup des idées fondatrices viennent d’un très petit nombre de tels exemples, en particulier le problème des trois corps et les équations du type Van der Pol.

Nous n’avons pas dans cet article la prétention de déployer exhaustivement la théorie, mais simplement d’introduire le lecteur à quelques points de vue et, dans une petite mesure, à quelques méthodes de ce beau domaine des mathématiques. Les premiers chapitres décrivent en assez grand détail des exemples tournant autour du pendule et introduisent des comportements asymptotiques de plus en plus complexes dont l’analyse nécessite les concepts plus abstraits de la dernière partie.

1. Le pendule sans frottement, un système hamiltonien

L’étude qualitative de l’ensemble des mouvements possibles d’un pendule rigide dans un plan vertical offre une illustration simple de la notion de système dynamique. Si l’expérience «montre» qu’il ne suffit pas de se donner la position initiale du pendule pour caractériser son mouvement, elle suggère par contre que la connaissance à un instant donné de sa position et de sa vitesse implique celle de son mouvement à tout instant ultérieur (et passé, mais il s’agit là d’expérience de «pensée»).

La position est repérée par un angle ou, ce qui revient au même, par un point x du cercle S¹ de rayon 1 centré sur le point de suspension et contenu dans le plan des oscillations; la vitesse angulaire est repérée par un nombre réel y (fig. 1). L’ensemble des couples (x , y ) constitue l’espace de phase du pendule: c’est un cylindre (fig. 2), produit topologique d’un cercle par une droite (on note ce produit S¹ 憐 R).

Chaque mouvement est représenté par une courbe orientée (éventuellement réduite à un point si le pendule reste immobile) tracée sur le cylindre et paramétrée par le temps qu’on suppose s’étendre de 漣 秊 à + 秊. Chaque courbe décrit une évolution au cours du temps de la position et de la vitesse du pendule, mais une même courbe géométrique correspond à une infinité de mouvements déduits les uns des autres par une translation du temps; on appelle courbes intégrales ou orbites ces courbes géométriques orientées: deux d’entre elles ne peuvent s’intersecter en vertu des remarques qui précèdent, c’est précisément ce qu’on a gagné à remplacer l’espace des positions (cercle) par l’espace de phase (cylindre).

La figure 3 représente quelques expériences et leur traduction dans l’espace de phase: la première rappelle que si, à un instant donné, x est égal à 0 (position la plus basse) ainsi que y (vitesse nulle) il en est de même à tout instant ultérieur (et passé): on dit que (0, 0) est un équilibre ; la deuxième concerne les petites oscillations autour de la position d’équilibre: le pendule étant supposé sans frottement, la conservation de l’énergie implique des oscillations entre la position de départ (à vitesse nulle) et la position symétrique par rapport au point de suspension. Dans l’espace de phase, ces oscillations correspondent à des orbites périodiques (courbes intégrales homéomorphes à un cercle). Dans la troisième, enfin, la vitesse initiale est suffisante pour que l’énergie totale dépasse l’énergie potentielle maximale: la vitesse ne s’annule jamais et, l’énergie se conservant, le pendule revient après chaque tour en un point donné avec une vitesse identique.

Il est maintenant possible d’esquisser le portrait de phase (ensemble des courbes intégrales); le dessin sur le cylindre étant difficile, on a considéré ce dernier comme un rouleau d’imprimerie, ce qui permet de représenter les courbes intégrales dans le plan: chaque figure tracée sur le cylindre donnera par impression une figure plane périodique en x (fig. 4; on notera que pour éviter les inconvénients de l’imprimerie on a rétabli le sens de l’axe des x en retournant la feuille). Pour les lecteurs mathématiciens, le cylindre n’est autre que le fibré tangent du cercle, et le déroulement représente le passage au revêtement universel.

Il semble difficile d’analyser par une expérience l’allure des courbes intégrales au voisinage du point ( 神, 0) du cylindre: en effet, ce point correspond à la position d’équilibre supérieure, qui est instable dans le sens où la moindre perturbation amène le pendule à occuper des positions très éloignées.

Le temps est venu de remarquer que la figure 4 n’est autre qu’une résolution presque complète de l’équation de Newton, づ = m ゔ, c’est-à-dire (fig. 5):

qu’on écrira:

On a noté l la distance du centre de gravité au point de suspension, m la masse, et déroulé le cercle (comme plus haut le cylindre) en considérant x comme une variable réelle.

L’équation (N¹) s’écrit encore comme équation différentielle autonome du premier ordre dans R²:

et s’interprète comme un champ de vecteurs X dans R² (indépendant du temps) périodique en x (de période 2 神), c’est-à-dire un champ de vecteurs tangents au cylindre (fig. 6).

Les solutions de (P¹) sont les courbes paramétrées (x (t ), y (t )) qui ont pour dérivée (vecteur vitesse) au temps t le vecteur (y (t ), face=F0019 漣 諸² sinx (t )): ce sont les courbes intégrales dont certaines ont été esquissées sur la figure 4. Les équilibres (courbes intégrales réduites à un point) sont encore appelés singularités : ce sont les points (x , y ) en lesquels le vecteur X(x , y ) s’annule. Le fait que les courbes intégrales remplissent le plan (le cylindre) et ne se rencontrent pas deux à deux est la traduction géométrique du théorème fondamental d’existence et d’unicité des solutions d’équations différentielles.

Il se trouve que l’absence supposée de frottement permet de trouver très facilement les courbes intégrales:

L’équation (P¹) s’écrit en effet sous la forme hamiltonienne :

où H(x , y ) = (y ²/2) + 諸²(1 漣 cosx ) est proportionnelle à l’énergie totale du pendule dans la position x , avec la vitesse y (H est périodique en x et définit donc une fonction sur le cylindre, appelée hamiltonien ).

On vérifie immédiatement que, le long d’une courbe intégrale (x (t ), y (t )), H reste constante (conservation de l’énergie): ces courbes sont donc contenues dans les courbes de niveau de H, ce qui permet de tracer la figure 7.

Il peut être agréable de remarquer que le plongement F du cylindre dans R³ défini par F(x , y ) = (sinx , 漣 y , (y ²/2) + 諸²(1 漣 cosx )) transforme H en la hauteur au-dessus du plan horizontal (fig. 7 bis où la forme des courbes d’énergie constante devient évidente). Pour une autre représentation, cf. fig. 33 bis .

On constate en particulier qu’il existe deux courbes intégrales issues du point ( 神, 0) sur le cylindre (le point lui-même n’appartient pas à ces courbes): celles-ci traduisent l’existence pour chaque position x 神 du pendule de deux vitesses initiales y ¹, y ² de signe contraire telles que lorsque le temps devient très grand le pendule tende vers la position d’équilibre supérieur (fig. 7). Il apparaît sur la figure 7 (et à l’expérience) que les deux équilibres (0, 0) et ( 神, 0) sont de nature différente: seul (0, 0) est stable au sens de Liapunov (le nom de Liapunov reviendra plusieurs fois dans cet article: ses contributions – dès sa thèse en 1892 – à la théorie de la stabilité sont fondamentales), ce qui signifie que si (x ⁰, y ⁰) est assez proche de (0, 0) la portion de la courbe intégrale passant par (x ⁰, y ⁰) et correspondant à des temps ultérieurs reste assez proche de (0, 0).

On aurait pu s’en rendre compte par une étude purement locale en linéarisant le champ de vecteurs (l’équation différentielle) P¹ au voisinage de chacun d’eux:

Au voisinage de (0, 0), l’équation (P¹) s’écrit:

où 見(x , y ) = 0 et 廓(x , y ) = 漣 諸²(sin x 漣 x ) sont du troisième ordre en (x , y ) et donc petits par rapport au premier terme lorsque (x , y ) est proche de (0, 0).

L’équation différentielle linéaire dans R² obtenue en supprimant 見 et 廓 est appelée linéarisée de (P¹) en (0, 0): remarquer qu’elle n’est plus périodique en x ; ses valeurs propres 梁 i 諸 sont imaginaires et ses courbes intégrales sont des ellipses (fig. 8 a; on a indiqué le spectre imaginaire de la matrice), courbes de niveau du hamiltonien:

qui, ce n’est pas un hasard, n’est autre que la partie quadratique du développement de Taylor de H en (0, 0). On dit que (0, 0) est un centre .

Au voisinage de ( 神, 0), l’équation (P¹) s’écrit:

où じ = x 漣 神, et où 見( じ , y ) = 0, 廓( じ , y ) = 諸² (sin じ 漣 じ ) sont du troisième ordre en ( じ , y ) et donc petits par rapport au premier terme lorsque (x , y ) est proche de ( 神, 0).

Les valeurs propres de la linéarisée de (P¹) en ( 神, 0) sont maintenant réelles (égales à 梁 諸) et les courbes intégrales sont des hyperboles (fig. 8 b; on a indiqué le spectre réel de la matrice), courbes de niveau du hamiltonien H^u ( じ , y ) = y ²2 漣 諸² じ ²2 qui, ce n’est toujours pas un hasard, coïncide avec la partie quadratique du développement de Taylor de H en ( 神, 0). On dit que ( 神, 0) est un col (encore appelé selle ). On notera le rôle particulier joué par la courbe de niveau dégénérée H^u ( じ , y ) = 0.

Dans les deux cas, on note une grande ressemblance qualitative entre le portrait de phase de (P¹) au voisinage d’un de ses équilibres et le portrait de phase des équations linéarisées en ces équilibres. Plus précisément, on déduit du lemme de Morse ^{[cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES (la théorie mathématique et ses applications)]} l’existence au voisinage de (0, 0) (resp. de ( 神, 0)) d’un changement de coordonnées différentiable transformant les lignes de niveau de H en les lignes de niveau de H^s (resp. de H^u ); il y a bien entendu des différences quantitatives: par exemple, dans le cas elliptique (valeurs propres purement imaginaires), la période des oscillations dépend de leur amplitude (c’est-à-dire de la courbe intégrale considérée) pour l’équation (P¹), mais n’en dépend pas pour l’équation linéarisée.

Nous verrons dans le chapitre suivant que cette ressemblance est un phénomène très général dans le cas hyperbolique (valeurs propres non purement imaginaires) et est au contraire forcée par le caractère hamiltonien de (P¹) dans le cas elliptique.

Au début de ce chapitre, nous avons dit que le point de vue des systèmes dynamiques consistait en la considération simultanée de l’ensemble des mouvements possibles, c’est-à-dire de l’ensemble des courbes intégrales. Pour formaliser cela, introduisons le flot de l’équation (P¹) (ou du champ de vecteurs X): c’est la famille de difféomorphismes ( 﨏^t )^t size=1^捻R du cylindre (bijections du cylindre différentiables ainsi que leurs inverses) définie par 﨏^t (x ⁰, y ⁰) = (x , y ) si x et y sont respectivement la position et la vitesse au temps t du pendule quand, à l’instant 0, celui-ci a la position x ⁰ et la vitesse y ⁰. On peut dire que ( 﨏^t )^t size=1^捻R est la solution de (P¹), puisque cette définition équivaut à d 﨏^t dt (x , y ) = X( 﨏^t (x , y )). Du caractère autonome de l’équation (champ de vecteurs indépendant du temps) on déduit immédiatement que ( 﨏^t )^t size=1^捻R est un groupe à un paramètre de difféomorphismes , c’est-à-dire vérifie 﨏^{s +t} = 﨏^s 獵 﨏^t quels que soient les réels s et t ; en particulier, 﨏^-t = ( 﨏^t )^-1. On dit encore que ( 﨏^t )^t size=1^捻R est le système dynamique associé à l’équation (ou au champ de vecteurs). Le fait que les 﨏^t soient différentiables est une conséquence des théorèmes classiques de dépendance des solutions d’équations différentielles en fonction des données initiales. Sur une variété compacte, il y a équivalence entre les notions de champ de vecteurs et de flot; dans le cas général, il faut éventuellement modifier la longueur des vecteurs (et non leur direction, ce qui revient à changer la loi du temps) pour que les solutions existent pour tout temps t et définissent donc un flot (et cela même pour des problèmes naturels: Gerver a donné récemment (1983) une démonstration heuristique de l’existence de solutions sans collision du «problème des cinq corps» telles que tous les corps partent à l’infini en un temps fini). Ici l’existence du hamiltonien H fournit une propriété supplémentaire: chacun des difféomorphismes 﨏^t conserve les aires (aire 﨏^t (U) = aire U, cf. fig. 9 où s et t sont positifs), ce qui équivaut au niveau différentiel à:

et explique que les matrices des équations linéarisées en (0, 0) et ( 神, 0) soient de trace nulle (le flot de l’équation linéarisée en un équilibre n’est autre que le flot des dérivées des 﨏^t en cet équilibre).

Finissons ce chapitre par une petite excursion non autonome, introduisant les variables action-angle et montrant l’invariance adiabatique de l’action: posons x = 1 諸 練 連2I 諸 sin 﨏, y = 連2I 諸 cos 﨏, c’est-à-dire:

﨏 = Arc tg 諸x y . L’équation (L): d ²x dt ² + 諸²x = 0, linéarisée en (0, 0) de (P¹), prend la forme:

particulièrement simple (on pourrait faire quelque chose d’analogue pour (P¹) mais les formules ne seraient plus aussi explicites).

Remarquons que la transformation (x , y ) 料 (I, 﨏) conserve les aires; les varia-bles (I, 﨏) sont les variables (action, angle) de la mécanique hamiltonienne ^{[cf. MÉCANIQUE ANALYTIQUE]}, en particulier l’action reste constante au cours d’un mouvement.

Supposons maintenant que la longueur du pendule (resp. le champ de gravitation) varie lentement en fonction du temps; plus précisément, considérons l’équation (face="EU Upmacr" 籠):

où 諸() est une fonction deux fois continûment différentiable de. On peut montrer que, de façon surprenante, si (x (t ), y (t )) est une solution de (face="EU Upmacr" 籠), et si

est le rapport de l’énergie à la fréquence (action) à l’instant t , I(t ) reste presque constant sur des intervalles de temps très grands (de l’ordre de 1/ 﨎); on dit que I est un invariant adiabatique de (L).

La démonstration repose sur la technique de moyennisation introduite par Van der Pol et intensivement utilisée (et perfectionnée) par Krilov et Bogolioubov. On écrit (face="EU Upmacr" 籠) dans les variables (I, 﨏):

et on la compare à l’équation moyennée en 﨏:

L’idée, que l’on retrouvera au chapitre 3 sous le déguisement des formes normales , est de simplifier l’équation (face="EU Upmacr" 籠) par un changement de variables de la forme P = I + 﨎k (I, 﨏), où k est une fonction périodique en 﨏; on obtient:

Il existe k périodique en 﨏 telle que le terme en 﨎 s’annule car g est de moyenne nulle en 﨎 (sinon il faudrait garder 﨎 凜 (I, )). Pour un tel choix, qui fournit un vrai changement de variables si 﨎 est assez petit, P ne varie que de O( 﨎) sur un temps de l’ordre de 1/ 﨎, ce qui entraîne le même résultat pour I. L’équation (face="EU Upmacr" 聾) décrit l’évolution systématique, le passage à (face="EU Upmacr" 籠) amène des oscillations autour de cette évolution moyenne (fig. 10).

Nous reviendrons sur cet exemple au chapitre 4 à propos de la résonance paramétrique . Ce procédé de moyennisation est très utilisé en mécanique céleste, en physique, en ingénierie, souvent dans des circonstances où il n’est pas justifié. Quant à l’invariance adiabatique de l’action, elle a été mise en avant en 1911 par Lorentz et Einstein pour expliquer la constance surprenante du rapport énergie-fréquence de radiation d’un atome, les champs électromagnétiques extérieurs variant lentement par rapport à la fréquence d’oscillation de l’atome.

2. Le pendule avec frottement linéaire, un système structurellement stable

Nous avons pu déterminer exactement les courbes intégrales de l’équation du pendule sans frottement en utilisant la conservation de la fonction H le long de ces courbes.

Si on suppose maintenant, ce qui est plus réaliste, que le pendule possède un certain frottement, la détermination explicite de ces courbes intégrales se complique énormément mais, au moins dans le cas d’un frottement dépendant linéairement de la vitesse, leur allure qualitative peut être encore facilement décrite en mêlant analyse locale par linéarisation aux équilibres et propriétés globales de la fonction H.

Plus précisément, si le frottement est proportionnel à la vitesse, l’équation de Newton devient:

ou encore:

L’allure des courbes intégrales est en grande partie dictée par l’observation suivante: la fonction H définie à propos de (P¹) est non croissante le long d’une courbe intégrale de (P²) (dissipation de l’énergie); elle est même strictement décroissante en tous les points où y 0 (on dit que H est une fonction de Liapunov pour (P²) dans le complémentaire de y = 0). En effet,

Corrélativement, la divergence du champ de vecteurs défini par (P²) est 漣 福: le flot associé «diminue» les aires.

D’autre part, les deux équilibres sont maintenant hyperboliques: au voisinage de (0, 0), (P²) s’écrit:

alors qu’au voisinage de ( 神, 0), dans les coordonnées じ , y introduites à propos de (P¹), elle s’écrit:

Le portrait de phase de l’équation linéarisée en (0, 0) dépend de la taille de 福: la figure 11 résume la situation en présentant simultanément le portrait de phase, le nom de la singularité et, dans un encart, le spectre de la matrice:

notons que les figures 8 et 11 étaient connues du célèbre mécanicien russe N. E. Joukovsky (elles apparaissent explicitement dans sa Thèse sur la cinématique des corps liquides publiée en 1876).

Au contraire, le portrait de phase de l’équation linéarisée en ( 神, 0) est qualitativement indépendant de la taille de 福 et ressemble à celui obtenu pour 福 = 0 (le spectre de la matrice:

est toujours formé de deux valeurs propres réelles de signes opposés).

On peut maintenant donner une idée assez précise du portrait de phase de (P²) suivant les valeurs de 福:

– Le théorème de linéarisation de Hartman (1963) sur lequel nous reviendrons affirme en effet l’existence d’un changement de coordonnées continûment dérivable défini au voisinage de (0, 0) (resp. de ( 神, 0)) transformant localement (P²) en l’équation linéarisée en (0, 0) (resp. en ( 神, 0)). Il faut noter qu’un changement de coordonnées qui ne serait que continu serait insuffisant pour distinguer entre les différentes situations représentées sur la figure 11; cependant, une telle distinction est nécessaire si l’on veut rendre compte de l’expérience; si 福 est petit, l’amortissement au voisinage de la position d’équilibre inférieure s’accompagne d’oscillations, si 福 est grand il n’y a pas d’oscillation (fig. 12). Dans les deux cas, on n’atteint la position d’équilibre qu’au bout d’un temps infini (ce qui, tout comme l’infinité d’oscillations, souffre d’un manque certain de réalisme, mais notre pendule est idéal...).

Nous venons en passant de prouver la stabilité asymptotique de l’équilibre (0, 0) de l’équation (P²): si (x ⁰, y ⁰) est assez proche de (0, 0) la portion de la courbe intégrale passant par (x ⁰, y ⁰) et correspondant à des temps ultérieurs reste assez proche de (0, 0) (stabilité au sens de Liapunov) et tend vers (0, 0) lorsque le temps tend vers l’infini ((0, 0) est un attracteur de l’équation (P²)).

En combinant ces informations locales avec celles qui sont données par la décroissance de H le long des courbes intégrales de (P²), on obtient les figures 13 (cas 0 麗 福 麗 2 諸), 14 (cas 福 = 2 諸) et 15 (cas 福 礪 2 諸).

Un mouvement typique du pendule amorti se lit sur ces figures: après un nombre fini de tours complets, d’autant plus grand que la vitesse initiale était grande, il y a amortissement au voisinage de la position d’équilibre inférieure; ce n’est que dans les cas (rares!) où le point initial (x ⁰, y ⁰) est situé sur la courbe W^s ( 神, 0) réunion de ( 神, 0) et des deux courbes intégrales tendant vers le point ( 神, 0) du cylindre que la position asymptotiquement atteinte sera la position d’équilibre supérieure. Le complémentaire de cette courbe est le bassin de l’attracteur (0, 0). Cette courbe et sa compagne définie de façon analogue en renversant le temps, appelées respectivement variété stable W^s ( 神, 0) et variété instable W^u ( 神, 0) (unstable en anglais) de l’équilibre hyperbolique ( 神, 0) ont une importance fondamentale pour l’étude géométrique de l’équation et nous les retrouverons, elles ou leurs avatars, tout au long de l’article. Elles sont définies par:

où ( 﨏^t )^t size=1^捻R est le flot de l’équation, défini au chapitre 1.

Pour l’équation linéarisée en ( 神, 0), dont le flot est le groupe à un paramètre de transformations linéaires:

dans les coordonnées ( じ , y ), ce sont les directions propres E^s et E^u de la matrice:

pour l’équation (P²) ce sont des courbes indéfiniment différentiables (en fait analytiques), même en ( 神, 0), où elles sont tangentes respectivement à E^s et E^u (fig. 16). Notons que cela ne découle pas du théorème cité de Hartman, puisque le changement de coordonnées fourni par ce dernier n’est en général pas deux fois différentiable (en dimension supérieure à deux, il n’est en général même pas différentiable mais seulement continu!).

L’existence de W^u ( 神, 0) manifeste l’instabilité asymptotique de l’équilibre supérieur: il existe des données initiales (x ⁰, y ⁰) aussi proches que l’on veut de ( 神, 0) telles que 﨏^t (x ⁰, y ⁰) sorte d’un voisinage fixé (assez petit) de ( 神, 0) pour tout t assez grand.

Revenant sur le théorème de linéarisation de Hartman, remarquons qu’il joue le rôle que tenait le lemme de Morse dans la situation hamiltonienne du chapitre 1. Si^s = (face=F0019 漣 1/2) ( 福 + 連 福² + 4 諸²) 麗 0 et^u = (1/2) (face=F0019 漣 福 + 連 福² + 4 諸²) 礪 0 sont les valeurs propres de (0 1 諸² 漣 福), un changement de coordonnées linéaire transforme en effet l’équation:

en:

c’est-à-dire en:

champ de gradient de la forme quadratique

De même que le lemme de Morse, le théorème de Hartman est valable en toutes dimensions: il affirme qu’en une singularité hyperbolique p d’un champ de vecteurs différentiable défini sur une variété différentiable, c’est-à-dire en un équilibre de l’équation associée tel que le spectre de l’équation linéarisée ne rencontre pas l’axe imaginaire, il existe un homéomorphisme local h qui, au voisinage de la singularité, transporte le flot ( 﨏^t )^t size=1^捻R de l’équation sur le flot (D 﨏^t (p ))^t size=1^捻R de l’équation linéarisée. On dit que h est une conjugaison topologique locale entre l’équation et sa linéarisée, ou encore une linéarisation topologique locale de l’équation. En dimension supérieure à deux, h n’est en général pas un difféomorphisme, comme le montre l’exemple:

au voisinage de (0, 0, 0). Cependant, s’il n’y a pas de résonance entre les valeurs propres¹, ...,ⁿ de l’équation linéarisée, c’est-à-dire pas de relation de la forme:

où les m ^j sont des entiers positifs ou nuls de somme supérieure à 1, et si le champ est indéfiniment différentiable, h est un difféomorphisme indéfiniment différentiable (théorème de Sternberg , 1958). Dans le cas de champs de vecteurs analytiques complexes, la démonstration de l’existence d’une linéarisation analytique remonte à la Thèse de Poincaré (1879) dans les cas non résonnants où le spectre est entièrement contenu dans un demi-plan ouvert de C, et à Siegel (1952) dans le cas contraire considérablement plus difficile où les conditions de non-résonance ne suffisent plus à cause des problèmes de «petits dénominateurs» que nous retrouverons dans un autre contexte au chapitre 4.

Le lemme de Morse étant à la base de la théorie de la stabilité des fonctions différentiables ^{[cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES (la théorie mathématique et ses applications)]}, on ne s’étonnera pas qu’il en soit de même du théorème de Hartman dans le monde des systèmes dynamiques: effectivement, deux champs de vecteurs linéaires hyperboliques suffisamment proches étant topologiquement conjugués (pas différentiablement, car les valeurs propres seraient des invariants), on voit que si l’on perturbe légèrement un champ de vecteurs X ayant une singularité hyperbolique en p , le nouveau champ Y aura également une singularité hyperbolique en un point q proche de p (cela découle du théorème des fonctions implicites) et il existera un homéomorphisme h d’un voisinage U de p sur un voisinage V de q transportant le flot de X sur le flot de Y, donc les courbes intégrales orientées de X sur les courbes intégrales orientées de Y (on dit que X est localement structurellement stable en p ). Pour que cela soit correct, il faut que la perturbation soit petite ainsi que toutes ses dérivées partielles d’ordre un.

Remarque. On peut se demander quel est ici l’analogue du théorème de stabilité locale des fonctions au voisinage d’un point où les

dérivées partielles ne sont pas toutes nulles, conséquence immédiate du théorème des fonctions implicites. C’est le théorème de la boîte de flot (flow-box ), également conséquence du théorème des fonctions implicites, qui affirme qu’au voisinage d’un point p où il ne s’annule pas, un champ de vecteurs est différentiablement conjugué à un champ de vecteurs constant (fig. 17) et est donc localement structurellement stable.

Nous pouvons introduire maintenant la notion globale de stabilité structurelle . Appelons deux équations (champs de vecteurs) topologiquement équivalentes s’il existe un homéomorphisme (global) transportant les courbes intégrales orientées (mais non paramétrées) de l’une sur les courbes intégrales orientées de l’autre (demander que les paramétrages soient préservés est trop exigeant, les périodes des orbites périodiques seraient par exemple des invariants). À partir de ce que nous venons de dire, il est possible de montrer que les diverses équations (P²) obtenues en faisant varier 福 礪 0 sont topologiquement équivalentes; cet exemple montre bien les limitations de cette notion, puisqu’une distinction telle qu’amortissement avec ou sans oscillation lui échappe. Appelons une équation structurellement stable si toute équation assez proche (au sens discuté plus haut dans les problèmes locaux) lui est topologiquement équivalente. Cette propriété concerne l’ensemble du portrait de phase (c’est-à-dire le système dynamique) et non pas une courbe intégrale (c’est-à-dire un mouvement) particulière. Sa définition remonte à Andronov et Pontryaguine qui parlaient de systèmes grossiers (1937). Les équations (P²) sont structurellement stables pour 福 0 alors que l’équation (P¹) ( 福 = 0) ne l’est pas. Lorsque dans une famille d’équations on en rencontre une qui n’est pas structurellement stable, on dit qu’il y a bifurcation : le portrait de phase change brusquement de nature; c’est le cas pour 福 = 0 dans la famille des équations (P²) où 福 passerait d’une valeur positive à une valeur négative.

On voit facilement qu’au niveau local une condition nécessaire de stabilité structurelle est que toutes les singularités de l’équation soient hyperboliques. Une première source de bifurcations est donc fournie par les singularités non hyperboliques (comme (0, 0) dans l’équation (P¹)): le problème du déploiement versel (cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES) se pose ici comme pour les fonctions, mais n’admet pas de solution aussi simple (voir les travaux de Takens, Arnold, Bogdanov), si ce n’est dans quelques cas particuliers (importants!) comme la bifurcation de Hopf que nous rencontrerons dans le chapitre suivant.

3. Bifurcation de cycles-limites et auto-oscillations

Si, au pendule avec frottement du chapitre précédent, on fournit un peu d’énergie lors de chaque passage par la position d’équilibre inférieure, on peut s’attendre à ce qu’une oscillation se maintienne sans amortissement.

Dans l’espace de phase, une telle oscillation correspond asymptotiquement à un cycle-limite stable, premier exemple d’un attracteur qui ne soit pas un équilibre. Il semble que ce soit lord Rayleigh qui ait, en 1877, écrit la première équation de ce type à propos de sa théorie des instruments à archet. L’équation (N³) correspondante est en rapport étroit avec la célèbre équation de Van der Pol (1920):

à la base de tant de développements dans la théorie des systèmes dynamiques; elle s’écrit:

c’est-à-dire:

Le (petit) paramètre 猪 est destiné à mettre en évidence l’apparition par bifurcation du cycle limite; quant à f size=1^猪 c’est une fonction périodique dont le graphe entre 漣 神 et + 神 est représenté sur la figure 18 pour diverses valeurs de 猪 proches de 0: en dehors d’un voisinage de 0, f size=1^猪(x ) = 福 comme dans le chapitre précédent; au voisinage de 0, on pourra supposer que f size=1^猪(x ) = x ² 漣 猪.

La figure 18 indique la forme du graphe de (qui coïncide avec le déploiement versel de la singularité (1/3)x ³ au voisinage de 0).

Si H est défini comme dans les chapitres précédents, on constate que:

le long d’une courbe intégrale (x (t ), y (t )) de (P³). Si 猪 諒 0, le comportement des courbes intégrales est donc qualitativement le même que dans le chapitre 2: l’équilibre (0, 0) est asymptotiquement stable; si 猪 礪 0, H ne décroît le long d’une courbe intégrale de (P³) qu’en dehors d’une étroite bande axée sur la génératrice x = 0 du cylindre.

Si | 猪| est assez petit, il est facile de montrer que le comportement qualitatif de la variété stable de l’équilibre ( 神, 0) ne change pas lorsque 猪 passe d’une valeur négative à une valeur positive. Pour étudier ce qui se passe au voisinage de (0, 0) il est commode de faire sur l’équation (P³) (déroulée dans R²) le changement de variables (non périodique en x ):

C’est l’analogue de celui qui permit à Liénard d’étudier l’équation de Van der Pol; on obtient l’équation:

dans laquelle on peut supposer que:

pourvu qu’on ne s’intéresse qu’à des morceaux de courbes intégrales pour lesquels |x | reste assez petit.

On notera (Q ³) l’équation ainsi obtenue:

.

Pour | 猪| 麗 2 諸 (ce qu’on suppose évidemment), les valeurs propres de l’équation linéarisée en (0, 0) sont:

et size=1^猪, associées par exemple aux vecteurs propres e size=1^猪 = ( size=1^猪, 漣 諸²) et 勒 size=1^猪. Dans la basee size=1^猪, 勒 size=1^猪 un vecteur réel (x , z ) a des coordonnées de la forme 﨡, 﨡, où:

Il suffit de connaître 﨡 pour repérer un tel vecteur, ce qui revient à choisir une identification de R² avec C. L’équation linéarisée s’écrit alors:

si on remplace 﨡 par 兀 = size=1^猪 﨡, l’équation (Q ³) devient:

Enfin on remarque, comme l’a fait Poincaré dans sa Thèse à propos du problème de linéarisation évoqué au chapitre 1, qu’un changement de variables bien choisi de la forme:

transforme cette équation en:

(c’est la présence de la résonance ⁰ = 2⁰ +⁰ qui empêche l’annulation du terme en 﨣| 﨣|²).

Si on oublie les termes d’ordre supérieur ou égal à quatre dans le développement de Taylor, on obtient une équation:

invariante par rotation (le remplacement de 﨣 par e ⁱ size=1^見 﨣 ne change pas la forme de l’équation); on dit que (face="EU Upmacr" 賂 ³) est une forme normale approchée de (Q ³). L’invariance rend évidente la détermination d’éventuelles orbites périodiques de (face="EU Upmacr" 賂 ³) entourant (0, 0): celles-ci ne peuvent être que des cercles centrés en (0, 0); en écrivant l’orthogonalité de 﨣 et d 﨣dt , on obtient | 﨣|² = 3 猪. Pour 猪 礪 0, (face="EU Upmacr" 賂 ³) possède donc une orbite périodique qui est un cercle de centre (0, 0) et de rayon 連3 猪 (dans le plan complexe de la variable 﨣).

La figure 19 indique les portraits de phase de (face="EU Upmacr" 賂 ³) pour diverses valeurs de 猪 proches de 0 (on a représenté ces portraits de phase dans l’espace produit R 憐 C pour faire apparaître la rapide croissance de la taille de l’orbite périodique en fonction du paramètre 﨣 et on a mis dans un encart le spectre de l’équation linéarisée en 﨣 = 0).

Pour 猪 麗 0, 﨣 = 0 est un attracteur ponctuel (équilibre asympotiquement stable); c’est encore vrai pour 猪 = 0, mais avec une attraction beaucoup moins forte (due seulement aux termes non linéaires du développement de Taylor); pour 猪 礪 0, 﨣 = 0 est devenu un répulseur , et c’est l’orbite périodique venant de naître qui attire asymptotiquement les courbes intégrales (attracteur périodique ou cycle-limite ).

Il se trouve qu’on peut montrer que les termes négligés en passant de (Q ³) à (face="EU Upmacr" 賂 ³) ne changent pas le comportement qualitatif des courbes intégrales dans un voisinage de (0, 0) (uniforme en 猪 pourvu que | 猪| soit assez petit): cela vient essentiellement de ce que la force de l’attraction de l’orbite périodique de (face="EU Upmacr" 賂 ³) pour 猪 礪 0 est assez grande par rapport à la perturbation qui fait passer de (face="EU Upmacr" 賂 ³) à (Q ³); techniquement, on peut regarder comment se transforment par le générateur 﨏¹ du flot ( 﨏^t )^t size=1^捻R de (Q ³) des courbes fermées voisines de l’orbite périodique de (face="EU Upmacr" 賂 ³) et montrer que la transformation ainsi définie est une contraction dans un espace de courbes régulières bien choisi; l’unique point fixe de cette contraction est une orbite périodique de (Q ³): cette méthode de transformées de graphes remonte à J. Hadamard et joue un rôle important dans la recherche de sous-variétés invariantes normalement hyperboliques d’un flot (c’est-à-dire préservées par tous les 﨏^t et telles que les répulsions ou attractions dans les directions normales l’emportent sur celles qui existent éventuellement dans les directions tangentielles).

Si 猪 礪 0 est assez petit, cette orbite périodique de (Q ³) se trouve tout entière dans une région où (Q ³) et (Q³) coïncident, ce qui montre l’existence d’une orbite périodique pour (Q³), donc pour (P³).

En travaillant un peu plus, on prouve que (Q³) (donc (P³)) ne possède pas d’autre orbite périodique (cela découle de ce qui précède au voisinage de (0, 0) et plus généralement d’un raisonnement de monotonie analogue à celui qu’on fait classiquement pour l’équation de Van der Pol, cf. Lefschetz, Differential Equations, Geometric Theory ).

On a esquissé sur la figure 20 le portrait de phase de l’équation (P³) pour 猪 礪 0 petit: pour presque toute donnée initiale, le mouvement asymptotique est une oscillation non amortie de période proche de 2 神諸 et dont l’amplitude est bien déterminée, contrairement à ce qui se passe pour le pendule sans frottement.

L’apparition d’un cycle limite par déstabilisation d’un équilibre est appelée dans la littérature bifurcation de Hopf (bifurcation d’Andronov chez certains auteurs russes); elle et ses généralisations (apparition de tores invariants sur lesquels les courbes intégrales sont presque quasi périodiques) jouent un grand rôle dans les modèles décrivant l’établissement de la turbulence hydrodynamique à partir d’un régime laminaire lorsque varie un certain paramètre (le nombre de Reynolds par exemple).

Remarquons qu’il s’agit d’un problème de perturbation, de caractère local, donc susceptible d’une approche analytique. Peut-être est-il bon de comprendre exactement ce qu’on perturbe et, pour cela, de contempler un instant la situation dégénérée dans laquelle f size=1^猪(x ) est remplacé au voisinage de 0 par la fonction constante égale à 漣 猪, tout en valant 福 hors d’un voisinage un peu plus grand. L’équation (face="EU Upmacr" 賂 ³) est alors remplacée par l’équation:

et la figure 19 est remplacée par la figure 21 (qui géométriquement ne se distingue pas de la figure correspondant à la famille d’équations linéaires d 﨣dt = size=1^猪 﨣).

On constate que le rôle de la partie réelle du coefficient de 﨣| 﨣|² dans (face="EU Upmacr" 賂 ³) a été de courber en un paraboloïde dans la direction des 猪 positifs le plan 猪 = 0 entièrement rempli d’orbites périodiques de l’équation:

Ce sont en quelque sorte les orbites périodiques de cette équation particulière (ou même de l’équation linéarisée d 﨣dt = i 諸﨣) qui donnent naissance à l’ensemble des orbites périodiques des diverses équations (Q ³) pour 猪 positif.

Le théorème de bifurcation de Hopf assure que la situation décrite sur l’exemple se produit «en général» (au renversement simultané du temps et du paramètre près) lorsque, en fonction d’un paramètre 猪, le spectre de l’équation linéarisée en un équilibre comporte un couple de valeurs propres imaginaires conjuguées traversant l’axe imaginaire. Cela vient de ce qu’«en général» la partie réelle des termes non linéaires significatifs de la forme normale est non nulle. Récemment, Alexander et Yorke ont montré que l’existence dans R 憐 R² d’une famille à un paramètre d’orbites périodiques des équations de la famille (fig. 19 et 21) ne dépend que de l’hypothèse sur les spectres et est de nature purement topologique (cf. infra , Deuxième excursion topologique ; ce ne sera plus le cas pour la bifurcation d’une orbite périodique en un tore invariant); curieusement, la démonstration fait intervenir la fibration de Heinz Hopf (celui de la bifurcation étant Eberhart). Notons enfin que la théorie des variétés centrales fournit des généralisations de ce théorème à Rⁿ et même à des espaces de dimension infinie, ce qui permet de l’appliquer à l’équation de Navier-Stokes (voir par exemple les travaux de Iooss).

Remarque 1 . En fournissant de l’énergie de manière constante, on peut entretenir un mouvement de rotation du pendule.

L’équation correspondante:

intervient dans les modèles de moteurs électriques synchrones.

On peut montrer par exemple que, si a 礪 0 et 福 礪 0 sont assez petits et a assez grand par rapport à 福, le portrait de phase contient un unique cycle limite faisant le tour du cylindre.

Remarque 2 . La famille d’équations:

très proches de l’équation de Van der Pol, montre bien le caractère local des affirmations du théorème de bifurcations de Hopf. Le changement de variables:

la transforme en:

Une bifurcation de Hopf se produit pour 猪 = 0 au point x = 1, z = 漣 2/3. Si 﨎 礪 0 est assez proche de 0, la croissance en 連 猪 de la taille de l’orbite périodique qui naît est complètement occultée par une croissance infiniment plus rapide qui se produit peu après et qui est due à la forme perturbation singulière des équations (à la limite 﨎 = 0 on obtient des équations dont l’ordre s’est abaissé). La figure 22 correspond au cas où 﨎 = 0,01; à gauche, 猪 = 0,001 259 5487 et, à droite, 猪 = 0,001 259 5488. Ce phénomène, baptisé effet canard , n’a pu être bien élucidé que grâce aux méthodes de l’analyse non standard qui permettent de supposer 﨎 infiniment petit et non nul (travaux de Benoît-Callot-Diener; cf. ANALYSE NON STANDARD).

Trois excursions topologiques

Première excursion topologique: les théorèmes de Poincaré-Bendixson et Poincaré-Hopf

Revenons à l’équation (Q³): il eût été possible de prouver l’existence d’une orbite périodique pour 猪 礪 0 par d’autres méthodes analytiques, par exemple la méthode de moyennisation rencontrée au chapitre 1 ou la réduction à l’étude locale d’un difféomorphisme de R au voisinage d’un point fixe, mais aussi par une méthode purement topologique en constatant l’existence, autour de (0, 0), d’un anneau A (homéomorphe au produit topologique S¹ 憐 I du cercle par un intervalle) dans lequel «entre» le champ de vecteurs (Q³) (fig. 23).

On constate que chaque courbe intégrale issue d’un point du rayon I = A 惡(x , z ) 捻 R², x = 0, z 礪 0 revient ultérieurement couper ce rayon, définissant une application de premier retour (ou application de Poincaré ) 﨏 de I dans I. Les théorèmes généraux sur les équations différentielles prouvent la continuité (et même la différentiabilité) de cette application. Si a = a ⁰ est sur le bord inférieur de A, 﨏(a ⁰) 礪 a ⁰; si a = a ¹ est sur le bord supérieur, 﨏(a ¹) 麗 a ¹; un argument de continuité montre alors l’existence d’au moins un point de I où 﨏( ) = , c’est-à-dire d’une orbite périodique de (Q³) entièrement située dans l’anneau A. Remarquons que cette partie de l’argument se généralise aux dimensions supérieures, toute application continue d’un disque de dimension n dans lui-même ayant au moins un point fixe (théorème de Brouwer). Bien entendu, il faut des renseignements plus précis de monotonie sur 﨏, liés à la dimension un, pour conclure à l’unicité de l’orbite périodique.

Ce qui précède est un exemple très particulier des raisonnements qui mènent au théorème de Poincaré-Bendixson , spécifique des équations dans le plan, et dont l’inexistence en dimension supérieure nous vaudra de rencontrer des attracteurs étranges .

Pour énoncer ce théorème, il faut introduire la notion d’ensemble limite d’un flot ( 﨏^t )^t size=1^{捻 R} dans R² (qui s’étend immédiatement au cas d’une variété quelconque): pour tout point p de R², on note:

L size=1^諸(p ) (resp. L size=1^見(p )) est appelé ensemble 諸-limite (resp. ensemble 見-limite ) du point p . Le théorème affirme qu’étant donné un flot de difféomorphismes continûment dérivables de R², tout ensemble limite compact non vide et ne contenant pas d’équilibre ne peut être qu’une orbite périodique.

Plus généralement, s’il n’y a qu’un nombre fini de singularités (équilibres) chaque ensemble limite est ou bien une orbite périodique ou bien l’union de singularités et de courbes intégrales 﨏^t (x ) telles que leurs limites lim 﨏^t (x ) soient des singularités (fig. 24).

t 轢梁秊

Ce théorème n’est vrai que sur le plan et la sphère; sur le tore, par exemple, il existe des champs de vecteurs sans singularité dont toutes les courbes intégrales sont denses. Enfin, si l’on remplace l’ensemble des courbes intégrales d’un flot sur R² par un feuilletage de codimension 1 avec singularités (décomposition du plan en union disjointe de courbes et de points [singularités] ressemblant à un flot en dehors des singularités), on obtient des comportements récurrents plus complexes pouvant comporter des courbes denses; introduits par Thurston dans ses études du groupe des difféomorphismes d’une surface, de tels feuilletages ont été systématiquement étudiés par Rosenberg sous le nom suggestif de labyrinthes . On trouvera dans la Thèse de Levitt (1983) un panorama de l’étude dynamique des feuilletages des surfaces.

D’autres contraintes topologiques très générales existent, fondées sur l’importante notion d’indice d’une singularité isolée d’un champ de vecteurs X, introduite implicitement par Poincaré dès 1880 (avec un signe opposé au nôtre). Cet indice est le nombre algébrique de tours (dans le sens trigonométrique) qu’effectue la direction orientée du vecteur X(x , y ) quand (x , y ) parcourt dans le sens trigonométrique un cercle centré sur la singularité et assez petit pour que le vecteur ne s’annule en aucun de ses points.

La figure 25 donne quelques exemples de calculs d’indice.

Un peu d’homotopie ^{[cf. TOPOLOGIE - Topologie algébrique]} permet de montrer que si, par déformation continue du champ, on passe d’un groupe de singularités isolées à un autre, la somme de leurs indices reste constante. Il n’est alors pas étonnant que les indices des singularités hyperboliques soient 梁 1, car toute singularité isolée peut être transformée en une réunion d’un nombre fini de singularités hyperboliques par une petite déformation (fig. 26, inspirée par le déploiement versel de l’ombilic elliptique, cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES); ces dernières sont en quelque sorte les briques élémentaires qui, par déformation et coalescence, permettent de construire les singularités plus complexes (donc moins stables).

Soit 阮 un compact connexe de R² bordé par un nombre fini de courbes fermées (sous-variété compacte à bord). Supposons que le champ X «entre» dans 阮 (en particulier, X ne s’annule pas sur le bord de 阮: fig. 27, relative à une seule singularité d’indice 2 漣 3 = 漣 1).

Le théorème de Poincaré-Hopf assure dans ce cas particulier que, si X ne s’annule dans 阮 qu’en des points isolés, la somme des indices de ses singularités dans 阮 est égale à 2 漣 n , où n est le nombre de composantes dans le bord de 阮; en particulier, si 阮 n’est pas un anneau, X doit avoir au moins une singularité dans 阮. La notion d’indice a un sens pour un champ de vecteurs sur une surface (et même sur une variété quelconque). Le théorème général, démontré par Poincaré en 1880 pour la sphère, en 1885 pour les surfaces (sous des hypothèses de «généricité»), et par Hopf en 1926 pour les variétés de dimension quelconque, affirme l’égalité entre la somme des indices des singularités (supposées toutes isolées) d’un champ de vecteurs sur une variété compacte M et la caractéristique d’Euler-Poincaré 﨑(M) de M (on comparera à l’énoncé analogue pour les fonctions de Morse). Si M est une surface, 﨑(M) est la somme alternée du nombre de sommets moins le nombre d’arêtes plus le nombre de faces d’une triangulation de M; par exemple 﨑(S²) = 2 pour la sphère (regarder un tétraèdre): tout champ de vecteurs sur S² a donc au moins une singularité (on ne peut pas peigner une sphère sans faire un épi). Par contre 﨑(T²) = 0 pour le tore et il existe effectivement des champs de vecteurs sans singularité sur le tore; nous les retrouverons au chapitre 6.

Deuxième excursion topologique: indice de Fuller, théorème de Seifert, et démonstrations topologiques des théorèmes de bifurcation d’orbites périodiques

De l’indice d’un point singulier isolé d’un champ de vecteurs X à l’indice d’un point fixe isolé p d’une application continue f d’une variété M dans elle-même, il n’y a qu’un pas: on se ramène à M = Rⁿ par le choix d’une carte locale et on calcule l’indice du champ de vecteurs qui, à q voisin de p , associe le vecteur f (q ) 漣 q 捻 Rⁿ . Une homotopie évidente montre que l’indice de p comme singularité d’un champ X est égal à l’indice de p comme point fixe de 﨏^{t o} , ( 﨏^t )^t size=1^{捻 R} étant le flot associé à X et t ⁰ un réel non nul quelconque. Le théorème de Poincaré-Hopf se généralise: si tous les points fixes d’une application f d’une variété compacte sont isolés, la somme de leurs indices est égale au nombre de Lefschetz de f , intersection homologique dans M 憐 M du graphe de f et de la diagonale (graphe de l’identité), et donc égale à la caractéristique d’Euler 﨑(M) si f est homotope à l’identité.

Considérons maintenant le couple (T, face=F0021 六) formé d’un réel positif T et d’une orbite périodique isolée de période T du champ de vecteurs X sur la variété M. Soit un petit morceau d’hypersurface coupant transversalement 六 au point p , et soit f le difféomorphisme local de qui à q associe le premier point où l’orbite de q rencontre à nouveau (comparer à la fig. 28). Si ⁰ est la période minimale de 六 et si T = n ⁰ (n entier positif, bien entendu), on définit l’indice de Fuller i (T, face=F0021 六) comme étant l’indice du point fixe p de f divisé par n . Étant donné un sous-ensemble compact 行 de ]0, + 秊[ 憐 M dont le bord ne contient aucun couple (t , x ) tel que l’orbite de x soit périodique de période t , et dont l’intérieur ne contient que des orbites périodiques (T, face=F0021 六) isolées, la somme des indices de Fuller i (T, face=F0021 六) correspondants peut être définie homologiquement et est invariante par homotopie du champ X tant que la condition au bord reste remplie (théorème de Fuller , 1967). La nécessité de travailler dans le produit de M par ]0, + 秊[ vient de ce que, si la période des orbites périodiques considérées n’est pas bornée, on ne peut espérer aucune invariance, les orbites pouvant disparaître par bifurcation de Hopf (période minimale devenant brutalement 0) ou dans le bleu du ciel (période minimale tendant vers l’infini).

Ce théorème a permis à Fuller de redonner une démonstration du théorème de Seifert (1950) affirmant que tout champ de vecteurs sur la sphère de dimension trois suffisamment proche du champ de Hopf possède au moins une orbite périodique (les orbites du champ de Hopf sont les cercles de Clifford, fibres de la fibration de Hopf H de S³ sphère unité de C² sur la sphère de Riemann S² = P¹(C), définie par H(z ¹, z ²) = [z ¹, z ²] en coordonnées homogènes, c’est-à-dire les intersections avec S³ des droites complexes de C²).

On ne sait toujours pas si un champ de vecteurs C size=1^秊 sur la sphère de dimension trois doit nécessairement posséder une orbite périodique (cf. infra , chap. 6).

Le théorème de Fuller a été récemment (1977) utilisé par Chow, Mallet-Paret et Yorke pour donner une nouvelle démonstration topologique du théorème de bifurcation de Hopf, améliorant celle d’Alexander et Yorke à laquelle on a fait allusion plus haut; dans son esprit, cette démonstration est analogue aux démonstrations qu’a données Rabinowitz de théorèmes de bifurcation de points singuliers. L’existence dans l’espace produit R 憐 Rⁿ (R espace du paramètre 猪) d’une famille à un paramètre d’orbites périodiques des équations:

définies dans Rⁿ , est assurée dès qu’il existe une valeur propre i 諸(0) de la dérivée Df ⁰(0) qui vérifie r ⁺ 漣 r ^- 0, où r ⁺(r ^-) est le nombre de valeurs propres de Df size=1^猪(0) traversant l’axe imaginaire en i 諸(0) de gauche à droite (droite à gauche) lorsque 猪 passe d’une valeur négative à une valeur positive.

Ce type de théorème est assez puissant pour impliquer par exemple les généralisations qu’ont données Weinstein et Moser du théorème de Liapunov sur l’existence d’orbites périodiques dans tout voisinage d’un point singulier d’un champ de vecteurs hamiltonien non dégénéré X⁰ dont les surfaces d’énergie sont compactes (pour les définitions, cf. infra , chap. 4): l’idée est de construire, en ajoutant un gradient, une famille à un paramètre de champs de vecteurs X satisfaisant aux conditions du théorème et ne pouvant avoir d’orbite périodique que pour 猪 = 0 (c’est l’exact opposé du théorème de Hopf générique où, l’origine étant un attracteur ou un répulseur, aucune orbite périodique proche de 0 ne peut exister pour 猪 = 0).

Troisième excursion topologique: structures de contact et équations différentielles implicites

Voici enfin un résultat d’inexistence d’orbites périodiques, conséquence d’un théorème remarquable démontré récemment (1982) par Bennequin.

Soit:

une équation différentielle «implicite» (par opposition aux équations dans lesquelles dx /dt est «explicité» en fonction de t et x ). La recherche de solutions x = 祥(t ) d’une telle équation conduit évidemment à des fonctions 祥 multiformes et il est naturel de lui substituer la notion plus géométrique de solution obtenue en l’interprétant comme la trace sur la surface d’équation f (t , x , y ) = 0 de la structure de contact naturelle 暈 de l’espace des jets J¹(R, R) = R^{3 [cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES (la théorie mathématique et ses applications)]}. Cette idée remonte à Lie et Poincaré; c’est d’ailleurs à ce propos que Poincaré a été amené à considérer des champs de vecteurs sur d’autres surfaces que le plan ou la sphère. Les coordonnées sont ici t (source), x (but), y (dérivée dx /dt ), et 暈 est le champ des noyaux de la 1-forme différentielle ^{[cf. VARIÉTÉS DIFFÉRENTIABLES]} 諸(t , x , y ) = y dt 漣 dx (penser à y = dx /dt !), c’est-à-dire le champ de plans qui au point (t ⁰, x ⁰, y ⁰) associe le plan d’équation y ⁰(t 漣 t ⁰) 漣 (x 漣 x ⁰) = 0. Ladite trace est un champ de droites sur ,

orientable car l’est, dont les «courbes intégrales» (bien définies à la paramétrisation près) sont par définition les «solutions» de l’équation f (t , x , dx /dt ) = 0. (De même, la considération d’hypersurfaces dans J¹(Rⁿ , R) conduit à la théorie géométrique des équations aux dérivés partielles du premier ordre.)

La forme 諸 est une forme de contact , autrement dit 諸 廬 d 諸 ne s’annule en aucun point (forme volume); d’après le théorème de Frobenius , cela implique que le champ de plans 暈 n’est au voisinage d’aucun point de R³ le champ des plans tangents aux feuilles d’un feuilletage de R³ par des surfaces: les seules variétés intégrales de 暈 (c’est-à-dire tangentes en tout point aux plans de 暈) sont des courbes, appelées courbes de Legendre de la structure, dont les solutions de l’équation proposée sont des exemples. Le théorème de Bennequin affirme que, si est une sphère ou un plan (surface de genre 0) et si une courbe de Legendre tracée sur est homéomorphe à un cercle, il existe un point de la courbe où le plan de 暈 est tangent à (c’est-à-dire qu’aucune solution de l’équation f (t , x , dx /dt ) = 0 n’est homéomorphe à un cercle). Ce résultat est assez fort pour impliquer l’existence de structures de contact «exotiques» sur R³, c’est-à-dire de champs de plans localement (analogue du théorème de Darboux sur les formes symplectiques) mais non globalement difféomorphes à la structure 暈. Plus précisément, si 見 est une 1-forme différentielle suffisamment générale sur R³, l’ensemble des points (t , x , y ) où 見 廬 d 見 s’annule est une hypersurface régulière S.

Si, en (t ⁰, x ⁰, y ⁰), le champ 遼 des noyaux de 見 est tangent à S, le champ de directions défini par 遼 sur S possède en ce point une singularité qui, au premier ordre, est un col ou un centre. Dans ce dernier cas, la famille à un paramètre de champs obtenus en prenant les traces de 遼 sur une famille de surfaces parallèles à S présente, au passage de S, une bifurcation de Hopf; on déduit alors de ce qui précède que la restriction de 遼 à l’une des composantes connexes (difféomorphes à R³) de l’intersection du complémentaire de S avec un voisinage de (t ⁰, x ⁰, y ⁰) est une structure de contact exotique.

4. Perturbations périodiques d’un pendule sans frottement et difféomorphismes du plan préservant les aires

L’application au pendule des chapitres 1 et 2 d’une perturbation dépendant périodiquement du temps introduit une grande complexité de comportement, celle précisément dont Poincaré disait à la fin du troisième tome de ses Méthodes nouvelles de la mécanique céleste (chap. XXXIII, p. 389, 1899): «On sera frappé de la complexité de cette figure que je ne cherche même pas à tracer. Rien n’est plus propre à nous donner une idée de la complication du problème des trois corps et en général de tous les problèmes de dynamique où il n’y a pas d’intégrale uniforme et où les séries de Bohlin sont divergentes.» Dans les quatre-vingts ans écoulés depuis la publication de ces lignes, de grands pas ont été faits vers la compréhension de ces situations, certains prévus par Poincaré, d’autres qui l’auraient certainement surpris (mais auraient fait grand plaisir à Weierstrass); les plus marquants sont le développement, à la suite de Hadamard (1898) et Morse (1921), des méthodes de dynamique symbolique qui ont conduit à l’analyse par Birkhoff, puis Smale, de la dynamique associée à un point homocline, la démonstration par Birkhoff en 1913 du dernier théorème géométrique de Poincaré assurant l’existence de mouvements périodiques sous-harmoniques, enfin, dans les années cinquante et soixante, le théorème de Kolmogorov, Arnold, Moser (K.A.M.), affiné plus tard par Rüssmann et Herman, qui assure l’existence de familles de mouvements quasi périodiques et résout ainsi le problème de la stabilité des mouvements périodiques elliptiques d’un système hamiltonien à deux degrés de liberté (malheureusement celui qui décrit le système solaire en a plus de deux). La démonstration de ce dernier théorème est rendue très difficile par la présence de petits dénominateurs qui font diverger les classiques séries de perturbation (le premier théorème surmontant une difficulté de ce type a été démontré par Siegel en 1942). Un nouveau progrès dans le problème des mouvements quasi périodiques a été accompli récemment par Aubry et Mather à l’aide de méthodes variationnelles inspirées par la physique et rejoignant celles qui ont été utilisées dès 1932 par Hedlund à propos des géodésiques sur un tore riemannien. Notons que, s’il aurait surpris Poincaré et Birkhoff, le théorème de K.A.M. a également surpris la plupart des physiciens, puisqu’il restreint beaucoup les possibilités de validité de l’hypothèse ergodique de Boltzmann.

Les équations étudiées dans ce chapitre et le suivant sont de la forme:

où X = (x , y ) varie dans R² (ou plutôt dans le cylindre S¹ 憐 R), et f : R² 憐 R²R² vérifie f (X, t + T) = f (X, t ) (T est la période ). Plutôt qu’assimiler une telle équation à un champ de vecteurs dépendant du temps, nous l’écrivons comme un champ de vecteurs autonome (indépendant du temps) sur R² 憐 S^1T ou (S¹ 憐 R) 憐 S^1T (coordonnés X, 精), où S^1T = R/TZ est un cercle de longueur T:

Nous ne considérons que des équations telles que la courbe intégrale issue d’un point (X⁰, 精⁰) recoupe au bout d’un temps T le plan (ou le cylindre) 精 = 精⁰ (ce n’est pas automatique, certaines courbes intégrales pourraient se perdre à l’infini); si ( 﨏^t )^t size=1^{捻 R} est le flot de (), 﨏^T induit donc un difféomorphisme P du plan (cylindre) 精 = 精⁰ sur lui-même, appelé application de Poincaré ou application de premier retour. Si on change 精⁰, on remplace P par une application conjuguée h ^-1 獵 P 獵 h (h difféomorphisme de R²).

La connaissance de P équivaut à la connaissance de (): à un point fixe de P (c’est-à-dire P(X⁰) = X⁰) correspond une orbite périodique de période T de (), à un point périodique de période q de P (c’est-à-dire P^q (X) = X) correspond une orbite périodique de période q T de (), appelée sous-harmonique (fig. 28 dans le cas où q = 3).

Plus généralement, une courbe intégrale de () correspond à une orbite de P (c’est-à-dire la suite des Pⁿ (X), n 捻 Z, pour un X dans R²): on est donc amené à considérer P comme un flot à temps discret (remplacement de ( 﨏^t )^t size=1^{捻 R} par 﨏ⁿ = Pⁿ , n 捻 Z).

Toutes les notions définies pour les flots gardent un sens pour les flots à temps discret, en particulier les notions de stabilité: le spectre de l’équation linéarisée en un équilibre devient le spectre de la dérivée du difféomorphisme en un point fixe; le rôle de l’axe imaginaire est tenu maintenant par le cercle unité: il y a évidemment une exponentielle entre les deux comme on le voit en se rappelant que la solution générale de l’équation différentielle d Xdt = AX est 﨏^t (X⁰) = e ^{t A} 練 X⁰. On peut donc parler de point fixe hyperbolique (le spectre de la dérivée ne rencontre pas le cercle unité) ou elliptique (le spectre est sur le cercle unité). Les variétés stable et instable d’un équilibre hyperbolique, introduites au chapitre 2, ont leurs analogues pour les points fixes hyperboliques; si P(X⁰) = X⁰ et si le spectre de DP(X⁰) est formé de deux valeurs propres réelles¹,² vérifiant 0 麗 |¹| 麗 1 麗 |²|, la variété stable:

et la variété instable:

sont deux courbes lisses injectivement immergées dans R² et laissées invariantes par P: la démonstration (valable également en dimension quelconque) se fait par exemple en utilisant la méthode des transformées de graphes évoquée au chapitre 3 à propos de la bifurcation de Hopf. Comme pour les flots, le théorème de linéarisation de Hartman fournit un homéomorphisme local h conjuguant P et DP(X⁰) au voisinage d’un point fixe hyperbolique (c’est-à-dire P = h ^-1 獵 DP(X⁰) 獵 h ); dans le cas de R², cet homéomorphisme local est continûment différentiable; dans le cas de Rⁿ il est indéfiniment différentiable pourvu que certaines conditions de non-résonance soient satisfaites par les valeurs propres de P (théorème de Sternberg , qui a été récemment étendu aux germes d’action hyperbolique de Z^p 憐 R^q , par Dumortier-Roussarie et Chaperon). Dans le cadre holomorphe, les théorèmes de linéarisation de Poincaré et de Siegel ont également leur contrepartie pour les difféomorphismes. On comprend donc très bien la dynamique de P au voisinage d’un tel point fixe; par contre, la dynamique globale associée aux intersections des variétés stables et instables, que rien n’interdit, va amener la situation incroyablement complexe dont parlait Poincaré. Pour les points fixes elliptiques, il n’y a pas de théorème de linéarisation, et le théorème de K.A.M. montrera qu’au voisinage d’un tel point la dynamique de P est en général beaucoup plus riche que celle de DP(X⁰) bien qu’elle en retienne quelques traits.

Le saut dans la complexité des phénomènes décrits dans ce chapitre par rapport à ceux des chapitres précédents vient du saut (de deux à trois) de la dimension de l’espace de phase, saut qui se traduit ici par le remplacement d’un flot sur une surface par un difféomorphisme sur cette même surface (R² ou S¹ 憐 R): il n’existe pas, par exemple, de théorème jouant le rôle pour les difféomorphismes du plan du théorème de Poincaré-Bendixson; il serait intéressant de savoir s’il existe une différence aussi marquée dans le passage de la dimension n à la dimension n + 1 pour n 閭 3.

Nous traitons dans ce chapitre des systèmes conservatifs (pas de frottement), pour lesquels P conserve les aires, et dans le chapitre suivant des systèmes dissipatifs.

Et, pour commencer, remplaçons notre vieux pendule par une balançoire sur laquelle une petite fille fléchit et étend périodiquement ses petites jambes; l’équation du mouvement est de la forme:

où 諸(t )² = g l (t ) est périodique, de période T.

L’équation des «petits mouvements», c’est-à-dire l’équation linéarisée en x = dx /dt = 0, n’est autre que l’équation de Hill :

qui intervient comme équation aux variations dans l’étude de la stabilité d’orbites périodiques dans sa théorie de la Lune.

Il nous faut tout d’abord étudier la stabilité de la position d’équilibre inférieure, ce qui mettra en évidence le phénomène de la résonance paramétrique. Comme il s’agit de petits mouvements, il semble légitime de commencer par étudier (H); nous le ferons sur le cas particulier de l’équation de Mathieu :

Les équations (N⁴), (H), (M) correspondent à des champs de vecteurs dans R² 憐 S^1T dont la divergence est nulle; leurs flots conservent donc le volume, d’où l’on déduit immédiatement que P conserve les aires (det P = 1). De plus, la linéarité de (H) (ou (M)) implique la linéarité de P.

Pour 﨎 = 0, l’application de Poincaré P définie par le plan 精 = 0 est l’application linéaire de R² (coordonnées x , y = dx /dt ):

.

Sa trace, égale à 2cos T 諸, est en valeur absolue inférieure à 2 tant que 諸 n’est pas un multiple entier de 神/T; dans ce cas, le spectre de P est un couple de valeurs propres distinctes imaginaires conjuguées situées sur le cercle unité de C (P est elliptique ). Remarquons que si une application linéaire Q conserve les aires et est proche de P, son spectre aura la même propriété. On en déduit que, si 諸 n’est pas un multiple entier de 神/T, l’équilibre inférieur est stable au sens de Liapunov pour l’équation (M) tant que | 﨎| est assez petit (fig. 29).

Au contraire, P est l’identité si 諸 = 2 k 神T, k entier, moins l’identité si 諸 = (2 k + 1) 神T, k entier. Pour 﨎 différent de 0, il peut apparaître des intervalles de valeurs de 諸 proches des multiples entiers de 神/T pour lesquelles l’équilibre inférieur devient instable (apparition de valeurs propres réelles, 1/ pour l’application de Poincaré qui devient hyperbolique ): pour le voir, on évalue P en résolvant l’équation par la méthode de variation des constantes et on étudie sa trace; on peut aussi chercher 諸 en fonction de 﨎 de façon que (M) possède une solution périodique non triviale de période T (qui correspond à l’existence d’une valeur propre 1 pour P) ou 2 T (valeur propre 漣 1), ce qui fournit la frontière entre les régions elliptiques et hyperboliques (fig. 30; les zones hachurées sont celles où P est hyperbolique). On cherche de telles solutions sous forme de développements en 﨎 par une méthode d’identification terme à terme.

L’existence de ces zones hyperboliques montre que l’invariance adiabatique de l’action évoquée à la fin du chapitre 1 n’a pas toujours lieu sur un temps infini dans le cas d’une équation linéaire de la forme:

見 petit, bien qu’elle ait lieu sur un temps de l’ordre de 1/ 見; en effet, si ( 諸, 﨎) est dans une telle zone, l’instabilité a pour conséquence une variation indéfinie de l’action au cours du temps. Si 見 est petit, ces zones sont très nombreuses, puisqu’elles sont centrées sur les points 諸 = m 見/2 (m entier), 﨎 = 0. Arnold a montré qu’une conséquence du théorème de K.A.M. est en fait l’invariance adiabatique de l’action sur un temps infini dès que la non-linéarité de l’équation rend dépendante de l’amplitude la période des oscillations.

Nous revenons maintenant à l’équation non linéaire:

dont nous étions partis, avec par exemple:

Nous pourrions alternativement considérer une perturbation périodique du pendule conduisant à une équation de la forme:

où par exemple:

Lorsque 﨎 = 0, l’application de Poincaré P = P⁰ des champs de vecteurs () associés à (N⁴) ou (N ⁴) dans (S¹ 憐 R) 憐 S^1T coïncide avec le difféomorphisme 﨏^T de S¹ 憐 R, où ( 﨏^t )^t size=1^{捻 R} est le flot de l’équation:

Nous supposerons qu’il n’y a pas résonance paramétrique, c’est-à-dire que DP⁰(0, 0) est elliptique avec deux valeurs propres distinctes.

Si | 﨎| est petit, l’application de Poincaré P size=1^﨎 associée à (N⁴) ou (N ⁴) est une perturbation de P⁰ qui possède encore deux points fixes, l’un elliptique, l’autre hyperbolique: pour (N⁴) c’est évident, car les points fixes sont les mêmes que pour P⁰, pour (N ⁴) c’est une conséquence du théorème des fonctions implicites. Puisque P size=1^﨎 conserve les aires, un argument remontant à Poincaré montre que, dans le cylindre S¹ 憐 R, la variété stable et la variété instable du point fixe hyperbolique X^h se rencontrent; on peut montrer qu’en général (par exemple dans les équations que nous avons choisies) leur intersection est transversale (contrairement à ce qui se passe pour P⁰ où, si elles ne se confondaient pas, elles devraient être disjointes): ce sera par exemple une conséquence des méthodes de Melnikov que nous exposerons dans le chapitre 5 (les premiers exemples analytiques ont été donnés par Cherry). Dans la situation déroulée sur le plan, ce sont les variétés stables et instables de deux points fixes hyperboliques consécutifs (X^h )¹ et (X^h )² qui s’intersectent (fig. 31).

Dès qu’existe un point d’intersection transversale H (point homocline ), il en existe une infinité: à cause de l’invariance des variétés stables et instables, P size=1^﨎n (H) est en effet un point homocline pour tout n , tendant vers X^h , (resp. vers (X^h )² ou (X^h )¹ dans la situation déroulée où le terme hétérocline serait plus adéquat) si n tend vers + 秊 ou 漣 秊. Ce fait oblige W^u (X^h ) et W^s (X^h ) à se replier indéfiniment, ce qui les rend bien difficiles à dessiner globalement (fig. 32, d’après Zehnder) et justifie les exclamations de Poincaré.

Nous allons en déduire, suivant Smale, qu’il existe alors un ensemble de Cantor K plongé dans S¹ 憐 R, contenant H et X^h et invariant par une certaine puissance P size=1^﨎N de P size=1^﨎, sur lequel la dynamique induite par P size=1^﨎N peut être explicitement décrite et est en rapport étroit avec le jeu de «pile ou face», donc très erratique; en particulier, tout point de K (donc le point H) est limite de points périodiques de P size=1^﨎N (donc de P size=1^﨎) contenus dans K, résultat déjà connu de Birkhoff.

Bien entendu, en itérant P size=1^﨎, on peut, dans la démonstration, remplacer H par un point homocline 論 arbitrairement proche de X^h , assez proche en particulier pour être dans une région où P size=1^﨎 est presque linéaire (c’est important car en itérant on reviendra «souvent» dans cette région). L’idée est de considérer la restriction à un petit rectangle R contenant X^h d’un itéré P size=1^﨎N de P size=1^﨎 bien choisi, en particulier de telle façon que l’intersection entre R et P size=1^﨎N(R) contienne au moins deux morceaux (deux composantes connexes) comme sur la figure 33 où, en plus du déroulement, on a représenté l’identification du cylindre S¹ 憐 R au plan privé de l’origine par le difféomorphisme qui à (e ^ix , y ) 捻 S¹ 憐 R associe e ^{y +ix} 捻 C 漣0 = R² 漣0, ce qui revient à transformer les génératrices du cylindre en demi-droites du plan passant par l’origine (la figure 33 bis correspond au cas où 﨎 = 0 : on y retrouve dans cette nouvelle présentation du cylindre le portrait de phase du pendule sans frottement).

Dans ces nouvelles coordonnées, l’action de P size=1^﨎N sur R se voit beaucoup plus facilement: on l’a représentée schématiquement sur la figure 34.

Si l’orbite (pⁿ )ⁿ size=1^{捻 Z} d’un point p ⁰ sous l’action du difféomorphisme P size=1^﨎N reste dans R, p ⁿ est pour tout n dans l’intersection P size=1^﨎-N(R) 惡 R, formée des deux «rectangles verticaux» ⁺⁰ et ⁺¹ (la figure 35 est une représentation schématique). On peut donc associer à p ⁰ une suite (i ⁿ (p ⁰))ⁿ size=1^{捻 Z} dont tous les termes sont égaux à 0 ou 1 par la règle p ⁿ 捻 ^{+i n (p 0)}. Cette suite décrit l’«histoire» de l’orbite de p ⁰. Étant donné une suite finie 見⁰, ..., 見^{n -1} de zéros et de uns, l’ensemble ⁺ size=1^見0 ... size=1^{見n -1} des points p pour lesquels p 捻 ⁺ size=1^見0, (P size=1^﨎N)(p ) 捻 ⁺ size=1^見1, ..., (P size=1^﨎N)^{n -1}(p ) 捻 ⁺ size=1^{見n -1} est un «rectangle vertical» dont l’épaisseur tend vers 0 lorsque n tend vers + 秊; de même, l’ensemble ^- size=1^見-n ... size=1^見-1 des p pour lesquels (P size=1^﨎N)^-1(p ) 捻 ⁺ size=1^見-1, ..., (P size=1^﨎N)^-n (p ) 捻 ⁺ size=1^見-n est un «rectangle horizontal» dont l’épaisseur tend vers 0 lorsque n tend vers + 秊. L’intersection size=1^見-n ... size=1^{見-1 ;} size=1^見0 ... size=1^{見n -1} = ^- size=1^{見n -1} ... size=1^見-1 惡 ⁺ size=1^見0 ... size=1^{見n -1} est donc un «carré» dont les côtés tendent vers 0 lorsque n tend vers + 秊 (fig. 35). On en déduit qu’un point p ⁰ dont l’orbite reste dans R est entièrement déterminé par la suite (i ⁿ (p ⁰))ⁿ size=1^{捻 Z} et que, réciproquement, à toute telle suite de zéros et de uns correspond un tel point p ⁰.

L’ensemble des suites ( 見ⁿ )ⁿ size=1^{捻 Z} de zéros et de uns, que l’on note0, 1^Z, peut être muni de la «topologie produit»: deux suites sont proches si elles coïncident sur un segment fini assez long. L’espace topologique obtenu est homéomorphe au classique ensemble de Cantor qu’on obtient en retranchant d’un intervalle le tiers médian et en recommençant cette opération indéfiniment sur les intervalles restant (fig. 36).

L’application 崙 de0, 1^Z dans R définie en associant à la suite ( 見ⁿ )ⁿ size=1^{捻 Z} l’unique point p ⁰ de R qui vérifie (P size=1^﨎N)ⁿ (p ⁰) 捻 ⁺ size=1^見n pour tout n dans Z est un homéomorphisme sur son image (l’inverse est p ⁰ 料 (i ⁿ (p ⁰)ⁿ size=1^{捻 Z})) qui conjugue P size=1^﨎N au décalage de virgule (Shift) 靖: face=F00190, 1^Z索0, 1^Z défini par 靖(( 見n )ⁿ size=1^{捻 Z}) = ( 廓ⁿ )ⁿ size=1^{捻 Z}, 廓ⁿ = 見^{n +1}. Autrement dit, 崙^-1 獵 P size=1^﨎N 獵 崙 = 靖; l’étude des orbites de P size=1^﨎N qui restent dans R se ramène donc à celle de 靖 qui est à la base de la formulation probabiliste du jeu de «pile ou face»: c’est ce qu’on appelle faire de la dynamique symbolique .

Pour montrer ce qu’on a gagné, notons seulement qu’un point 神 = ( 見n )ⁿ size=1^{捻 Z} de0, 1^Z est périodique (face=F0019 說k , 靖^k ( 神) = 神) si et seulement si la suite 見ⁿ est périodique ( 見^{n +k} = 見ⁿ pour tout n ); on en déduit l’existence d’une infinité d’orbites périodiques de P size=1^﨎N, à savoir celles des points de la forme 崙( 神) (la méthode de Melnikov, exposée dans le chapitre suivant, peut également servir à prouver l’existence d’orbites périodiques de longues périodes de P size=1^﨎 dans cette région de S¹ 憐 R).

Cette étude a été formalisée par Smale en 1965 dans un exemple presque identique, devenu célèbre sous le nom de fer à cheval , et c’est ce nom que nous utiliserons pour signifier une situation analogue à celle qui est décrite par la figure 34; à l’origine, il y a les travaux de Cartwright-Littlewood et Levinson dont nous parlerons au chapitre suivant, premiers exemples d’équations différentielles structurellement stables ayant une infinité d’orbites périodiques. La figure 37 est une (difficile!) tentative de représenter l’ensemble K = 崙(face=F00190, 1^Z) 說 R (qui est invariant par P size=1^﨎N et contient X^h et 論). Elle montre que K, donc0, 1^Z, est homéomorphe à l’ensemble produit de0, 1^Z par lui-même; cette propriété de l’ensemble de Cantor d’ignorer les dimensions (il est homéomorphe à son carré) est la raison profonde de l’existence de courbes de Peano (courbes continues remplissant tout le plan).

Étudions maintenant la dynamique de P size=1^﨎 au voisinage de son point fixe elliptique X^e , proche de (0, 0): celle-ci n’est pas facile à deviner; en particulier il n’est pas clair a priori qu’il y ait stabilité au sens de Liapunov (définition analogue pour les difféomorphismes à celle qui a été donnée pour les flots au chapitre 1). Comme dans le chapitre 3, on commence par écrire l’équation sous une forme où la géométrie sous-jacente apparaît mieux; lorsque 﨎 = 0, (N⁴) et (N ⁴) se réduisent à (N¹) et il suffit de reprendre l’équation (Q ³) dans le cas où F size=1^猪 est identique à 0 (donc z = y ) et 猪 = 0, c’est-à-dire:

Cette équation équivaut (variation des constantes) à l’équation intégrale:

L’application P⁰, définie par intégration de (Q⁴) pendant le temps T, peut donc s’écrire:

Si 﨎 est assez petit, l’application P size=1^﨎 aura une forme analogue:

dans des coordonnées centrées sur son point fixe elliptique X^e . Si on oublie les termes d’ordre supérieur ou égal à quatre, on obtient la forme normale approchée 刺 size=1^﨎( 﨣) = 﨣e ^{i (b 0(} size=1^{﨎)+b 1(} size=1^﨎)| size=1^﨣|2) qui, comme le fait P⁰, laisse invariant chaque cercle centré en 﨣 = 0, et qui sur chacun de ces cercles induit une rotation dont l’angle varie (décroît ici) lorsque le rayon croît (propriété de torsion (twist ) typiquement non linéaire (fig. 38, que l’on comparera à la fig. 9).

La description de la dynamique de P size=1^﨎, obtenue en tenant compte des termes O (| 﨣|⁴) supprimés, est l’objet de théorèmes très remarquables, le théorème géométrique de Poincaré-Birkhoff et le théorème de Kolmogorov-Arnold-Moser (plus exactement ici le théorème des courbes invariantes ).

La résultat est (approximativement) décrit sur la figure 39, que nous allons commenter.

Tout d’abord, cette figure ne décrit qu’un «petit» voisinage 六 du point fixe elliptique. Dans ce voisinage, P size=1^﨎 laisse invariantes une infinité de courbes fermées proches de cercles centrés sur le point fixe mais, contrairement à ce qui se passe pour 刺 size=1^﨎, ces courbes fermées ne remplissent pas le voisinage 六: leur intersection avec un rayon issu du point fixe est homéomorphe à l’ensemble de Cantor déjà rencontré; cependant, elles sont de plus en plus denses au voisinage du point fixe (ce dernier est un point de densité, au sens où la mesure de Lebesgue relative de l’intersection avec un rayon tend vers 1 lorsque la longueur du rayon tend vers 0).

Les zones intermédiaires, nommées domaines d’instabilité par Birkhoff, sont bordées par des courbes invariantes lipschitziennes (théorie de Birkhoff ); elles contiennent des points périodiques elliptiques et hyperboliques dont la période est d’autant plus grande que l’on est proche du point fixe (un point périodique de période q est elliptique si c’est un point fixe elliptique de P size=1^﨎q ; on parle de même d’orbite périodique elliptique). Autour des points périodiques elliptiques de période q , la dynamique de P size=1^﨎q est analogue à celle de P size=1^﨎 au voisinage de son point fixe (cela explique en particulier que beaucoup de chercheurs essaient d’appliquer à ce problème les méthodes du groupe de renormalisation ). Quant aux points périodiques hyperboliques de période q , leurs variétés stables et instables (définies en les considérant comme points fixes de P size=1^﨎q ) s’intersectent transversalement suivant une figure complexe, qui est l’exact analogue de la figure 32.

Les théorèmes cités affirment en fait que cette description correspond à la dynamique, au voisinage d’un point fixe elliptique, d’un difféomorphisme du plan conservant les aires suffisamment général (on dit plutôt générique , cf. chap. 9).

L’importance de l’existence de courbes fermées invariantes qui l’entourent dans tout voisinage du point fixe tient à ce qu’elle implique immédiatement la stabilité de ce dernier (au sens de Liapunov dans le futur et le passé): l’orbite d’un point proche est en effet piégée par ces courbes qui séparent chacune l’espace en deux morceaux invariants par P size=1^﨎 (cela ne sera plus vrai en dimension supérieure).

La démonstration en est rendue difficile par la présence de petits dénominateurs qui font diverger les changements de coordonées formels (c’est-à-dire définis par des séries éventuellement divergentes) mettant P size=1^﨎 sous forme normale, même si P size=1^﨎 est analytique.

Par exemple, si b ⁰( 﨎)/2 神 n’est pas rationnel, il existe un changement de coordonnées formel conservant les aires qui transforme P size=1^﨎 en une application de la forme 﨣 料 﨣ei 﨏(| 﨣|²) laissant invariant chaque cercle de centre O (forme normale de Birkhoff ). Ces changements de coordonnées peuvent être obtenus en composant une infinité de changements de coordonnées analogues à celui qui, dans l’équation (Q ³) du chapitre 3, fait passer de 兀 à 﨣, c’est-à-dire:

leur divergence vient de ce que les coefficients 塚^kl sont de la forme:

et comportent donc des dénominateurs pouvant devenir arbitrairement petits quand (k 漣 l 漣 1) est assez grand; elle a donc pour cause la présence, dans tout voisinage de P size=1^﨎, d’applications ayant un point fixe résonnant (c’est-à-dire tel que la dérivée en ce point soit une rotation périodique) pour lesquelles certains de ces dénominateurs s’annulent.

A posteriori, cette divergence n’est pas pour nous étonner puisque dans ces coordonnées la dynamique de P size=1^﨎 apparaît comme beaucoup trop simple!

C’est à Kolmogorov (1954) que revient l’idée de localiser le problème au voisinage d’un cercle invariant de la forme normale tronquée 刺 size=1^﨎 et non plus au voisinage du point fixe elliptique lui-même. Plus précisément, Kolmogorov choisit un cercle sur lequel 刺 size=1^﨎 induit une rotation d’angle 2 神諸, où 諸 est mal approché par les rationnels, c’est-àdire: il existe des constantes C 礪 0 et 廓 閭 0 telles que, pour tout rationnel p /q , on ait:

ce choix est naturel dans la mesure où l’on s’attend à ce que les cercles sur lesquels 刺 size=1^﨎 induit une rotation d’angle 2 神 p/q , donc entièrement formés de points périodiques de période q , soient détruits par la perturbation qui fait passer de 刺 size=1^﨎 à P size=1^﨎. La deuxième idée est de chercher dans ce voisinage une courbe fermée invariante par P size=1^﨎 sur laquelle P size=1^﨎 induit un difféomorphisme conjugué à la rotation d’angle 2 神諸 (ce qui signifie qu’il existe un difféomorphisme de cette courbe sur un cercle, transformant la restriction de P size=1^﨎 en cette rotation). La troisième idée est d’utiliser à cette fin une méthode d’itération rapide (type méthode de Newton). Les premières démonstrations complètes de théorèmes de ce type ont été données par Arnold (1961) dans le cas analytique et Moser (1962) dans le cas différentiable. Une version très intéressante de ce théorème a été donnée par Rüssmann (1970): on n’exige plus que P size=1^﨎 préserve les aires et, au lieu de trouver des courbes invariantes, on trouve des courbes dont l’image (globale, pas celle de chaque point) est obtenue par un déplacement constant le long des rayons (courbes translatées ); ce théorème de Rüssmann joue un rôle important dans l’étude des bifurcations des difféomorphismes non conservatifs de R² au voisinage d’un point fixe elliptique; il montre que, dans ce problème, la conservation de l’aire joue le rôle d’une condition de codimension 1 (annuler la translation).

Si on compare au théorème de bifurcation de Hopf décrit au chapitre 3, ou plutôt à son analogue pour les familles de difféomorphismes, la difficulté est ici beaucoup plus grande car les cercles laissés invariants par la forme normale approchée n’ont a priori aucune robustesse vis-à-vis des perturbations (ils ne sont pas normalement hyperboliques). Cela se retrouve au niveau des démonstrations: théorème usuel des contractions pour Hopf, théorème raffiné de fonctions implicites dans un espace vectoriel métrique (non normé) complet (théorème de Nash-Moser sous la forme donnée par Sergeraert ou Hamilton qu’utilise Herman dans la démonstration – sans doute la plus transparente de toutes –qu’il a donnée du théorème de Rüssmann).

Lorsque le difféomorphisme n’est que peu différentiable, ces théorèmes deviennent faux; des contre-exemples ont été donnés par Takens (en classe C¹) et Herman (en classe C³ size=1^漣 size=1^﨎). Il est clair que la méthode ne fonctionne pas, puisqu’on ne peut même pas écrire de forme normale approchée. Enfin, même dans le cas indéfiniment différentiable C 秊, il existe des difféomorphismes conservant les aires qui possèdent une orbite dense dans un disque centré sur un point fixe elliptique: une telle orbite passe au voisinage de chaque point du disque et sa présence exclut la possibilité d’une courbe fermée invariante entourant le point fixe. Ces exemples non génériques sont dus à Anosov et Katok; ils sont rendus possibles par le viol de la propriété de torsion: leur forme normale formelle est linéaire.

Le destin des cercles invariants de 刺 size=1^﨎 sur lesquels 刺 size=1^﨎 induit une rotation d’angle 2 神 p/q est d’être en général transformés en un nombre fini d’orbites périodiques de période q de P size=1^﨎.

Cela est vrai pour tout p/q 麗 b ⁰( 﨎)/2 神 et assez proche de b ⁰( 﨎)/2 神 mais, si q n’est pas trop grand par rapport à:

la démonstration en est locale et aisée alors que, si q est trop grand, elle devient globale et d’une grande subtilité. Cela vient de ce que P^q ne peut être considéré comme une perturbation de 刺 size=1^﨎q au voisinage du cercle incriminé que si q n’est pas trop grand par rapport au rayon de ce cercle. Dans ce cas, il est facile de prouver à l’aide du théorème des fonctions implicites usuel l’existence d’une courbe fermée proche d’un cercle transformée «radialement» par P size=1^﨎q (fig. 40); la conservation de l’aire implique qu’une telle courbe recoupe son image: les points d’intersection sont des points fixes de P size=1^﨎q , c’est-à-dire des points périodiques de période q de P size=1^﨎.

Zehnder a pu montrer qu’en général on obtient ainsi deux orbites périodiques, l’une elliptique, l’autre hyperbolique. En considérant ces orbites périodiques comme des «souvenirs de résonances», c’est-à-dire comme obtenues par perturbation à partir d’un point fixe elliptique d’un difféomorphisme proche dont la dérivée en ce point est 﨣 料 e ^{i 2} size=1^神p/q 﨣, il a pu montrer également qu’il y a «en général» intersection transversale des variétés stables et instables des orbites périodiques hyperboliques comme sur la figure 39.

Si q est grand, on ne peut prouver l’existence de ces orbites périodiques qu’en utilisant des arguments plus globaux. Considérons la région située entre deux courbes invariantes correspondant à des angles 2 神諸¹ et 2 神諸² tels que 諸¹ 麗 p /q 麗 諸²; en itérant q fois P size=1^﨎, on obtient un difféomorphisme de cet anneau qui conserve les aires et qui «fait tourner les bords» dans deux sens différents (fig. 41, à comparer avec la fig. 33 bis ; cela n’a vraiment de sens que si on déroule l’anneau...).

Par une interprétation très subtile de la conservation de l’aire (existence de «presque-orbites» traversant l’anneau d’un bord à l’autre) et un argument plus classique de variation d’angle, Birkhoff a démontré en 1912 que, comme l’avait conjecturé Poincaré la même année, juste avant sa mort, un tel homéomorphisme a au moins deux points fixes (la démonstration de l’existence du deuxième date en fait de 1925). Un tel théorème semblait jusqu’ici échapper aux méthodes habituelles de la topologie algébrique, mais Conley et Zehnder viennent (déc. 1982) d’en donner une démonstration très simple fondée sur la très classique approche variationnelle, qui a le mérite immense de se généraliser aux dimensions supérieures (cf. infra , Excursion variationnelle ). Notons que de tels difféomorphismes de l’anneau se présentent naturellement dans le fameux problème du billard , étudié par Birkhoff, auquel toute cette théorie s’applique (cf. chap. 8).

Revenons maintenant à la taille du voisinage dans lequel l’analyse qui précède est valable, problème important pour les applications (stabilisation de plasmas, dit-on, etc.) qui est actuellement l’objet de nombreuses recherches fondées sur la technique du «groupe de renormalisation». Si les démonstrations initiales donnaient un voisinage beaucoup plus petit que ne le faisaient les simulations numériques (par exemple celles de Henon et Heiles sur le mouvement d’une étoile dans une galaxie cylindrique), une nouvelle démonstration de Herman fondée sur le théorème de point fixe de Leray-Shauder donne, dans le cas d’une classe de nombres 諸 spécialement mal approchés par les rationnels, des estimations plus raisonnables. C’est d’autant plus important qu’il est bien difficile de distinguer numériquement une courbe d’un ensemble de Cantor à petits trous. Cette remarque n’est pas gratuite dans la mesure où Aubry et Mather ont montré récemment que les cercles invariants de 刺 size=1^﨎 sur lesquels 刺 size=1^﨎 induit une rotation d’angle 2 神諸 se transforment, lorsque 諸 est trop bien approché par les rationnels, en ensembles de Cantor invariants par P size=1^諸, intermédiaires entre des orbites périodiques et des courbes fermées invariantes.

Nous finirons ce chapitre par quelques commentaires hamiltoniens, volontairement restreints, puisque nous avons choisi de ne pas développer ici tout ce formalisme.

Soit H(x ¹, ..., x ⁿ , y ¹, ..., y ⁿ ) une fonction différentiable définie sur R²ⁿ . Le système dynamique hamiltonien associé à H s’écrit (comparer à l’écriture de (P¹) au chap. 1):

et le flot associé laisse évidemment invariante chaque hypersurface de niveau ^c de H (hypersurfaces d’énergie constante H = c ).

Si n est égal à 2 (systèmes à deux degrés de liberté, par exemple le problème restreint des trois corps ), ces hypersurfaces sont de dimension trois. L’étude de la restriction de (H) à une hypersurface ^c au voisinage d’une orbite périodique se ramène à l’étude précédente en considérant l’application de Poincaré P définie par le flot sur un petit morceau de surface S transverse à l’orbite périodique (fig. 42). Celle-ci conserve les aires et la seule différence avec les exemples précédents est que le temps de retour sur le plan n’est plus constant: la dynamique de P ne renseigne plus que sur la géométrie du portrait de phase de (H) dans ^c et non sur les paramétrages des courbes intégrales.

Cependant, sur un tore invariant de (H) correspondant à une courbe invariante de P sur laquelle P est conjuguée à une rotation d’angle 2 神諸, Kolmogorov a montré que, si 諸 est mal approché par les rationnels, les courbes intégrales sont quasi périodiques : cela signifie qu’il existe un système de coordonnées sur le tore (c’est-à-dire un système de coordonnées périodiques en les deux variables sur R²) dans lequel la restriction au tore du champ de vecteurs 輻^H défini par H est un champ constant de pente irrationnelle. La démonstration, qui revient à rendre constant le temps de premier retour sur S, se ramène à la résolution d’une équation aux différences :

(on cherche 﨡 Z-périodique connaissant 﨏 Z-périodique telle que 咽10 﨏( ) d = 0).

La résolution formelle par séries de Fourier montre immédiatement qu’il y a un problème de petits dénominateurs, ici très simple.

La figure 42 représente une partie d’une telle courbe intégrale dans une hypersurface d’énergie constante (cette courbe est dense sur le tore). Comme précédemment, l’existence de tels tores dans tout voisinage de l’orbite périodique prouve sa stabilité au sens de Liapunov.

Si maintenant on fait varier l’énergie, on est amenés à considérer 輻^H comme une famille à un paramètre ( 輻^H| ^c )^c de champs de vecteurs. Les bifurcations qui se produisent sont très riches (bifurcation d’orbites périodiques à partir de résonances, etc.); étudiées par Deprit et Henrard, Meyer, Palmore, Abraham, Takens, elles sont joliment dessinées dans l’ouvrage Foundations of Mechanics d’Abraham et Marsden.

Lorsque le nombre de degrés de liberté est supérieur à deux, les choses se compliquent: les théorèmes de K.A.M. fournissent des tores invariants Tⁿ pour des systèmes voisins de systèmes complètement intégrables, mais de tels tores ne séparent plus en deux composantes les surfaces d’énergie constante, qui sont de dimension 2n 漣 1, et leur existence, si elle suffit à ruiner l’hypothèse ergodique de Boltzmann, ne suffit plus à assurer des résultats de stabilité: c’est le phénomène de diffusion d’Arnold dont des illustrations numériques ont été obtenues par Froeschle en couplant légèrement deux difféomorphismes de R² ayant un point fixe elliptique. En conclusion, notre connaissance du portrait de phase d’équations telles que (N⁴) ou (N ⁴) est loin d’être complète. Ce long chapitre a donné le premier rôle aux orbites périodiques (et quasi périodiques) et c’est d’elles encore dont il sera question dans l’excursion variationnelle qui suit. Dans l’étude des points homoclines, nous avons vu que les points périodiques de P size=1^﨎 étaient denses dans l’ensemble de Cantor invariant associé à un tel point; dans la situation de la figure 39, il n’est pas déraisonnable de conjecturer avec Poincaré que les points périodiques sont denses au voisinage du point fixe, c’est-à-dire que les orbites périodiques (de longues périodes) sont denses au voisinage d’une orbite périodique elliptique d’un système hamiltonien à deux degrés de liberté.

L’importance de ces orbites pour la description des mouvements généraux est donc très grande. Finissons comme nous avons commencé, en citant Poincaré (Méthodes nouvelles de la mécanique céleste , t. I, chap. III, p. 82): «D’ailleurs, ce qui nous rend ces solutions périodiques si précieuses, c’est qu’elles sont pour ainsi dire la seule brêche par où nous puissions essayer de pénétrer dans une place jusqu’ici réputée inabordable.»

Excursion variationnelle

On sait, depuis les travaux de Maupertuis, Euler, Lagrange, Hamilton, Jacobi, que le principe de moindre action est l’un des fondements possibles de la mécanique, les équations du mouvement apparaissant comme les équations d’Euler associées à l’intégrale d’action (cf. CALCUL DES VARIATIONS et MÉCANIQUE ANALYTIQUE, ainsi que Arnold, Méthodes mathématiques de la mécanique classique ); cela reste vrai pour des équations hamiltoniennes du type (H) (cf. supra ) et ce même si H dépend du temps.

En particulier, une orbite périodique de période T de (H) n’est autre qu’un point critique de la fonctionnelle:

définie sur un espace de fonctions 塚: [0, T]R²ⁿ suffisamment régulières (par exemple l’espace de Sobolev H¹) telles que 塚(0) = 塚(T). Le lecteur pourra formellement vérifier que l’accroissement 嗀 崙( 塚) = 崙( 塚 + 嗀塚) 漣 崙( 塚) correspondant à un accroissement 嗀塚 = ( 嗀x ¹, ..., 嗀x ⁿ , 嗀y ¹, ..., 嗀y ⁿ ) de 塚 se met sous la forme:

Cette approche a donné de nombreux résultats, en particulier lorsque les surfaces d’énergie constante ont de bonnes propriétés de convexité (travaux de Weinstein, Moser, Rabinowitz, Clarke, Ekeland, Lasry...), mais les orbites périodiques relevant de l’application du théorème géométrique de Poincaré-Birkhoff dont il a été question dans ce chapitre semblaient lui échapper. On peut s’en étonner aujourd’hui au vu de la grande simplicité de la démonstration que Conley et Zehnder viennent de donner (déc. 1982) du théorème suivant, conjecturé du bout des lèvres par Arnold, qui implique le théorème de Poincaré-Birkhoff: Soit f : T²T² un difféomorphisme du tore ayant un relèvement à R² de la forme face="EU Caron" キ (x , y ) = (x + 﨏(x , y ), y + 祥(x , y )), où 﨏 et 祥 sont périodiques de période 1 en x et y . Si f préserve les aires et le centre de gravité (c’est-à-dire:

ce qui évite toute «rotation» systématique, cf. infra , Excursion homologique , in chap. 5), f possède au moins trois points fixes (et en général au moins quatre). Les hypothèses permettent d’écrire f comme l’application de premier retour (au bout d’une période) d’une équation différentielle du type (H) dont le hamiltonien est périodique en t de période 1. La méthode variationnelle s’applique sans problème à la fonctionnelle 崙 (ou plutôt à celle, analogue mais symétrique, obtenue en remplaçant y dx par (1/2) (y dx 漣 x dy )); un argument classique du type Liapunov-Schmidt (troncature des séries de Fourier) ramène le problème en dimension finie, et un argument de topologie algébrique permet de conclure à l’existence d’au moins trois points critiques (le nombre 3 est le nombre minimal de points critiques que possède une fonction différentiable sur T²). L’intérêt de la méthode est qu’elle se généralise aux difféomorphismes symplectiques du tore T²ⁿ qui conservent le centre de gravité: ceux-ci doivent posséder au moins (n + 1) points fixes (et en général au moins 2ⁿ ); l’inconvénient est qu’elle nécessite des hypothèses rigides de type métrique alors que des hypothèses topologiques plus souples seraient peut-être (?) suffisantes.

Un autre succès récent de l’approche variationnelle fut la démonstration par Aubry, puis Mather, de l’existence dans les zones d’instabilité d’ensembles de Cantor invariants de nombre de rotation donné.

Disons simplement que la méthode généralise la remarque suivante, due à Birkhoff: parmi tous les polygones à q côtés inscrits dans un convexe du plan, et ayant un ordre fixé de leurs sommets, ceux dont la somme des longueurs des côtés est maximale (resp. minimax) définissent des trajectoires périodiques de période q du jeu de billard correspondant (pour la définition, cf. chap. 8).

5. Perturbations périodiques d’auto-oscillations et «courbes de Birkhoff»

Nous reprenons l’étude du chapitre précédent en présence d’un frottement ou d’une excitation s’annulant avec la vitesse. Plus précisément, nous considérons des équations du type suivant, généralisant (N ⁴):

qu’on écrira:

ou:

Nous envisagerons tout d’abord le cas où |f (x )| et |g (t )| sont petits (plus précisement les petites perturbations périodiques des équations (N²) et (N³)), occasion de discuter la fonction de Melnikov, la bifurcation de Hopf des difféomorphismes, et la belle théorie des difféomorphismes du cercle, puis le cas où |f (x )| et |g (t )| sont très grands, qui se réduit pratiquement à l’équation de Van der Pol forcée étudiée par Cartwright-Littlewood, Levinson, Levi, exemple dont l’importance historique et paradigmatique est considérable.

L’équation:

avec par exemple p (t ) = cos 2 神t T, | 﨎| petit, peut être considérée, lorsque | 福| est petit, comme une perturbation de l’équation (N ⁴).

Si 福 = 0, les variétés stables et instables du point fixe hyperbolique X^h 捻 S¹ 憐 R de l’application de Poincaré P size=1^﨎 = P^0, size=1^﨎 ont une intersection transverse, donc insensible aux petites perturbations.

Il en est donc de même tant que | 福| est assez petit (c’est-à-dire qu’il existe un point fixe hyperbolique de l’application de Poincaré P size=1^福, size=1^﨎 associée à (N⁵) et ses variétés stables et instables se coupent transversalement). Cependant, si 﨎 = 0, on déduit du chapitre 2 qu’il n’y a plus d’intersection dès que 福 0 (fig. 13, 14, 15 si 福 礪 0).

Afin d’évaluer la région du plan ( 福, 﨎) correspondant à des paramètres pour lesquels l’application de Poincaré P size=1^福, size=1^﨎 associée à (N ⁵) possède des points homoclines, nous introduisons la fonction de Melnikov , de portée assez générale, mais surtout utile dans l’étude des petites perturbations de systèmes conservatifs. Nous raisonnerons sur la figure 43, faite dans l’espace de phase (S¹ 憐 R) 憐 S^1T (déroulé en R² 憐 [0, T]) de l’équation (P⁵).

Dans le cas où 福 = 﨎 = 0, cette figure devient une partie du portrait de phase de l’équation:

(fig. 44).

Au point fixe hyperbolique X^h de l’application de Poincaré P size=1^福, size=1^﨎 correspond une orbite périodique «hyperbolique» 六^h de l’équation différentielle; aux variétés stables et instables de X^h correspondent les variétés stables et instables W^s (face=F0021 六^h ) et W^u (face=F0021 六^h ) de 六^h , définies de façon analogue (ensemble des courbes intégrales positivement ou négativement asymptotiques à 六^h ). Lorsque 福 = 﨎 = 0, ces variétés invariantes se confondent; la remarque évidente qui est à la base de la méthode de Melnikov est la transversalité de l’unique intersection des courbes intégrales situées sur ces variétés avec le plan M d’équation x = 0 (fig. 44). Cette transversalité est encore réalisée pour 福, 﨎 assez petits et il existe donc pour tout une unique courbe intégrale t塚^s (t , ) contenue dans W^s (face=F0021 六^h ) telle que 塚^s ( , ) 捻 M; de même, il existe une unique courbe intégrale t塚^u (t , ) contenue dans W^u (face=F0021 六^h ) telle que 塚^u ( , ) 捻 M (on paramètre les courbes intégrales de façon que 精(t ) = t ). Les deux points 塚^u ( , ) et 塚^s ( , ) sont situés sur la droite «verticale» intersection de M et de size=1 =(x , y , 精)| 精 = . On note ( ) la distance 塚^s ( , ) 漣 塚^u ( , ) mesurée suivant cette droite (parallèle à l’axe des y ; fig. 43). Il n’est pas dificile d’effectuer un calcul approché de ( ): Soit t(x ^h (t ), y^h (t )) la paramétrisation de la variété instable (stable) de l’équilibre hyperbolique de l’équation (P¹) (pendule sans frottement) telle que x ^h (0) = 0; on note 嗀(t , ) l’aire algébrique du parallélogramme de R² 令 ^t engendré par les vecteurs 塚^s (t , ) 漣 塚^u (t , ) et (y ^h (t 漣 ), face=F0019 漣 諸²sinx ^h (t 漣 )) (fig. 45). Lorsque t = , ces deux vecteurs sont orthogonaux, et 嗀( , ) = y ^h (0) 練 ( ) = 2 諸( ) (fig. 46).

En dérivant 嗀(t , ) par rapport au temps et en tenant compte des comportements asymptotiques de 塚^s (t , ) et 塚^u (t , ), on montre que:

Il est pratiquement clair sur la figure 43 qu’il existe ⁰ tel que ( ⁰) = 0 si et seulement si W^s (face=F0021 六^h ) et W^u (face=F0021 六^h ) se rencontrent; on montre également que l’intersection est transversale si et seulement si ( ) s’annule transversalement en ⁰ (c’est-à-dire, avec une dérivée non nulle). On obtient alors facilement la figure 47 qui décrit, suivant les valeurs de ( 福, 﨎) assez proches de (0, 0), le comportement des variétés stables et instables du point fixe hyperbolique de l’application de Poincaré P size=1^福, size=1^﨎 (ou plus exactement d’une moitié de celles-ci). Comme d’habitude, on a déroulé S¹ 憐 R sur R².

Ce type d’étude est possible chaque fois qu’on connaît explicitement la solution non perturbée (x ^h (t ), y ^h (t )). De nombreux calculs, faisant intervenir des fonctions spéciales, ont été faits dans des cas concrets, notamment par Holmes. Par la même méthode, on peut décrire les bifurcations des orbites périodiques sous-harmoniques de larges périodes, voisines de ces variétés stables et instables (voir les travaux de Chow, Hale, Mallet-Paret). L’existence de ces orbites périodiques est beaucoup plus facile à mettre en évidence que celle des orbites périodiques de larges périodes voisines d’un point fixe elliptique (cf. chap. 4), au moins pour celles qui ne passent que deux fois au voisinage du point fixe hyperbolique. Les courbes C¹ et C² de la figure 47 correspondent aux cas où les variétés stables et instables du point fixe hyperbolique de P size=1^福, size=1^﨎 sont tangentes. Les travaux de Newhouse montrent qu’au voisinage d’une telle situation il doit «génériquement» exister une infinité d’orbites périodiques asymptotiquement stables (ou asymptotiquement stables pour P^-1 size=1^福, size=1^﨎), ce qui donne de beaux exemples de systèmes dynamiques ayant de manière stable une infinité d’attracteurs.

Pour étudier P size=1^福, size=1^﨎 au voisinage de (0,0) lorsque 福 0, il est commode de la considérer comme une perturbation de P size=1^福,0 (qui n’est autre que 﨏^T, où ( 﨏^t )^t size=1^{捻 R} est le flot de l’équation (P²)). Le phénomène de résonance paramétrique ne se produit pas ici puisque, si 福 0, le point (0, 0) est un point fixe hyperbolique de P size=1^福,0 (le spectre de DP size=1^福,0(0, 0) est strictement à l’intérieur du disque unité si 福 礪 0 et strictement à l’extérieur si 福 麗 0). L’analogue pour 福 礪 0 de la figure 30 est représenté sur la figure 48: les zones hachurées correspondent aux cas où le point fixe de P size=1^福, size=1^﨎 proche de (0, 0) cesse d’être un attracteur.

Si par contre on remplace le frottement linéaire 福dx/dt par le terme non linéaire f size=1^猪(x ) dx/dt introduit au chapitre 3, une phénoménologie très riche apparaît, due en particulier à la possibilité de résonances entre la période T des oscillations forcées et la période propre, voisine de 2 神/ 諸, de l’auto-oscillation qui naît par «bifurcation de Hopf» pour 猪 礪 0 lorsque 﨎 = 0. Autrement dit, pour retrouver ce phénomène de résonance paramétrique, on compense la non-conservativité du système par la présence d’un paramètre faisant passer l’équilibre considéré de l’attraction à la répulsion. Nous ne décrirons pas en détail ces phénomènes, nous contentant du «cas général» où T ne résonne «pas trop» avec 2 神/ 諸; c’est une nouvelle occasion de rencontrer la bifurcation de Hopf, cette fois pour des difféomorphismes (apparition de courbes fermées invariantes attractantes par déstabilisation d’un point fixe dans une famille de difféomorphismes ou, ce qui revient au même, apparition de tores invariants attractants par déstabilisation d’une orbite périodique dans une famille d’équations différentielles).

Nous noterons P size=1^猪, size=1^﨎 l’application de Poincaré de l’équation:

Comme au chapitre 4, celle-ci peut être évaluée à l’aide de la formulation intégrale de (N⁵) (méthode de variation des constantes).

En particulier, dans les coordonnées locales 﨣 du chapitre 3, l’intégration pendant le temps T de (N³) fournit:

que l’on comparera à la formule (R⁴) du chapitre 4.

Pour 猪 礪 0, l’orbite périodique de (N³) apparue par bifurcation de Hopf (cf. chap. 3) est laissée invariante par P size=1^猪,0; on retrouve sur (R⁵) que cette courbe est proche du cercle | 﨣| = 連3 猪.

Lorsque T n’est pas un multiple entier de 2 神/ 諸, le théorème des fonctions implicites assure l’existence pour | 﨎| petit d’un unique point fixe X size=1^猪, size=1^﨎 de P size=1^猪, size=1^﨎, proche de (0, 0) et dépendant différentiablement de 﨎 et 猪.

La distance du spectre size=1^猪, size=1^﨎, size=1^猪, size=1^﨎 de DP size=1^猪, size=1^﨎(X size=1^猪, size=1^﨎) à celui de DP size=1^猪,0(0, 0) est de l’ordre de 﨎; la figure 49 représente ces spectres dans le cas général: on remarquera que les deux paramètres réels 猪, 﨎 peuvent être remplacés par le paramètre complexe size=1^猪, size=1^﨎 (on suppose toujours | 猪| et | 﨎| assez petits).

Si q T n’est pas un multiple entier de 2 神/ 諸 pour q = 2, 3 ou 4 (c’est-à-dire^{q 0,0} 1, q = 2, 3 ou 4), il existe un système de coordonnées (encore appelé 﨣) centré sur X size=1^猪, size=1^﨎 et dépendant différentiablement de 猪, 﨎 dans lequel P size=1^猪, size=1^﨎 s’écrit, à des termes d’ordre 4 près, comme une forme normale approchée :

et, bien entendu,

La démonstration est la même que dans l’établissement des formes normales approchées du chapitre 3; l’exclusion des résonances fortes^{q 0,0} = 1, q = 2, 3 ou 4, vient de ce qu’une relation de la forme^{q 0,0} = 1 rend impossible l’élimination des termes ( 﨣)^{q -1} dans le développement de Taylor de P size=1^猪, size=1^﨎 en X size=1^猪, size=1^﨎 par un changement de variables dépendant continûment de 猪, 﨎 au point 猪 = 﨎 = 0 (présence dans le changement de variables de termes en:

qui «explosent» pour 猪 = 﨎 = 0 et ce même si, dans notre cas, le développement de Taylor de P^0,0( 﨣) ne comporte pas lui-même de terme ( 﨣)^{q -1}).

Lorsque q = 2, 3 ou 4, ces termes sont de degré inférieur ou égal à 3 et concurrencent le terme en 﨣| 﨣|² dans la détermination de la dynamique des formes normales approchées au voisinage de 﨣 = 0.

L’étude de (R⁵), fondée sur la méthode des transformées de graphes , est en tout point parallèle à l’étude de (Q ³) faite au chapitre 3: lorsque a ⁰( 猪, 﨎)/a ¹( 猪, 﨎) est négatif, c’est-à-dire a ⁰( 猪, 﨎) 礪 0, P size=1^猪, size=1^﨎 laisse invariante une courbe fermée 暈 猪, 﨎 entourant X size=1^猪, size=1^﨎, proche dans les coordonnées 﨣 du cercle:

et attirant tout point 﨣 0 voisin de 0 (c’està-dire lim [distance (Pⁿ size=1^猪, size=1^﨎( 﨣), face=F0021 暈 猪, 﨎)] = 0). Le n 轢+ 秊passage de( 猪, 﨎) de l’intérieur à l’extérieur du cercle unité de C correspond donc au remplacement d’un attracteur ponctuel de P size=1^猪, size=1^﨎 par un attracteur homéomorphe à un cercle. Ce théorème de bifurcation de Hopf des difféomorphismes, bien que fortement analogue à son cousin du chapitre 3, en diffère assez profondément; tout d’abord, ce n’est pas un théorème de nature topologique ne concernant que le spectre: si dans une famille F size=1^益 à un paramètre de difféomorphismes locaux de R² fixant 0, le spectre de DF size=1^益(0) traverse le cercle unité de C lorsque 益 varie, il se peut qu’aucun des F size=1^益 ne possède de courbe fermée invariante entourant 0, et ce même s’il n’y a pas de résonance. L’existence de telles courbes dépend de conditions non linéaires comme l’«attractivité» du point fixe 0 de P^0,0. D’autre part, le passage d’un difféomorphisme à sa forme normale approchée simplifie exagérément la dynamique sur les courbes invariantes, remplaçant en particulier par des rotations périodiques les difféomorphismes de ces courbes ayant des orbites périodiques isolées et masquant ainsi le phénomène connu des physiciens sous le nom d’accrochage de fréquences (qui implique également une perte de régularité de la courbe invariante lorsque 猪 croît, suggérée par le passage de la figure 15 à la figure 13).

La figure 50 indique les zones de l’espace des paramètres (T, 猪) où se produit ce phénomène lorsqu’on fixe une valeur de 﨎; il faut bien entendu que 﨎 soit non nul, car P size=1^猪,0 faisant partie du groupe à un paramètre associé à l’équation (N³), sa restriction à une courbe fermée invariante est toujours différentiablement conjuguée à une rotation et n’a donc pas d’orbite périodique isolée. Sur cette figure 50, C désigne (dans le cas 﨎 = 0) la courbe formée des couples (T, 猪) tels que P size=1^猪,0 ait une courbe fermée invariante 暈 猪,0 sur laquelle il induit un difféomorphisme conjugué à la rotation d’angle 2p 神/q 諸; L désigne (pour 﨎 0) la languette formée des couples (T, 猪) tels que P size=1^猪, size=1^﨎 ait une courbe fermée invariante 暈 猪, 﨎 sur laquelle il induit un difféomorphisme ayant un nombre pair d’orbites périodiques isolées de période q (faisant p tours avant de se refermer), une moitié étant attractante et l’autre répulsive (ces orbites s’éliminent par paires sur le bord de la languette).

La figure 51 indique la dynamique de la restriction, notée f , de P size=1^猪, size=1^﨎 à 暈 猪, 﨎 lorsque (T, 猪) est dans la languette représentée sur la figure 50 (on a choisi p = 2, q = 5).

On voit que, lorsque 猪 décroît, T = ⁰ et 﨎 0 étant fixés, la restriction de P size=1^猪, size=1^﨎 à sa courbe fermée 暈 猪, 﨎 passe par une infinité de bifurcations consistant en l’apparition et la disparition de paires d’orbites périodiques dont la période tend vers l’infini quand 猪 se rapproche de la valeur pour laquelle se produit la bifurcation de Hopf. Sur la famille d’équations différentielles autonomes () dans R² 憐 S^1T associée à (N⁵), cela correspond à l’apparition et à la disparition d’orbites périodiques de période q T situées sur un tore invariant et le parcourant q fois dans une direction (celle de S^1T) et p fois dans l’autre. En général, contrairement à ce qui se passe pour 﨎 = 0, chacune de ces paires d’orbites périodiques existe stablement sur tout un intervalle de valeurs de 猪: c’est de là que vient la dénomination accrochage de fréquences .

Le phénomène que nous venons de décrire est aujourd’hui bien compris: son analyse repose sur le nombre de rotation d’un homéomorphisme du cercle, introduit par Poincaré en 1885. L’importance de cette notion et la beauté de la théorie des difféomorphismes du cercle qu’elle a suscitée (théorèmes de Denjoy, Arnold, Herman) justifient qu’on s’y arrête quelque peu.

Un homéomorphisme f du cercle de rayon 1 de R² dans lui-même conservant l’orientation peut toujours se relever (se «dérouler») en un homéomorphisme face="EU Caron" キ de R dans R de la forme face="EU Caron" キ ( ) = + 﨏( ), où 﨏 est une fonction 2 神-périodique (c’est-à-dire 﨏( + 2 神) = 﨏( ) pour tout ); f est alors défini en interprétant comme un angle, c’est-à-dire f ( modulo 2 神) = face="EU Caron" キ ( ) modulo 2 神. Poincaré définit le nombre de rotation 福( face="EU Caron" キ ) de face="EU Caron" キ comme la limite (dont il montre qu’elle est uniforme et indépendante de ):

On se convainc facilement de ce que la classe 福(f ) de 福( face="EU Caron" キ ) modulo 2 神 est indépendante du choix du relèvement face="EU Caron" キ et représente une rotation moyenne de f (regarder l’exemple d’une rotation: les face="EU Caron" キ sont des translations face="EU Caron" キ ( ) = + 見 + k , k 捻 2 神Z).

Les propriétés de 福 sont les suivantes:

1. 福 ne dépend que de la classe de conjugaison de f dans le groupe des homéomorphismes du cercle (c’est-à-dire 福(f ) = 福(h ^-1 獵 f 獵 h )).

2. 福 dépend continûment de f (uniformité de la limite).

3. Théorème de Poincaré : 福(f ) 捻 2 神Q/2 神Z (c’est-à-dire 福(f )/2 神 est rationnel modulo 1) si et seulement si f a un point périodique. Si 福(f ) = 2 神 p /q modulo 2 神 avec p et q premiers entre eux, tous les points périodiques de f ont pour période q .

4. Théorème de Denjoy (1932 et 1946): si f est un difféomorphisme de classe C² et si 福(f ) 殮 2 神Q/2 神Z, f est conjugué par un homéomorphisme à la rotation d’angle 福(f ).

5. Théorème d’Herman (1976): il existe un sous-ensemble A de R 漣 Q dont le complémentaire a une mesure de Lebesgue égale à 0 tel que, si f est de classe C size=1^秊 et si 福(f ) 捻 2 神A/2 神Z, l’homéomorphisme conjuguant f à la rotation d’angle 福(f ) soit de classe C size=1^秊.

Poincaré eût été surpris par ces deux derniers théorèmes: par analogie avec la théorie des groupes fuchsiens, il croyait possible l’existence de difféomorphismes analytiques f du cercle vérifiant 福(f ) 殮 2 神Q/2 神Z et laissant invariant un ensemble de Cantor dans le cercle (et donc non topologiquement conjugués à une rotation). De tels difféomorphismes existent en classe C¹ (contre-exemples de Denjoy obtenus en «gonflant» une orbite d’une rotation non périodique) mais pas en classe C^k , k 閭 2. Quant au théorème d’Herman, remarquable contrepoint entre les propriétés arithmétiques du nombre de rotation et les propriétés de régularité de l’homéomorphisme conjuguant f à une rotation, il est la version globale (conjecturée par Arnold) d’un théorème antérieur d’Arnold qui montre le résultat lorsque f est analytique et suffisamment proche de la rotation d’angle 福(f ) en utilisant la méthode d’itération rapide citée au chapitre 4.

La figure 50 apparaît maintenant comme la division du plan des paramètres (T, 猪) suivant la valeur de 福(P size=1^猪, size=1^﨎| size=1^暈 size=1^猪, size=1^﨎). Si T = ⁰ et 﨎 sont fixés, l’application 猪 料 福(P size=1^猪, size=1^﨎| size=1^暈 size=1^猪, size=1^﨎) est non croissante au voisinage de la valeur de 猪 pour laquelle apparaît 暈 猪, 﨎.

Si 﨎 = 0, cette application est différentiable; lorsque 﨎 0, au contraire, cette application est localement constante sur une infinité d’intervalles (dont la réunion des intérieurs est localement dense) où elle prend des valeurs de la forme 2 神 p /q (modulo 2 神Z). Son graphe, analogue à celui de la fonction singulière de Lebesgue, est un escalier du diable (fig. 52, pour T = ⁰ et 﨎 fixés).

Ce comportement est générique pour des familles à un paramètre de difféomorphismes du cercle: il existe en général un ouvert dense de valeurs du paramètre pour lesquelles le difféomorphisme possède un nombre fini d’orbites périodiques attractantes et un nombre égal d’orbites périodiques répulsives comme sur la figure 51, où l’on a noté f = P size=1^猪, size=1^﨎| size=1^暈 size=1^猪, size=1^﨎. De tels difféomorphismes sont structurellement stables (une petite perturbation conduit à un difféomorphisme topologiquement conjugué), et l’affirmation ci-dessus rentre dans le cadre du théorème de stabilité de Peixoto que nous évoquerons au chapitre 9.

Un complément apporté par Herman à son théorème (et déjà observé par Arnold dans le cas local) nous permet de préciser ce qui précède. En étudiant les propriétés de l’application f 料 福(f ) en un point où f est différentiablement conjuguée à une rotation, il montre que dans toute famille à un paramètre (C¹ par rapport au paramètre) de difféomorphismes f size=1^猪 de classe C³ du cercle où le nombre de rotations 福(f size=1^猪) varie effectivement, l’ensemble des valeurs du paramètre 猪 pour lesquelles f size=1^猪 est différentiablement conjuguée à une rotation non périodique est de mesure de Lebesgue positive. Ce résultat s’applique dans notre cas: les intervalles de valeurs de 猪 sur lesquels 福(P size=1^猪, size=1^﨎| size=1^暈 size=1^猪, size=1^﨎) est constant sont denses au voisinage de la bifurcation de Hopf, mais leur complémentaire est un ensemble de Cantor de mesure de Lebesgue positive.

Les courbes invariantes de P size=1^猪, size=1^﨎 que nous venons de voir apparaître par bifurcation de Hopf ne persistent, quand on continue à faire varier les paramètres, que si leur attraction reste appréciable (hyperbolicité normale), et les modalités de leur disparition sont loin d’être bien comprises, sauf bien sûr si cette dernière correspond à une bifurcation de Hopf inversée.

Ce problème n’est pas étranger à celui, évoqué au chapitre 4, de la taille du domaine de validité du théorème de K.A.M., ce qui montre sa difficulté.

Le cas où la disparition coïncide avec l’existence sur la courbe d’un nombre fini d’orbites périodiques a fait l’objet d’études numériques de Curry et Yorke (1978) et de travaux plus théoriques de Aronson, Chory, Mac Gehee, Hall (1982). Lorsque l’attraction devient insuffisante, une dynamique complexe peut apparaître avec, en particulier, création de points homoclines: le difféomorphisme s’aperçoit en quelque sorte que les variétés stables et instables de points périodiques hyperboliques consécutifs n’ont plus aucune raison de coïncider pour former une portion de courbe différentiable, et peuvent au contraire s’intersecter transversalement; si simultanément les orbites périodiques coalescent par paires et sont sur le point de disparaître, la courbe peut perdre sa différentiabilité par une suite infinie de repliements s’accumulant au voisinage des points périodiques. Newhouse, Palis, et Takens ont montré qu’après qu’elle a disparu, un itéré de f (d’autant plus grand qu’on est plus proche du point de disparition) possède un fer à cheval. Des attracteurs étranges «à la Henon» (cf. chap. 6) semblent également surgir: c’est l’un des modèles de ces routes vers la turbulence qui sont depuis une dizaine d’années l’objet d’un riche dialogue entre mathématiciens, numériciens, physiciens théoriciens et physiciens expérimentateurs (on en reparlera à l’occasion dans les chapitres suivants).

La diminution de l’attraction qu’exerce une courbe invariante peut être liée à la proximité d’une courbe invariante répulsive, situation qui dans le cas plus simple de points fixes d’un endomorphisme de l’intervalle a été proposée par Pommeau et Manneville comme modèle de l’intermittence (bouffées de turbulence). Dans le cas des flots de R², la situation générale est également très simple, puisqu’elle se ramène à l’étude des points fixes d’une application de Poincaré: une orbite périodique attractante et une orbite périodique répulsive se rapprochent, se confondent (l’unique orbite survivante est alors répulsive d’un côté, attractante de l’autre), et disparaissent. Dans le cas des difféomorphismes de R² que nous considérons ici, la dynamique devient très complexe au moment de la disparition sur toute une plage de valeurs des paramètres. Dans l’espace produit de R² par l’espace des paramètres, cette complexité est de type hamiltonien (cf. chap. 4, en particulier fig. 39) et peut conduire à des courbes de Birkhoff (souvent appelées attracteurs de Birkhoff) analogues à celles que nous rencontrerons dans la deuxième partie de ce chapitre.

Nous considérons maintenant le cas où, dans l’équation (N⁵), les termes d’excitation 漣 f (x ) dx /dt et g (t ) sont très grands.

On supposera que f (x ) = 猪a 﨏(x /a ) pour x compris entre 漣 神 et + 神, où 﨏 est une fonction paire dérivable très proche de la fonction valant 漣 1 entre 漣 1 et + 1 et + 1 ailleurs, a est assez petit, et 猪 est très grand. De même, on choisira g (t ) = 猪塚p (t ), où p est une fonction impaire dérivable T-périodique, très proche de la fonction T-périodique qui vaut 漣 1 entre 漣 T/2 et 0 et + 1 entre 0 et T/2 (fig. 53).

Le changement de variables:

transforme (Q⁵) en l’équation:

où: 淋(X) = 咽X0 﨏(t ) dt est représentée sur la figure 54 (pour des valeurs de X comprises entre 漣 神/a et 神/a ).

Si a est assez petit et si 諸 = 1 (ce qu’on peut supposer par un changement de mesure du temps), l’équation (S⁵) est pratiquement identique à l’équation suivante, étudiée par Levinson dans les années quarante à la suite des travaux de Cartwright et Littlewood sur l’équation de Van der Pol, et plus récemment par Levi:

où 﨎 礪 0 est petit; cette équation équivaut encore à:

Lorsque, 﨎 étant fixé, b varie dans un certain intervalle [b ¹, b ²] (auquel nous nous limiterons), l’application de Poincaré P^b de (L) semble avoir beaucoup en commun avec l’application P size=1^猪, size=1^﨎 de l’équation (N⁵) pour 猪 礪 0: il existe en effet un point fixe répulsif X^r (spectre de la dérivée DP^b (X^r ) tout entier hors du disque unité de C) entouré d’un anneau invariant A (c’est-à-dire P^b (A) 說 A) qui attire en particulier tous les points «repoussés» par X^r (fig. 55, très schématique, où les flèches indiquent très grossièrement l’action de P^b ).

Mais, alors que dans le cas de l’équation (N⁵) l’intersection:

des itérés positifs de A est (pour 猪 礪 0 assez petit) une courbe fermée très différentiable, il n’en est plus de même pour l’équation (L): trois types de comportement (A, B, et g ) se partagent l’intervalle [b ¹, b ²] (fig. 56, où la somme des longueurs des g ⁱ tend vers 0 avec 﨎).

Lorsque b appartient à l’un des intervalles Aⁱ , il ne se passe rien d’autre que ce qu’on attendrait dans (N⁵) lorsque 猪 grandit: les accrochages de fréquences font perdre de plus en plus de régularité à la courbe invariante, jusqu’à ce que l’intersection des itérés de l’anneau ressemble à la figure 57, que l’on comparera à la figure 51, mais également à la figure 13. Sur la figure 57, l’intersection:

coïncide avec l’union de l’orbite périodique attractante et de la variété instable de l’orbite périodique de type col (union des variétés instables pour P^{b 2n i +1} des différents points de l’orbite).

L’unique attracteur de P^b est dans ce cas l’orbite périodique de type foyer; la période commune des deux orbites périodiques est grande ((2 n ⁱ + 1) avec n ⁱ de l’ordre de 1/ 﨎), constante sur Aⁱ , et décroît quand on passe de Aⁱ à A^{i +1} (remplacement de n ⁱ par n ^{i +1} = n ⁱ 漣 1). Nous ne parlerons pas des g ⁱ , zones complexes de bifurcation.

Lorsque b appartient à l’un des intervalles Bⁱ , une nouvelle paire (col-foyer) d’orbites périodiques apparaît, de période (2 n ⁱ 漣 1). L’existence dans l’anneau d’orbites périodiques de périodes minimales 2 n ⁱ 漣 1 et 2 n ⁱ + 1 est l’indice d’une grande complexité topologique de

par exemple, cette intersection n’est sûrement pas homéomorphe à un cercle, car le nombre de rotation de la restriction de P^b serait un rationnel de dénominateur égal à la fois à 2 n ⁱ 漣 1 et à 2 n ⁱ + 1!

Par analogie (mais seulement par analogie) avec une construction faite par Birkhoff en 1932, nous appellerons cette intersection une courbe de Birkhoff ; ce n’est pas un «attracteur» au sens actuel du mot (cf. chap. 6), même si ce sens est encore mouvant: pour presque tous les points z , les itérés Pⁿ (z ) tendent, lorsque n tend vers + 秊, vers l’une des deux orbites périodiques attractantes.

Les points de

dont les itérés ne tendent pas vers l’une de ces deux orbites forment un ensemble de Cantor invariant C (qui contient les deux orbites périodiques de types col associées), et la courbe de Birkhoff

n’est autre que la réunion des deux orbites périodiques attractantes et de la variété instable de C (définie chap. 7: C est un exemple d’ensemble invariant hyperbolique): c’est une courbe étrange qui, bien que topologiquement distincte d’un cercle, sépare le plan de la même façon que le ferait un cercle, ce qui ne l’empêche pas localement de ressembler au produit d’un ensemble de Cantor par un intervalle. Nous n’essaierons pas de tracer ce joli monstre (le lecteur déçu comprendra en regardant les dessins de Marie Charpentier dans le Journal de mathématiques pures et appliquées , 1935); la figure 58 (prise dans l’article de Levi) représente simplement ses sous-ensembles invariants (qu’on a reproduits également dans le cas où b est dans Aⁱ , ce qui permet la comparaison avec la fig. 57).

L’existence de C et la dynamique de P^b |C (en particulier l’existence d’une infinité d’orbites périodiques, et d’une infinité d’orbites chaotiques ), mais aussi la stabilité structurelle de P^b se prouvent par des méthodes de dynamique symbolique analogues à celles qu’on a utilisées au chapitre 4 pour élucider la structure associée aux points homoclines. Puisque c’est en réfléchissant à cet exemple que Smale inventa le fer à cheval, cela ne doit pas nous étonner. En fait, l’application de ces méthodes demande un très gros travail, et (L) est une des rares équations explicites en dimension trois dont on soit capable d’exhiber précisément le comportement géométrique.

L’idée exploitée par Levi est la suivante: il commence par montrer la présence dans A d’un petit rectangle r (fig. 59) ayant la propriété que, pour tout point z de A, il existe un entier n (z ) 礪 0 tel que P^{b n (z )}(z ) soit dans r ; il montre alors qu’avant de revenir dans r les images de r par les itérés de P^b sont étirées et pliées comme sur la figure 59, ce qui commence à faire penser au fer à cheval...

Le phénomène vient de ce qu’au voisinage du point V (fig. 59) le champ de vecteurs dépendant du temps associé à (L) a un comportement complexe dû à la présence du segment UV séparant une zone où le champ est dirigé vers la droite (et très grand) d’une zone où il est dirigé vers la gauche (et également très grand). Cela est analogue à ce qui se passe dans l’effet canard discuté au chapitre 3 et peut être analysé à l’aide de l’analyse non standard ^{[cf. ANALYSE NON STANDARD]}.

Pour finir, notons que, au bout de suffisamment d’itérations, l’image de A ressemble beaucoup à une courbe fermée; on peut alors penser à P^b comme à une application du cercle sur le cercle. La différence avec (N⁵) est que cette application n’est plus bijective; en particulier, l’unique nombre de rotation est remplacé par un intervalle de nombres de rotation, ce qui autorise des orbites périodiques de périodes différentes.

La réduction (au moins heuristique) d’une équation différentielle en dimension trois à une application non injective (endomorphisme) en dimension un, éventuellement discontinue, est assez fréquente; nous la rencontrerons encore au chapitre 6 à propos de l’attracteur de Lorenz (à la différence près que dans ce cas on ne sait pas vraiment montrer que le modèle géométrique reflète bien les propriétés de l’équation).

Excursion homologique: nombre de rotation et cycles asymptotiques

Ainsi que l’a remarqué Herman, le nombre de rotation 福(f ) d’un homéomorphisme f du cercle S¹ conservant l’orientation peut être défini par 福( face="EU Caron" キ ) = 猪( 﨏), où face="EU Caron" キ = Id + 﨏 est un relèvement de f en un difféomorphisme de R, et 猪 une mesure de probabilité sur le cercle (mesure de Radon de masse totale égale à 1) invariante par f ; une telle mesure existe toujours, mais n’est unique que si 12 神 練 福( face="EU Caron" キ ) est irrationnel.

Si X est un espace métrique connexe compact ayant, par exemple, le type d’homotopie d’une variété, et f un homéomorphisme de X isotope à l’identité qui préserve la mesure de probabilité 猪, on peut leur associer de façon analogue une classe de rotation 福 size=1^猪(f ) appartenant au quotient du premier groupe d’homologie H¹(X, R) de X par un certain sous-groupe 臨 contenu dans H¹(X, Z) (et donc discret). Par exemple, lorsque X est le tore Tⁿ de dimension n , 臨 coïncide avec H¹(Tⁿ , Z) 黎= Zⁿ 說 Rⁿ 黎= H¹(Tⁿ , R), et la classe de rotation 福 size=1^猪(f ) peut être naturellement considérée comme un élément de Tⁿ = Rⁿ /Zⁿ ; si 猪 définit le volume sur Tⁿ la nullité de 福 size=1^猪(f ) s’interprète comme la conservation par f du centre de gravité au sens de l’excursion variationnelle du chapitre 4.

Plus généralement, si X est une variété différentiable orientable de dimension n et si 猪 définit le volume sur X (mesure de Lebesgue normalisée), la dualité de Poincaré associe à 福 size=1^猪(f ) l’élément de H^{n -1}(X, R) qui fait correspondre à chaque (n 漣 1)-cycle le volume total balayé par ses images par une isotopie de f à l’identité.

Ces classes de rotation, définies sous leur forme la plus générale par Fathi et Visetti, avaient été primitivement construites par Schwartzman dans le cas des flots sous le nom de cycles asymptotiques et utilisées par ce dernier pour caractériser les flots possédant une section globale (hypersurface rencontrant transversalement chaque courbe intégrale). On peut dans ce cas les décrire simplement en termes des courants géométriques de Ruelle et Sullivan.

6. Couplage d’auto-oscillations et attracteurs étranges

La simplicité des attracteurs que nous avons rencontrés jusqu’ici (sous-espaces homéomorphes à des points, des cercles ou des tores) nous a dispensés de donner de ce terme une définition plus précise qu’«ensemble des mouvements asymptotiques dans un système dissipatif». Même dans le pendule forcé que nous avons étudié à la fin du chapitre 5, la plupart des courbes intégrales de la région considérée tendaient asymptotiquement vers l’une ou l’autre de deux orbites périodiques, et la complexité de l’ensemble invariant que nous avons trouvé pour l’application de Poincaré (localement le produit d’une droite par un ensemble de Cantor) reflétait seulement l’intrication des bassins d’attraction de ces orbites. Dans les exemples de ce chapitre, c’est l’attracteur lui-même qui sera responsable de la complexité et un essai de définition vaut d’être tenté: on appellera attracteur d’un système dynamique (à temps continu ou discret) un sous-ensemble fermé invariant (c’est-à-dire une réunion d’orbites) de l’espace de phase, attirant tous (presque tous serait plus raisonnable) les points d’un voisinage (le bassin de l’attracteur ) lorsque le temps croît, et ayant une certaine forme d’indécomposabilité, par exemple la possession d’une orbite dense (c’est-à-dire une orbite venant arbitrairement près de chaque point du sous-ensemble). Si, comme pour les courbes de Birkhoff, cette dernière condition n’est pas remplie, on parlera d’«ensemble invariant attractant». Notons que cette quasi-homogénéité implique qu’un attracteur est, au voisinage de chacun de ses points, ou bien très simple (par exemple une variété), ou bien très compliqué (au mieux un ensemble de Cantor). Suivant une tradition dont l’origine est mal établie, l’adjectif étrange sera accolé à tout attracteur ne ressemblant pas à ceux qu’on a déjà rencontrés dans les chapitres précédents. Cela implique souvent une dépendance sensible des conditions initiales (croissance exponentielle avec le temps des erreurs sur les données), propriété quelquefois exigée (avec raison) comme certificat d’étrangeté et souvent assimilée (à tort) à la notion de comportement chaotique (comme le dit Thom, quoi de moins chaotique qu’un difféomorphisme d’Anosov linéaire sur le tore, cf. infra . 7).

Peu de résultats théoriques concernent la structure et la classification des attracteurs, mais chaque jour apporte de nouvelles équations susceptibles de posséder un attracteur étrange, issues de la physique (turbulence), de la chimie (réactions oscillantes), de l’économie... Une énorme masse d’expériences numériques commence à s’accumuler, expériences dont l’interprétation même est souvent sujette à caution: comment distinguer par exemple un honnête attracteur étrange d’une perfide orbite périodique de très longue période?

Avant de passer en revue quelques exemples, nous allons montrer que tous peuvent se rencontrer dans le système formé par n pendules non linéaires faiblement couplés, c’est-à-dire dans les petites perturbations du système d’équations:

où f size=1^猪 est défini comme dans l’équation (N³) du chapitre 3.

L’espace de phase est maintenant le produit de n exemplaires de S¹ 憐 R, c’est-à-dire Tⁿ 憐 Rⁿ (Tⁿ = (S¹)ⁿ = tore de dimension n ) qu’on peut dérouler (dans n directions) sur R²ⁿ . Si les 猪ⁱ sont positifs, l’étude du chapitre 3 montre que le champ de vecteurs défini par (N³ⁿ ) possède un ensemble invariant attractant homéomorphe au tore Tⁿ ; d’autre part, pour des 猪ⁱ bien choisis, ce tore est un attracteur (toutes ses orbites sont denses et les mouvements asymptotiques sont quasi périodiques). Si nous introduisons un couplage suffisamment faible entre les pendules, par exemple en les reliant par des ressorts très lâches, un tore attractant homéomorphe à Tⁿ continue à exister mais il est possible que sur ce tore le flot possède maintenant un attracteur étrange qui attire asymptotiquement presque toutes les courbes intégrales dans l’espace de phase. Le comportement asymptotique de ces dernières semblera erratique et il sera bien difficile de se persuader que le système observé obéit encore aux règles d’un strict déterminisme, ce qu’il fait cependant (il ne faut pas confondre déterminisme et possibilité pratique de prévision). Un tel exemple a été proposé pour la première fois en 1971 par Ruelle et Takens comme un modèle possible d’apparition de la turbulence hydrodynamique, à un moment où la pertinence même des équations de Navier-Stokes (cf. équations aux DÉRIVÉES PARTIELLES) dans les phénomènes turbulents était très discutée: le fait que des équations en dimension finie (et pas trop grande) puissent sous-tendre un comportement asymptotique aussi complexe peut faire espérer qu’une équation d’évolution aux dérivées partielles (qu’il faut penser heuristiquement comme une équation différentielle sur un espace de phase de dimension infinie) soit suffisamment riche pour contenir les phénomènes de turbulence.

Nous nous contenterons, suivant Ruelle, Takens et Newhouse, de montrer que par une perturbation arbitrairement petite (ainsi que ses dérivées) d’un champ de vecteurs quasi périodique sur Tⁿ (n 閭 4) (par exemple celui qui est défini par (N³ⁿ )), on peut introduire pratiquement n’importe quel attracteur étrange. Le champ de vecteurs sur ce tore est, pour un bon choix de coordonnées ¹, ..., ⁿ , de la forme:

(on considère les ⁱ comme des réels, en déroulant Tⁿ sur Rⁿ : un champ de vecteurs sur Tⁿ correspond simplement à un champ de vecteurs sur Rⁿ 2 神-périodique par rapport à chaque variable). Au prix d’une perturbation arbitrairement petite, on peut supposer les réels 諸ⁱ de la forme 諸ⁱ = ck ⁱ , c 捻 R, k ⁱ 捻 Z, i = 1, 2, ..., n , les k ⁱ étant de surcroît linéairement indépendants. Un difféomorphisme (changement de coordonnées) du tore, défini par un isomorphisme linéaire unimodulaire de Rⁿ , transforme alors ce champ de vecteurs en un champ de la forme:

dont l’application de Poincaré P: T^{n -1} T^{n -1}, définie en regardant le premier retour des courbes intégrales sur le tore 﨡ⁿ = 0, est l’identité. Il suffit donc de savoir approcher l’identité du tore T^{n -1} par un difféomorphisme 磊 ayant un attracteur étrange. Cela peut être fait à partir d’un difféomorphisme f du tore T^{n -2} isotope à l’identité , c’est-à-dire joint à l’identité par une famille à un paramètre (f ^t )^t size=1^{捻 [0,1]} de difféomorphismes, f ⁰ = Id, f ¹ = f , et ayant un attracteur étrange: soit ( 﨏^t )^t size=1^{捻 R} l’unique flot sur T^{n -2} 憐 S¹ = T^{n -1} tel que 﨏^t ( 﨡¹, ..., 﨡^{n -2}, 0) = (f ^t ( 﨡¹, ..., 﨡^{n -2}), 2 神t ) (pour éviter des problèmes de différentiabilité, on suppose f ^t = f ⁰ pour t voisin de 0 et f ^t = f ¹ pour t voisin de 1); son application de Poincaré sur T^{n -2} 憐0 n’est autre que la restriction de 﨏² size=1^神, c’est-à-dire f . Si 﨎 est petit, le difféomorphisme 﨏 size=1^﨎 est proche de l’identité et commute (via la projection de T^{n -2} 憐 S¹ sur S¹) avec la rotation 﨡^{n -1} 料 﨡^{n -1} + 2 神﨎 de S¹.

On obtient 磊 en choisissant 﨎 = 1/N, N entier très grand, et en remplaçant cette rotation par un difféomorphisme h proche n’ayant que deux orbites périodiques de période N, l’une attractante, l’autre répulsive (face="EU Upmacr" 磊 devient un produit gauche au-dessus de h ). Le difféomorphisme ainsi construit est très proche de l’identité si N est grand, et possède un attracteur homéomorphe à la réunion disjointe de N copies de l’attracteur de f .

Dans la suite du chapitre, nous donnons quelques exemples classiques de difféomorphismes et de flots ayant (ou semblant avoir) un attracteur étrange. Ces exemples sont décrits dans un espace euclidien mais peuvent être plongés dans n’importe quelle variété de même dimension. Au lecteur de reconnaître si les difféomorphismes correspondant sont ou non isotopes à l’identité.

Exemple 1: le solénoïde de Smale . Il s’agit d’un difféomorphisme de R³ qui applique un tore solide dans lui-même après l’avoir étiré et replié (fig. 60, où l’on a représenté une coupe z = C ^te ). Plus précisément, si on identifie le cercle S¹ à l’ensemble des nombres complexes z de module égal à 1 et le disque D² à l’ensemble des nombres complexes w de module inférieur ou égal à 1, la restriction de f au tore solide S¹ 憐 D² peut être définie par

C’est un produit gauche au-dessus de l’application du cercle sur lui-même définie par zz ² (doublement des angles).

Par récurrence, on voit que f ⁿ (S¹ 憐 D²) s’enroule 2ⁿ fois autour du tore, coupant chaque disque d’équation z = constante transversalement en 2ⁿ disques de rayon (1/4)ⁿ . L’attracteur

intersecte chaque disque suivant un ensemble de Cantor et est localement homéomorphe au produit de l’ensemble de Cantor par un intervalle (il est en fait réunion d’une famille non dénombrable de courbes). La figure 61, tirée d’un article de Ruelle, indique l’orbite d’un point X⁰ arbitraire du tore solide: les itérés f ⁿ (X⁰) se rapprochent de l’attracteur et, pour n assez grand, appartiennent pratiquement à l’attracteur (bien entendu ils ne peuvent y appartenir en toute rigueur que si l’on part d’un point X⁰ 捻 A). La figure 61 représente f ⁿ (X⁰) pour 51 諒 n 諒 5 050.

Le solénoïde A est un bon exemple d’ensemble invariant hyperbolique d’un difféomorphisme: en chaque point il existe une direction (celle du disque) dans laquelle f contracte et une direction (celle du cercle, le long de A) où f dilate, et ces directions dépendent continûment du point (cf. chap. 7); c’est aussi un exemple idéalement simple d’attracteur étrange: il attire tout un voisinage (son bassin), il possède une orbite dense, et il présente le phénomène de «dépendance sensible des conditions initiales»: même si X⁰ et Y⁰ sont initialement très proches, les destins de f ⁿ (X⁰) et f ⁿ (Y⁰) n’ont en général rien en commun, ce qui rend toute prédiction à long terme impossible. Ces deux derniers points peuvent se voir en utilisant une méthode de dynamique symbolique tout à fait analogue à celle qui a permis, au chapitre 4, d’analyser le fer à cheval (qui, lui, n’est pas un attracteur): on découpe le tore solide en deux parties ⁺⁰ =(e ⁱ size=1, w ), 0 諒 諒 神 et ⁺¹ =(e ⁱ size=1, w ), 神 諒 諒 2 神 et on associe à un point X⁰ la suite (iⁿ )ⁿ size=1^{捻 Z} définie par i ⁿ = 0 si f ⁿ (X⁰) 捻 ⁺⁰, i ⁿ = 1 si f ⁿ (X⁰) 捻 ⁺¹; il y a bien sûr ambiguïté si, par hasard, f ⁿ (X⁰) 捻 ⁺⁰ 惡 ⁺¹, ce qui ne pouvait arriver dans le cas du fer à cheval.

Étant donné une suite finie 見⁰, ..., 見^{n -1} de zéros et de uns, l’ensemble ⁺ size=1^見0... size=1^{見n -1} des X⁰ 捻 S¹ 憐 D² tels que f ⁱ (X⁰) 捻 intérieur de ⁺ size=1^見i , i = 0, ..., n 漣 1, est un secteur d’angle 2 神/2ⁿ ; l’ensemble ^- size=1^見-n... size=1^見-1 des X⁰ 捻 S¹ 憐 D² tels que f ⁱ (X⁰) 捻 intérieur ⁺ size=1^見i , i = 漣 1, 漣 2, ..., 漣 n , est un tube de rayon (1/4)ⁿ faisant une fois le tour de S¹ 憐 D². L’intersection size=1^見-n... size=1^見-1; size=1^見0... size=1^{見n -1} = ^- size=1^見-n... size=1^見-1 惡 ⁺ size=1^見0... size=1^{見n -1} définit donc, lorsque n tend vers l’infini, un unique point (z , w ) qui appartient à A (en particulier, l’argument de z s’écrit (0, 見⁰, 見¹...) 練 2 神 en base 2).

On vient de définir une application 崙 de0, 1^Z sur A qui, comme au chapitre 4, vérifie f 獵 崙 = 崙 獵 靖, où 靖 est l’homéomorphisme de décalage de virgule . Cette application est surjective mais non injective (elle l’était pour le fer à cheval); cependant, l’image réciproque d’un point de A ne contient jamais plus de deux éléments de0, 1^Z et 崙 rend pratiquement les mêmes services que si elle était bijective pour montrer l’existence d’orbites périodiques, d’orbites denses, etc.

Nous quitterons le solénoïde en signalant à son propos la théorie de Williams qui permet de comprendre la topologie et la dynamique des attracteurs hyperboliques ayant, comme A, la propriété que la restriction de f à l’attracteur est une dilatation: soit p ¹ la projection de S¹ 憐 D² sur S¹ et 神 l’application de A dans le produit

d’une infinité de cercles définie par 神(X) = (p ¹ 獵 f ^-i (X))⁺ size=1^{秊i = 0}; on peut montrer que 神(A) coïncide avec l’ensemble:

et que 神 est un homéomorphisme de A sur L qui conjugue f |A à l’homéomorphisme g de L défini par g ((z ⁱ )⁺ size=1^{秊i = 0}) = (z ^{i 2})⁺ size=1^{秊i = 0}.

Pour les mathématiciens, L n’est autre que la limite projective L = 良lim (S¹ 﨏 良 S¹ 﨏 良 S¹ 良...), avec 﨏(z ) = z ².

Dans le cas général, 﨏 : S¹S¹ est remplacé par une application continue dilatante 﨏: KK d’une variété branchue K, et la même construction de limite projective permet de retrouver l’attracteur et la dynamique sur l’attracteur (fig. 69 et 70).

Exemple 2: l’attracteur de Lorenz . En 1963, le météorologue Edward Lorenz publia un article qui, à cette époque, échappa complètement aux mathématiciens: il y écrivait la famille d’équations différentielles dans R³ ci-dessous, obtenue par une troncation «sauvage» à partir des équations de la convection ( 靖 correspond à un nombre de Prandtl, r à un nombre de Rayleigh, et b à un facteur de forme):

et montrait numériquement qu’elle semblait, pour 靖 = 10, r = 28, b = 8/3, posséder un attracteur étrange (le mot n’était pas encore employé). La figure 62 (obtenue par Lanford) montre un morceau de la courbe intégrale passant par l’origine: cette courbe fait un certain nombre de tours 見 sur la boucle de droite, puis un certain nombre 廓 à gauche, et ainsi de suite; si on change un tant soit peu la donnée initiale, on obtient une courbe intégrale analogue, mais la suite de nombres 見, 廓, ... est complètement différente, manifestant la dépendance sensible des conditions initiales (ce que Lorenz appelait l’effet aile de papillon , arguant que le battement des ailes d’un papillon pouvait suffire à modifier complètement les conditions météorologiques quelque temps après).

On ne sait toujours pas montrer que l’équation (M) possède un attracteur étrange, mais Guckenheimer et Williams ont construit des modèles géométriques analogues sur lesquels tout peut être démontré. La structure du flot est très riche; il y a par exemple une infinité d’orbites périodiques nouées, c’est-à-dire plongées dans R³ de façon non triviale (Birman et Williams); de plus, le type topologique de l’attracteur peut varier continûment, ce qui montre sa complexité. Ces attracteurs de Lorenz géométriques sont obtenus par une construction de limite inverse analogue à celle qui fournit le solénoïde. Leur étude se ramène en deux étapes à celle d’un endomorphisme (application non inversible) discontinu d’un intervalle: existence sur un plan bien choisi d’une application de Poincaré P (non partout définie), puis quotient de P par un feuilletage invariant. Remarquons enfin que de nombreux exemples expérimentaux d’attracteurs étranges sont fondés sur une telle suite infinie de choix (ici, de quel côté tourner?).

Exemple 3: l’attracteur de Hénon et quelques autres . Si on cherche à transformer l’exemple du fer à cheval en un attracteur, il est naturel de considérer des applications du plan dans le plan appliquant un carré dans lui-même après étirement et repliement. Un exemple de ce type (au renversement de l’orientation près) a été proposé par Hénon en 1976 à la suite de calculs numériques de Pomeau et Ibanez sur les équations de Lorenz (pour des valeurs des paramètres autres que celles que nous avons citées). Il s’agit de l’application, maintenant célèbre,

La figure 63 (tirée de l’article de Hénon) montre 10 000 points de l’orbite de (0,631 354 48; 0,189 406 34) lorsque a = 1,4, b = 0,3, ainsi que trois agrandissements successifs de portions de cette orbite.

L’orbite semble parcourir un attracteur étrange localement homéomorphe au produit d’un ensemble de Cantor par un intervalle, mais cela n’a pas pu être démontré pour le moment: Misiurewicz et Szenc ont seulement montré l’existence d’un point homocline, et donc d’un ensemble de Cantor invariant; Misiurewicz a réussi par contre à prouver l’existence d’un attracteur étrange dans un modèle linéaire par morceaux de ce difféomorphisme, proposé par Lozi. Des doutes ont été émis sur l’existence effective d’un attracteur étrange, en particulier par Newhouse qui avance la possibilité d’une orbite périodique attractante de très longue période; remarquons qu’en changeant un peu a et b on peut obtenir comme attracteur des orbites périodiques de petites périodes.

D’autres exemples, généralisant celui de Henon, ont été proposés dans le même esprit, en particulier les difféomorphismes de la forme (A.C.T.), (x , y ) 料 (y , f size=1^猪(y ) + bx ), étudiés par Arneodo, Coullet et Tresser (si f size=1^猪(y ) = 猪y (l 漣 y ), on retrouve l’exemple de Henon). Dans la limite de contraction infinie (b = 0), l’image de R² devient le graphe de l’endomorphisme f size=1^猪. Le comportement dynamique étonnamment complexe de cette famille à un paramètre d’endomorphismes (que l’on peut restreindre à l’intervalle [0, 1]) a été étudié récemment par une foule d’auteurs (Sharkovski, Feigenbaum, Coullet, Tresser, Misiurewicz, Collet, Eckman, Lanford, Milnor, Thurston, Bowen, Guckenheimer...).

Le problème de l’itération des endomorphismes polynomiaux (et des fonctions rationnelles) de C avait été étudié par Julia et Fatou au début du siècle et vient d’être repris par Douady, Hubbard, Sullivan, Mañe, Sad (fig. 65). Nous ne ferons cependant qu’effleurer le sujet dans le cas réel. Tant que 猪 est inférieur ou égal à 1, l’origine est un attracteur global: pour tout x dans [0, 1], lim f ⁿ (x ) = 0. Quand 猪 augmente, n 轢+ 秊l’origine devient un point fixe répulsif et un nouveau point fixe attractant apparaît pour 猪 = 猪¹; puis ce nouveau point fixe devient à son tour répulsif et est remplacé, lorsque 猪 dépasse la valeur 猪², par une orbite périodique attractante de période 2 qui elle-même devient répulsive lorsque 猪 dépasse 猪³ et est remplacée par une orbite périodique attractante de période 2², etc. jusqu’à ce que, pour une valeur 猪 size=1^秊 = lim 猪ⁿ = 3,569 94, f size=1^猪 size=1^秊 ait n 轢秊

pour attracteur un ensemble de Cantor (adhérence d’une «orbite périodique de période 2 秊»).

Cet exemple montre bien la difficulté qu’il y a à définir le terme attracteur: l’ensemble de Cantor n’attire certainement pas tous les points d’un voisinage, puisqu’il est limite d’orbites périodiques répulsives. Des essais de définition, adaptés au comportement d’un système sous l’effet de petites perturbations aléatoires, ont été récemment proposés par Conley et Ruelle.

La figure 64 illustre cette suite de bifurcations dans l’espace produit R 憐 [0, 1] (les pointillés indiquent les répulseurs, les pleins indiquent les attracteurs; le dessin n’est pas à l’échelle).

La suite 猪ⁿ converge vers 猪 size=1^秊 comme une série géométrique dont la raison est universelle parmi les familles du même type (ayant en particulier un unique maximum non dégénéré), phénomène dont l’explication est fournie par des arguments de renormalisation analogues à ceux qui, en physique, prouvent l’universalité des exposants critiques. De même, la dynamique de f size=1^猪 size=1^秊 sur son attracteur est universelle: c’est l’adding machine des ergodiciens (addition de 1 sur les entiers dyadiques); on peut en déduire que f size=1^猪 size=1^秊 est topologiquement conjugué à la restriction de f ² size=1^猪 size=1^秊 à un petit intervalle. Lorsque 猪 dépasse la valeur 猪 size=1^秊, des orbites périodiques de période impaire apparaissent, et la dynamique devient beaucoup plus compliquée. L’ordre dans lequel les différentes périodes apparaissent en fonction de 猪 est lui aussi universel (ordre de Sharkovski); par exemple, s’il existe une orbite périodique de période 3, il existe des orbites périodiques de toutes les périodes (comparer avec les homéomorphismes du cercle où une seule période existe à la fois!): ce dernier théorème a été retrouvé par Li et Yorke, et baptisé du nom suggestif de Period three implies chaos . Les méthodes utilisées sont de type dynamique symbolique (invariants de pétrissage de Milnor et Thurston). Récemment, Jakobson a montré l’existence d’un ensemble de mesure positive de valeurs de 猪 pour lesquelles f size=1^猪 possède une mesure invariante absolument continue par rapport à la mesure de Lebesgue, ce qui implique un comportement asymptotique très complexe.

Des exemples de flots fondés sur ces endomorphismes ont également été étudiés, et tout une zoologie d’attracteurs étranges numériques provient d’équations décrivant des oscillations forcées avec feedback de la forme:

c’est-à-dire

Nous rencontrerons au chapitre 7 d’autres exemples d’attracteurs, en particulier les difféomorphismes d’Anosov sur le tore pour lesquels le tore entier est un attracteur. Il ne faudrait pas croire que les seules variétés pouvant être un attracteur soient des tores: une condition nécessaire et suffisante est l’existence d’un flot ou d’un difféomorphisme sur la variété ayant une orbite dense; nous avons déjà signalé les exemples d’Anosov et Katok de difféomorphismes du disque; mentionnons également le théorème de Fathi et Herman montrant l’existence sur la sphère S³ de dimension trois d’un difféomorphisme minimal (c’est-à-dire dont toutes les orbites sont denses), en contraste avec la conjecture de Seifert selon laquelle tout champ de vecteurs C size=1^秊 sur S³ aurait au moins une orbite périodique (un contre-exemple en classe C¹ a été donné par Schweitzer à partir des contre-exemples de Denjoy déjà cités).

Finissons ce chapitre avec un peu de théorie en mentionnant deux outils topologiques permettant de mesurer la complexité d’un système dynamique, la dimension de Hausdorff d’un sous-ensemble invariant et l’entropie topologique .

Étant donné un compact K de Rⁿ et deux réels positifs d , 﨎, on note:

où I parcourt l’ensemble des familles finies de boules de rayons r ⁱ 諒 﨎, i 捻 I, dont la réunion recouvre K. La limite 猪(K, d ) de la fonction décroissante 﨎 料 猪(K, d , 﨎) lorsque 﨎 tend vers 0 est la mesure de Hausdorff d-dimensionnelle de K. Il existe un unique réel dim (K) 閭 0, appelé dimension de Hausdorff de K, tel que 猪(K, d ) = + 秊 si d 麗 dim (K) et 猪(K, d ) = 0 si d 礪 dim (K). Lorsque K est une sous-variété de Rⁿ , d coïncide avec sa dimension usuelle et est donc entier, mais la réciproque n’est pas vraie: par exemple, un sous-ensemble de R de mesure de Lebesgue positive a une dimension de Hausdorff égale à 1. L’ensemble de Cantor triadique décrit au chapitre 4 (fig. 36) a une dimension de Hausdorff égale à Log 2/Log 3, et le solénoïde de l’exemple 1 a une dimension de Hausdorff égale à 3/2. Sous le nom d’objets fractals , les compacts de R² de dimension de Hausdorff strictement comprise entre un et deux ont fait récemment l’objet d’une grande campagne de promotion initiée par Mandelbrot. En dehors des attracteurs et des courbes de Sierpinski, on en rencontre de très beaux dans les frontières de bassins d’attraction, en particulier dans l’étude dynamique des polynômes d’une variable complexe (fig. 65, due à Maryvonne Daguenet-Teissier).

Un théorème démontré par Douady et Oesterlé en 1980, et généralisé par Ledrappier, permet de borner la dimension de Hausdorff des ensembles invariants attractants d’un difféomorphisme f en fonction des exposants caractéristiques de f définis au chapitre 8. Ce théorème est valable en dimension infinie et concerne donc certaines équations aux dérivées partielles.

Quant à l’entropie topologique d’une application continue f d’un espace métrique M, elle mesure en quelque sorte le taux de croissance asymptotique du nombre d’orbites de f qu’on distingue par une observation de plus en plus fine; plus précisément, appelons s (n , 﨎) le nombre maximal d’éléments d’un sous-ensemble E de M ayant la propriété que deux points x et y quelconques de E vérifient: distance (f ^k (x ), f ^k (y )) 礪 﨎 pour au moins un entier positif k inférieur à n , et définissons:

L’entropie topologique de f est la limite h (f ) = lim h (f , 﨎). C’est un invariant de 﨎 轢0conjugaison topologique de f (si g est un homéomorphisme de M, h (g 獵 f 獵 g ^-1) = h (f )), dont on peut montrer que c’est la borne supérieure des entropies métriques h size=1^猪(f ) (cf. théorie ERGODIQUE) pour toutes les mesures de probabilité 猪 sur M invariantes par f . Si la dimension de M est un ou deux, il y a une forte relation entre la condition de positivité de h (f ) et l’abondance des points périodiques: théorème de Bowen et Franks montrant la positivité de h (f ) pour tout endomorphisme f de l’intervalle ayant un point périodique de période impaire, théorème de Katok affirmant qu’un difféomorphisme f de classe C² d’une surface vérifiant h (f ) 礪 0 possède nécessairement des points périodiques hyperboliques ayant des points homoclines transversaux, théorème de Franks et Williams affirmant qu’un flot ( 﨏^t )^t size=1^{捻 R} sur la sphère de dimension trois vérifiant h ( 﨏¹ ) 礪 0 possède nécessairement une infinité d’orbites périodiques nouées d’une infinité de façons (comparer à l’attracteur de Lorenz). Par contre, en dimension supérieure, un difféomorphisme peut avoir une entropie topologique positive sans posséder la moindre orbite périodique (Herman a même construit de tels f qui sont minimaux, c’est-à-dire dont chaque orbite est dense).

Enfin de nombreux travaux tournent autour de l’entropy conjecture de Shub qui concerne les liens existant entre l’entropie topologique de f et sa topologie algébrique. Un exemple en est le théorème de Manning qui, étant donné une application continue f d’une variété compacte M dans elle-même, minore h (f ) par le logarithme de la plus grande valeur propre de l’application induite par f sur le premier groupe d’homologie de M.

Un autre est le théorème de Ruelle et Sullivan, application de la théorie des courants géométriques , montrant que, sous une condition minime d’orientabilité, un difféomorphisme d’Anosov f (cf. chap. 7) a la propriété que e ^{h (b )} est une valeur propre de l’application qu’il induit sur le groupe d’homologie correspondant à sa dimension dilatante («instable»).

7. Hyperbolicité

Les exemples des chapitres précédents ont fait apparaître la possibilité dès la dimension trois (dimension deux pour les difféomorphismes, un pour les endomorphismes) d’un comportement asymptotique complexe de certains ensembles de courbes intégrales (orbites); de plus, ce comportement, loin d’être exceptionnel, pouvait comme dans l’exemple 1 posséder une certaine stabilité vis-à-vis des perturbations des équations.

Nous donnons dans ce chapitre un aperçu de la structure de ces mouvements asymptotiques pour les systèmes dynamiques vérifiant l’axiome A de Smale, seule grande classe dont existe une théorie globale satisfaisante: la stabilité y est due en dernier ressort à la stabilité des points fixes hyperboliques des difféomorphismes d’un espace vectoriel normé complet. Ce dernier fait est très général: si l’attracteur d’un difféomorphisme ne change pas sous de petites perturbations éventuellement différentes à chaque itération, il est hyperbolique au sens défini ci-dessous (théorème de Franks).

Cette théorie a son origine dans les travaux d’Anosov sur le flot géodésique des variétés riemanniennes à courbure négative, flot relié en particulier au mouvement d’une masse ponctuelle sur la variété; nous décrirons succinctement le cas des surfaces de courbure négative constante à la fin du chapitre, mais nous parlerons surtout des systèmes dynamiques à temps discret (difféomorphismes), le cas des flots étant parallèle. Les références accessibles sont l’article de Smale dans le Bulletin de l’A.M.S. , le petit livre de Bowen, et le cours de Lanford aux Houches; un exposé plus technique se trouve dans le livre de Shub.

Commençons par préciser la notion de comportement asymptotique en définissant l’ensemble non errant 行(f ) d’un difféomorphisme f d’une variété M comme le complémentaire du sous-ensemble de M formé des points possédant un voisinage U tel que f ^k (U) 惡 U soit vide pour tout entier k (points errants; dans le cas d’un flot, la condition devient 﨏^t (U) 惡 U vide pour t 礪 t ⁰). Il est clair que 行(f ) est fermé, invariant par f (c’est-à-dire f ( 行) = 行), et contient le sous-ensemble Per(f ) des points périodiques de f , donc aussi son adhérence Per(f ); il contient également l’ensemble limite L(f ), réunion des ensembles 諸-limite L size=1^諸(x ) et 見-limite L size=1^見(x ) de tous les points x de M (pour leur définition, cf. supra , chap. 3): en particulier, si M est compacte et si U est un voisinage de 行(f ), il existe pour tout x un entier N 礪 0 tel que f ⁿ (x ) appartienne à U dès que |n | est supérieur à N. Enfin, presque par définition, tout point de 行(f ) est 嗀-pseudopériodique pour tout 嗀 礪 0 dans le sens suivant: appelons 嗀-pseudo-orbite de f une suite (finie ou infinie) de points ..., x ⁿ , x ^{n +1}, x ^{n +2}, ... de M tels que d (f (x ⁿ ), x ^{n +1}) 諒 嗀 pour tout n (c’est une (portion d’) orbite si on se permet des erreurs de l’ordre de 嗀); un point x est dit 嗀-pseudopériodique s’il existe une 嗀-pseudo-orbite de la forme x ⁰ = x , x ¹, ..., x ^p = x .

Identifions les ensembles non errants des divers flots ou difféomorphismes rencontrés jusqu’ici: pour le pendule du chapitre 1, tout l’espace de phase est formé de points non errants, mais l’introduction du frottement au chapitre 2 réduit brutalement l’ensemble non errant à un couple de points (les deux équilibres); au chapitre 3, la non-linéarité permet à une orbite périodique de s’ajouter à ces deux points. Dans les deux chapitres suivants, l’ensemble non errant du flot se reconstitue par suspension à partir de celui du difféomorphisme de Poincaré P: dans la situation conservative du chapitre 4, le théorème de récurrence de Poincaré implique que 行(P) contienne tout un voisinage du point fixe elliptique; dans la situation dissipative de la fin du chapitre 5, 行(P) est évidemment contenu dans la courbe de Birkhoff. Bien entendu, tous les attracteurs du chapitre 6 font partie de l’ensemble non errant du difféomorphisme ou du flot correspondant.

Dans tout ce qui suit, nous supposerons que 行(f ) est compact. Introduisons maintenant la notion d’ensemble hyperbolique , déjà rencontrée sous les traits de point fixe hyperbolique, orbite périodique (normalement) hyperbolique, ensemble de Cantor invariant hyperbolique (dans le fer à cheval), attracteur hyperbolique (le solénoïde): un sous-ensemble 炙 de M invariant par f est dit hyperbolique si, pour tout point x de 炙, le comportement, par rapport à l’orbite de x , des orbites des points y de 炙 proches de x est analogue au comportement des orbites au voisinage d’un point fixe hyperbolique: l’espace tangent en x à M se décompose en somme directe d’un sous-espace stable E^{s x} correspondant aux directions dans lesquelles se trouvent des y proches dont l’orbite converge exponentiellement vers l’orbite de x , et d’un sous-espace instable E^{u s} (unstable en anglais) correspondant aux directions dans lesquelles se trouvent des y proches dont l’orbite diverge exponentiellement de l’orbite de x , au moins tant qu’elle en est assez proche; on demande de plus que cette décomposition dépende continûment de x , et que le taux de convergence ou de divergence vérifie des estimations uniformes.

Techniquement, on demande que la restriction à 炙 du fibré tangent à M admette une décomposition en somme directe topologique de deux sous-fibrés E^s et E^u invariants par l’application Tf tangente à f , et qu’il existe deux constantes C 礪 0, et 0 麗 麗 1, telles que, pour tout n 閭 0,

Cela suppose choisie sur M une métrique riemannienne; si 炙 est compact, cette condition est indépendante du choix de la métrique; il est d’ailleurs possible (et nous le ferons désormais) de choisir cette dernière de façon que C = 1.

La définition est la même pour un flot, à l’exception de l’existence d’une direction neutre (celle du champ de vecteurs).

Particulièrement important du point de vue historique (flots géodésiques) est le cas où la variété M elle-même est hyperbolique (on dit que f est un difféomorphisme d’Anosov ou que 﨏^t est un flot d’Anosov ).

En voici un exemple simple, dû à R. Thom: il s’agit du difféomorphisme f du tore T² = R²/Z² (produit de deux cercles de longueur 1) défini à partir de l’automorphisme linéaire unimodulaire de R² de matrice:

Remarquons que, si a/b , c/d est un point de R² à coordonnées rationnelles, le point x correspondant sur T² est périodique pour f ; en effet, l’image par A ⁿ de a/b , c/d est de forme ( 見/bd , 廓/bd ), 見 et 廓 entiers, et il n’y a qu’un nombre fini de tels points à des translations entières près. Les points périodiques de f sont donc denses dans T², ce qui montre que 行(f ) = T².

Les valeurs propres de A sont:

et les sous-espaces propres correspondants ont pour équation

(correspondant à^s ) et:

(correspondant à^u ); on voit immédiatement que T² est hyperbolique pour f : l’espace tangent en chaque point s’identifie canoniquement à R², E^{s (x 0,y 0)} est la droite:

et E^{u (x 0,y 0)} la droite:

respectivement parallèles aux sous-espaces propres de A . Pour la métrique usuelle de R², on peut prendre C = 1 et =^s = 1/^u . Les deux feuilletages (décomposition en union de courbes) définis dans R² par ces droites sont invariants par A (l’image d’une de ces droites par A est une droite parallèle); les pentes étant irrationnelles, ils se projettent dans T² en deux feuilletages transverses (ici orthogonaux) invariants par f (feuilletage stable et feuilletage instable) dont chaque feuille est dense dans T² (fig. 66, à interpréter comme les figures des premiers chapitres: R² est le déroulement (dans deux directions) de T²).

En chaque point du tore passent deux courbes (la feuille de chacun des feuilletages) respectivement tangentes au sous-espace stable et au sous-espace instable; la feuille stable (instable) de x est envoyée par f sur la feuille stable (instable) de f (x ) et la restriction de f à la feuille stable (instable) est une contraction uniforme (resp. une dilatation uniforme).

Cette description globale de la dynamique est permise ici par la linéarité de A mais, de même qu’en un point fixe hyperbolique d’un difféomorphisme non linéaire existent des sous-variétés invariantes respectivement tangentes au sous-espace propre stable et au sous-espace propre instable de la dérivée, on peut démontrer (Hirsh, Pugh, Shub) l’existence de telles variétés stables et instables pour un ensemble hyperbolique général (avec une démonstration du même type bien que plus compliquée): pour toute partie X de M, notons:

où d désigne la distance sur M définie par la métrique riemannienne. Si 炙 est un fermé invariant hyperbolique pour f , W^s (x ) et W^u (x ) sont des sous-variétés injectivement immergées pour tout point x de 炙, tangentes en x respectivement à E^{s x} et E^{u x} . Si 﨎 礪 0 est assez petit,

sont des disques plongés respectivement dans W^s (x ) et W^u (x ) et varient continûment avec x dans 炙. La figure 67 représente ces feuilletages invariants dans le cas très simple du solénoïde. Bien entendu, W^s (x ) et W^s (y ) ne peuvent se couper si x est différent de y ; de même pour W^u (x ) et W^u (y ); ainsi, lorsque f est un difféomorphisme d’Anosov, on retrouve comme dans l’exemple précédent deux feuilletages (décompositions de la variété en union de sous-variétés «parallèles») transverses invariants, le feuilletage stable sur lequel f contracte et le feuilletage instable sur lequel f dilate.

Les propriétés très fortes de stabilité associées aux ensembles invariants hyperboliques découlent toutes du théorème fondamental suivant: Tout morphisme approché d’un système dynamique dans un système dynamique hyperbolique est proche d’un unique vrai morphisme . Plus précisément, si 炙 est un ensemble hyperbolique pour f , il existe des nombres positifs 見, K, et un voisinage U size=1^炙 de 炙 dans M, tels que, pour tout homéomorphisme h d’un espace topologique X dans lui-même et toute application continue i de X dans U size=1^炙 telle que le diagramme:

commute à 見 près (d (f 獵 i , i 獵 h ) 諒 見), il existe une unique application continue j de X dans M telle que d (i , j ) 諒 Kd (f 獵 i , i 獵 h ) rendant commutatif le diagramme:

c’est-à-dire f 獵 j = j 獵 h .

La démonstration se fait en étudiant l’application k 料 f 獵 k 獵 h ^-1 de l’espace des applications continues k : XM dans lui-même: au voisinage de i , cette application est proche d’un opérateur linéaire hyperbolique et possède donc un unique point fixe qui est l’application j cherchée. On voit d’ailleurs que rien ne change si on remplace f par un difféomorphisme g assez proche de f dans la topologie C¹.

Si on choisit X = Z, h (i ) = i + 1, on obtient la propriété d’approximation de toute 﨎-pseudo-orbite (au sens défini plus haut) par une vraie orbite; en choisissant X = M, h proche de f , on peut en déduire la stabilité structurelle des difféomorphismes d’Anosov: si g est assez proche du difféomorphisme d’Anosov f dans la topologie C¹, il existe un homéomorphisme j proche de l’identité, qui conjugue g à f (g = j ^-1 獵 f 獵 j ). Le comportement dynamique complexe de ces difféomorphismes résiste donc à de petites perturbations. Bien entendu, des théorèmes analogues valent pour les flots.

Introduisons enfin l’axiome A de Smale: un difféomorphisme f sera dit axiome A si 行(f ) est hyperbolique et si de plus 行(f ) = Per(f ). (c’est-à-dire si tout point de 行(f ) est limite de points périodiques). Dans le cas des flots, on demande que 行( 﨏^t ) = 行 聆 F (union disjointe), où F est l’ensemble (fini) des équilibres, tous hyperboliques, et 行 est hyperbolique et coïncide avec l’adhérence de l’ensemble des orbites périodiques du flot.

Ces conditions sont assez fortes pour qu’il existe un bon théorème de structure, dû à Smale, le théorème de décomposition spectrale , assurant (dans le cas d’un difféomorphisme) l’existence d’une décomposition 行(f ) = 行¹ 聆 ... 聆 行^s de 行(f ) en ensembles hyperboliques basiques 行ⁱ , c’est-à-dire en fermés invariants disjoints ayant les propriétés suivantes: les points périodiques de la restriction f | 行ⁱ sont denses dans 行ⁱ , f | 行ⁱ est topologiquement transitif ( 行: est l’adhérence d’une orbite), et il existe un voisinage U de 行ⁱ tel que

Le nom du théorème vient du cas particulier du difféomorphisme induit sur l’espace projectif complexe P^{n -1}(C) par un élément du groupe linéaire SL(n , C). Dans la démonstration, on explicite les 行ⁱ comme adhérences des classes d’équivalence de points périodiques pour la relation d’équivalence p 黎 q si les intersections W^u (orbite de p ) 惡 W^s (orbite de q ) et W^s (orbite de p ) 惡 W^u (orbite de q ) sont non vides et transverses.

Si Per(f ) est hyperbolique et si une suite de points périodiques p ⁱ converge vers x , p ⁱ est équivalent à p ^j pour i et j assez grands à cause de la transversalité de W^u size=1^﨎(x ) et W^s size=1^﨎(x ) et de la continuité en x des variétés stables et instables: il n’y a donc qu’un nombre fini de classes d’équivalence.

Dans le même ordre d’idées, Per(f ) possède une structure de produit local dès qu’il est hyperbolique: cela signifie que si 﨎 et 嗀 sont assez petits, si x et y appartiennent à Per(f ) et sont distants de moins de 嗀, l’intersection de W^s size=1^﨎(x ) et W^u size=1^﨎(y ) est un unique point [x , y ] qui appartient à Per(f ) (fig. 68).

Cette structure de produit local permet de préciser le théorème d’approximation des 﨎-pseudo-orbites par de vraies orbites:

Pour tout 兀 礪 0, il existe un 﨎 礪 0 tel que toute 﨎-pseudo-orbite (x ⁱ ) dans Per(f ) soit « 兀-pistée» par l’orbite d’un point x de Per(f ) (au sens où d (f ⁱ (x ), x ⁱ ) est inférieur à 兀 pour tout i ); un tel x est forcément unique si 兀 est assez petit (expansivité de f |Per(f )). Ce théorème est l’outil fondamental d’étude de la dynamique d’un difféomorphisme axiome A; il implique par exemple que, pour chaque ensemble basique 行ⁱ , on ait

ce qui donne un bon contrôle de la dynamique dans M tout entière en fonction de celle dans 行(f ) puisque

Il suffirait d’ailleurs d’étudier les attracteurs (ensembles basiques 行ⁱ tels que W^s ( 行ⁱ ) soit ouvert): on voit en effet que presque tout point de M appartient au bassin d’un attracteur. En fait, la dynamique de f sur chacun des ensembles basiques 行ⁱ est très bien connue grâce aux méthodes de la dynamique symbolique, déjà rencontrées à propos du fer à cheval et du solénoïde.

La possibilité de coder les orbites de la restriction de f à 行ⁱ vient de l’existence de partitions de Markov de 行ⁱ : découvertes par Adler et Weiss pour les automorphismes linéaires du tore, par Sinaï pour les difféomorphismes d’Anosov, et par Bowen pour les ensembles basiques généraux, ce sont des décompositions de 行ⁱ en union de Q sous-ensembles R^j , j = 0, 1, ..., Q 漣 1, les rectangles de la partition (ainsi nommés car contenant [x , y ] dès qu’ils contiennent x et y , (fig. 68)), ayant par rapport à f de «bonnes» propriétés que nous n’expliciterons pas.

L’exemple le plus simple est la décomposition de l’ensemble de Cantor K du fer à cheval du chapitre 4 en (K 惡 ⁺⁰) 聆 (K 惡 ⁺¹).

Exactement comme au chapitre 4 (ou au chap. 6 pour le solénoïde) on définit une application surjective 崙: ^A行ⁱ en associant à une suite ( 見ⁿ )ⁿ size=1^{捻 Z} d’entiers égaux à 0, 1, ..., ou Q 漣 1, l’unique point x de l’intersection

c’est-à-dire l’unique point x dont l’orbite ait l’«histoire» décrite par la suite ( 見ⁿ )ⁿ size=1^{捻 Z} (x 捻 R size=1^見0, f (x ) 捻 R size=1^見1, f ^-1(x ) 捻 R size=1^見-1, etc.). ^A est maintenant le sous-espace de0, 1, ..., Q 漣 1^Z formé des suites ( 見ⁿ )ⁿ size=1^{捻 Z} telles que A size=1^{見n +1} size=1^見n = 1 pour tout n dans Z, où A = (A ^jk ) est la matrice Q 憐 Q définie par A ^jk = 1 si l’intérieur de R^k 惡 f ^-1(R^j ) est non vide dans 行ⁱ , A ^jk = 0 autrement (matrice de transition ). L’application 崙 vérifie 崙 獵 靖^A = f 獵 崙, où 靖^A : ^AA est le décalage de virgule défini comme au chapitre 4 (le système dynamique défini par 靖^A est appelé un sous-shift de type fini ). De plus, l’image réciproque 崙^-1(x ) d’un point de 行ⁱ n’a jamais plus de Q² éléments et est réduite à un seul élément pour la plupart des points de 行ⁱ , ce qui permet de déduire de l’étude de 靖^A beaucoup de renseignements sur la dynamique de f| 行ⁱ . Par exemple, 見 捻 ^A est un point périodique de 靖^A si et seulement si 崙( 見) 捻 行ⁱ est un point périodique de f , ce qui a permis à Manning de prouver la rationalité de la fonction 﨣^{f |} size=1^行i de f | 行ⁱ , où:

et Nⁿ (f ) est le nombre de points périodiques de période n de f . De même, les entropies topologiques h ( 靖^A ) et h (f | 行ⁱ ) sont égales, ce qui permet de prouver l’entropy conjecture pour f | 行ⁱ (Shub, Williams).

De nombreux attracteurs étranges vérifiant l’axiome A ont été étudiés: en plus du solénoïde (chap. 6), citons les difféomorphismes D.A. (dérivés d’Anosov) de Smale obtenus à partir d’un automorphisme linéaire du tore en «gonflant» la variété instable de (0, 0) suivant un procédé rappelant la construction de «contre-exemples de Denjoy» dans la théorie des difféomorphismes du cercle (fig. 69, que l’on comparera à la figure 66); l’attracteur est ici:

où 嗀 est un petit disque centré en (0, 0) et ² est la réunion de l’attracteur (qui coupe localement les feuilles du feuilletage stable suivant un ensemble de Cantor) et de la variété instable de (0, 0), qui est un ouvert dense dont on a représenté une partie en grisé. Citons également les attracteurs de R² de la forme:

(cf. fig 70, où l’on a représenté l’image par f du disque troué D), introduits par Plykin, qui sont les premiers exemples d’attracteurs hyperboliques dans R² (ce que n’est sûrement pas l’attracteur de Hénon s’il existe!).

On a représenté sur les figures 69 et 70 le feuilletage stable, invariant par f et transverse à l’attracteur; la restriction de f à ce feuilletage est une contraction alors que, dans les deux cas, sa restriction à l’attracteur est une dilatation, ce qui permet d’appliquer la théorie de Williams décrite à propos du solénoïde: l’attracteur est homéomorphe à la limite projective 良lim (K 﨏 良 K 﨏 良 K ...), où K est la «variété branchue» de dimension 1 obtenue en faisant le quotient d’un voisinage de l’attracteur (respectivement le complémentaire d’un petit disque 嗀 centré sur (0, 0) et un disque D percé de trois trous) par le feuilletage stable (on identifie à un point chaque feuille) et 﨏 l’application dilatante induite par f (voir les figures; sur la figure 69, on remarquera que, en homologie, la matrice de 﨏 est (2 11 1), ce qui n’est pas un hasard). Rappelons également que la dynamique sur l’attracteur est l’extension naturelle de 﨏 à la limite projective.

Un sous-produit de la démonstration du théorème de décomposition spectrale est le fait qu’un difféomorphisme d’Anosov vérifie l’axiome A. Nous avons indiqué plus haut la «stabilité structurelle» de ces difféomorphismes; en ce qui concerne les difféomorphismes axiome A généraux, deux théorèmes fondamentaux de stabilité existent: disons que f possède un cycle s’il existe i ⁰, ..., i ^t compris entre 1 et s tels qu’il existe une orbite issue (asymptotiquement) de 行^{i 0} et aboutissant asymptotiquement à 行^{i t} , une orbite issue de 行^{i t} et aboutissant à 行^{i t -1}, ..., une orbite issue de 行^{i 1} et aboutissant à 行^{i 0}, c’est-à-dire tels que:

Le théorème de 行-stabilité de Smale affirme qu’un difféomorphisme f axiome A qui n’a pas de cycle est 行-stable (si g est proche de f dans la topologie C¹, f | 行(f ) et g | 行(g ) sont conjugués par un homéomorphisme de 行(f ) sur 行(g )).

L’absence de cycle permet de construire une filtration M¹ 說 M² 說 ... 說 M^s de M par des sous-variétés à bords invariantes telles que:

ce qui ramène essentiellement l’étude à celle de la stabilité de f | 行ⁱ pour chaque i . Cette dernière se montre comme pour les difféomorphismes d’Anosov en utilisant la propriété d’approximation des 﨎-orbites.

Enfin, comme beaucoup de théorèmes analogues ^{[cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES (la théorie mathématique et ses applications)]}, le théorème de stabilité structurelle de Robbin-Robinson est fondé sur une propriété de transversalité. Rappelons que deux sous-variétés sont transverses si en chaque point d’intersection leurs espaces tangents engendrent l’espace tangent à la variété ambiante. Le théorème affirme qu’un difféomorphisme f axiome A qui vérifie la condition de transversalité forte (W^s (x ) et W^u (y ) transverses quels que soient x et y dans 行(f )) est structurellement stable: si g est C¹-proche de f , il lui est conjugué par un homéomorphisme.

Réciproquement, il est possible de montrer que certaines formes très fortes de stabilité impliquent l’axiome A, mais on ne sait pas si la stabilité structurelle suffit à l’impliquer.

Signalons que, contrairement à ce qui se passe pour les orbites périodiques (théorème de Kupka-Smale), la non-transversalité des variétés stables et instables d’un ensemble invariant hyperbolique 炙 peut être un phénomène persistant, même en dimension deux: à l’aide de la notion d’épaisseur d’un ensemble de Cantor, Newhouse a montré qu’un difféomorphisme de classe C² d’une surface peut présenter d’une manière stable (insensible aux perturbations) la configuration de la figure 71, qui représente un ensemble invariant hyperbolique «sauvage» 炙, avec un ensemble de Cantor de points de tangence. Un tel difféomorphisme n’est pas structurellement stable (il n’est même pas 行-stable) et il en est de même des difféomorphismes suffisamment voisins, ce qui montre que la stabilité n’est pas un phénomène aussi courant qu’on aurait pu l’espérer. Ces exemples sont importants car, sous certaines conditions, l’ensemble hyperbolique 炙 en question peut être limite d’une infinité d’orbites périodiques attractantes de périodes arbitrairement grandes, porte ouverte à d’éventuelles confusions avec un attracteur étrange (exemple de l’attracteur de Henon?)

Nous terminons ce chapitre par une courte discussion de l’exemple qui est à la base de la théorie des systèmes dynamiques hyperboliques, le flot géodésique sur une surface riemannienne compacte de courbure négative constante égale à 漣 1. Une telle surface s’écrit comme le quotient du demi-plan H =z = x + iy 捻 C|y 礪 0 muni de la métrique «hyperbolique»:

par un sous-groupe discret 臨 du groupe PSL(2, R) des isométries de H conservant l’orientation; H est le plan de Lobatchevski , ou demi-plan de Poincaré et les éléments de PSL(2, R) sont représentés par des matrices:

à coefficients réels, de déterminant égal à 1, agissant sur H par:

Le flot géodésique ( 﨏^t )^t size=1^{捻 R} est défini sur la variété de dimension trois ¹ dont les points sont les vecteurs de longueur 1 tangents à la surface . Si X 捻 ¹ est un tel vecteur tangent, 﨏^t (X) est le vecteur tangent obtenu en suivant sur la longueur t dans la direction de X l’unique géodésique tangente à X (fig. 72).

Une courbe intégrale du flot correspond donc à une géodésique orientée de (les géodésiques sont solution d’une équation différentielle du second ordre; par chaque point il en passe une infinité et le passage de à ¹ , analogue au passage de (N¹) à (P¹) au chapitre 1, permet de les disjoindre deux à deux).

La variété ¹ étant le quotient de ¹H par l’action naturelle de 臨, la description locale du flot est la même pour ¹ et ¹H.

Les géodésiques de H sont les demi-droites perpendiculaires à l’axe des x et les demi-cercles centrés sur l’axe des x . Étant donné une géodésique orientée 塚, l’ensemble des géodésiques orientées qui lui sont positivement (resp. négativement) asymptotiques forme dans ¹H une sous-variété de dimension deux, la variété stable (resp. instable) de la courbe intégrale 塚 correspondant à 塚 (fig. 73 a dans H et fig. 73 b dans ¹H). Chacune de ces variétés est feuilletée par des courbes transportées les unes sur les autres par le flot, les variétés stables (resp. instables) des points de 塚; ces courbes sont en fait les courbes intégrales de deux autres flots sur ¹H (ou T¹ ), les flots horocycliques , définis comme le flot géodésique mais à partir des horocycles orientés (cercles tangents à l’axe des x , orthogonaux aux géodésiques) (fig. 73).

Pour donner une idée de la richesse de ces flots sur ¹ , indiquons que chaque courbe intégrale d’un flot horocyclique est dense dans ¹ (Hedlund, 1936) et que le flot géodésique sur ¹ est ergodique (Hopf, 1939). Une description algébrique très simple peut en être donnée en remarquant que ¹H s’identifie naturellement à PSL(2, R) (étant donné un vecteur X de longueur 1 tangent à H, il existe un unique élément de PSL(2, R) transportant sur X le vecteur (0, 1) tangent au point i ). ¹ s’identifie donc à l’espace homogène PSL(2, R)/ 臨; le flot géodésique est défini par:

et les flots horocycliques par:

Les problèmes relatifs à ces flots deviennent alors très algébriques: par exemple, Ratner a montré que le simple isomorphisme mesurable des flots horocycliques sur PSL(2, R)/ 臨¹ et PSL(2, R)/ 臨² suffit à impliquer l’existence de g 捻 PSL(2, R) tel que 臨² = g ^-1 臨¹g .

Remarquons qu’en intervertissant l’ordre dans lequel on fait agir R (via multiplication à gauche par les matrices: (e ^t 00 e ^-t )) et 臨 sur PSL²(R), on est amené à étudier l’action de 臨 sur l’«espace» des géodésiques de H; en remplaçant H par l’intérieur du disque unité, on voit que cet «espace» des géodésiques s’identifie naturellement au complémentaire de la diagonale dans le tore S¹ 憐 S¹ (une géodésique est déterminée par ses deux points à l’infini). Les propriétés du flot géodésique sur = H/ 臨 sont reflétées dans celles de l’action «diagonale» de 臨 sur S¹ 憐 S¹. On rejoint ainsi la théorie des groupes fuchsiens , développée à partir de 1882 par Poincaré à propos de ses recherches sur les équations différentielles linéaires du second ordre dans le domaine complexe, et très étudiée aujourd’hui (en particulier du point de vue dynamique, par Sullivan dans une suite de travaux fascinants); dans cette optique, la théorie des systèmes dynamiques devient en fait la théorie des actions différentiables des groupes sur les variétés et pas seulement celle des actions de Z ou de R.

8. Mesures invariantes et exposants de Liapunov

Nous abordons dans ce chapitre quelques aspects de la théorie ergodique des systèmes dynamiques qui, sous l’impulsion initiale de Birkhoff et von Neumann (démonstration en 1931 du théorème ergodique) puis de l’école russe (introduction de l’entropie métrique par Kolmogorov en 1958, travaux de Sinaï sur l’hypothèse ergodique à partir de 1962), s’est beaucoup développée ces dernières années: lorsque la description géométrique d’un portrait de phase s’avère trop difficile, on la remplace par une description de type statistique; comme en thermodynamique, on cherche des «mesures d’équilibre» qui décrivent le comportement asymptotique de «presque toutes» les courbes intégrales (orbites) du système. Dans le cas hamiltonien, il existe une mesure invariante lisse, la mesure de Liouville; en 1885, Boltzmann émet sa fameuse hypothèse à propos de la théorie cinétique des gaz: en termes d’aujourd’hui, il suppose que la situation «générale» pour un système hamiltonien est l’ergodicité dans les hypersurfaces E d’énergie constante, c’est-à-dire que tout sous-ensemble mesurable A invariant par le flot et situé dans une hypersurface E vérifie(A) = 0 ou(E 漣 A) = 0 (il supposait en fait que chaque orbite pouvait parcourir toute l’hypersurface E, ce qui est impossible). Le théorème de K.A.M. est venu infirmer cette hypothèse en montrant l’existence de familles de tores invariants de mesure totale positive pour les systèmes hamiltoniens suffisamment proches (dans la topologie C^k , k assez grand) d’un système complètement intégrable.

Cependant, pour certains systèmes dynamiques possédant une mesure invariante lisse, l’ergodicité a pu être démontrée. On a cité à la fin du chapitre 7 le flot géodésique sur une surface de courbure négative. Le billard sur un tore troué (fig. 74 a), dans lequel la boule se déplace sur le tore suivant des droites (après déroulement sur R²) et se réfléchit sur le bord du trou avec un angle de réflexion égal à l’angle d’incidence (fig. 74 b: tennis dans une forêt), est un cas limite d’un tel flot: si on imagine que la trajectoire de la boule change de face chaque fois qu’il y a réflexion, on est amené à remplacer le tore troué par son double, isomorphe à une surface de genre deux (fig. 74 c). Ce billard est un modèle du système composé de deux particules se déplaçant sur un tore et ayant des chocs élastiques et son ergodicité pour la mesure invariante naturelle a été prouvée par Sinaï. (C’est Kolmogorov qui, après avoir montré que l’hypothèse ergodique n’était pas vérifiée à cause des orbites périodiques elliptiques – théorème K.A.M. –, aurait suggéré qu’elle pouvait être rétablie par la prise en compte des chocs.) Le flot géodésique sur un ellipsoïde est au contraire «complètement intégrable»; le billard correspondant est un billard elliptique , lui aussi complètement intégrable dans le sens où il existe une famille à un paramètre de caustiques enveloppant les trajectoires (théorème de Poncelet, fig. 75, où l’on a représenté une telle caustique). Si on déforme un peu le bord, le théorème de K.A.M. montre que cette famille à un paramètre est remplacée par un ensemble de Cantor de caustiques. Par contre, de grandes déformations conduisent aux stades dont l’ergodicité a pu être démontrée par Bunimovitch (fig. 76, où l’on a représenté un morceau de trajectoire du stade).

Dans le cas dissipatif, on s’intéresse essentiellement aux mesures invariantes supportées par les «attracteurs» du système, c’est-à-dire par celles qui décrivent la statistique de «suffisamment» d’orbites.

Il est bien naturel de commencer par des moyennes de mesures de Dirac 嗀^y le long des orbites, c’est-à-dire (dans le cas d’un difféomorphisme f auquel nous nous limiterons) par les «limites vagues» de suites de mesures de la forme:

où x est un point de l’espace de phase. Si nous ne considérons que des orbites bornées, dire que la limite vague 猪^x existe signifie que pour toute fonction 﨏 continue sur l’espace de phase, la limite

existe.

Lorsqu’elle existe, une telle mesure asymptotique est nécessairement invariante ( 猪^x ( 﨏 獵 f ) = 猪^x ( 﨏)); d’autre part, la masse totale 猪^x (1) est égale à 1 (on dit que 猪^x est une mesure de probabilité ). Ces deux propriétés impliquent, grâce au théorème de récurrence de Poincaré (cf. théorie ERGODIQUE) que le support de 猪^x est contenu dans l’ensemble non errant 行(f ) de f (et donc que 猪^x (X 漣 行(f )) = 0). Dorénavant, nous supposerons que 行(f ) est compact. Pour ce qui concerne l’existence de ces mesures 猪^x , tournons-nous vers le théorème ergodique de Birkhoff (cf. théorie ERGODIQUE): si f : XX laisse invariante une mesure de probabilité 猪, il existe un sous-ensemble Y de X de 猪-mesure nulle tel que 猪^x existe pour tout point x de X 漣 Y; si de plus (f , 猪) est ergodique (c’est-à-dire que les seuls ensembles mesurables invariants par f sont de mesure 0 ou 1), on peut choisir Y de façon que, pour tout x dans X 漣 Y, 猪^x soit égale à 猪. Remarquons qu’a priori cette affirmation ne dit rien sur les points errants puisque, comme nous venons de le voir, X 漣 行(f ) est de 猪-mesure nulle.

D’un autre côté, nous voyons que, quelle que soit la mesure de probabilité ergodique invariante par f , il existe un point x dans X tel que 猪^x soit définie et égale à 猪; ce résultat prend tout son sel quand on se rend compte que, bien souvent, il existe énormément de mesures de probabilité ergodiques invariantes; dans le cas de l’homéomorphisme 靖 : face=F00190, 1^Z索0, 1^Z de décalage (shift) sur deux symboles introduit à propos du fer à cheval ou du solénoïde, on en construit très facilement une infinité, paramétrées par p , 0 麗 p 麗 1, en donnant la mesure p size=1^見(1 漣 p )^k size=1^漣 size=1^見 au sous-ensemble de0, 1^Z formé des suites (a ^l )^l size=1^{捻 Z} telles que a ^{i 1} = a ^{0i 1}, ..., a ^{i k} = a ^{0i k} , où 見 des a ^{0i j} sont des 0 (jeu de pile ou face biaisé encore appelé schéma de Bernoulli en souvenir des travaux de J. Bernoulli sur ce jeu). Toutes ces mesures sont bien sûr étrangères l’une à l’autre. Cependant, les expériences numériques suggèrent assez souvent qu’une mesure asymptotique bien définie puisse être associée à un «attracteur».

Dans le cas des difféomorphismes axiome A, l’explication est donnée par le théorème de Ruelle (1975, étendu par Ruelle et Bowen aux flots axiome A). Rappelons que les attracteurs 行ⁱ d’un difféomorphisme axiome A sont caractérisés parmi les ensembles basiques par une variété stable W^s ( 行ⁱ ) (bassin) de mesure de Lebesgue positive, ce qui implique que, pour presque tout (au sens de la mesure de Lebesgue) point x de X, l’orbite f ⁿ (x ) tend vers un attracteur lorsque n tend vers + 秊. Le théorème de Ruelle affirme qu’il existe dans ces conditions un sous-ensemble Zⁱ de W^s ( 行ⁱ ), de mesure de Lebesgue nulle dans W^s ( 行ⁱ ), et une mesure de probabilité ergodique 猪⁺ size=1^行i sur 行ⁱ invariante par f telle que, pour tout x dans W^s ( 行ⁱ ) 漣 Zⁱ , 猪^x existe et soit égale à 猪⁺ size=1^行i . La mesure 猪⁺ size=1^行i ainsi canoniquement attachée à un attracteur axiome A peut être légitimement appelée mesure d’équilibre de cet attracteur. Le comportement statistique de presque tout x (au sens de la mesure de Lebesgue) est donc décrit, dans le cas d’un difféomorphisme axiome A, par la donnée d’un nombre fini de mesures ergodiques.

La mesure 猪⁺ size=1^行i a d’autres caractérisations, très intéressantes: c’est l’unique mesure ayant la propriété d’absolue continuité des mesures conditionnelles sur les variétés instables de l’attracteur par rapport à la mesure de Lebesgue sur ces variétés instables, propriété elle-même reliée au fait d’être limite de mesures stationnaires sous l’action de petites perturbations stochastiques (cette dernière propriété explique, d’après Ruelle, les résultats numériques); c’est aussi l’unique mesure qui maximise (c’est-à-dire rend égale à 0) l’expression:

différence entre l’entropie métrique de 猪 (cf. théorie ERGODIQUE) et la moyenne de la somme (avec multiplicité) des exposants caractéristiques positifs (définis ci-dessous; on peut d’ailleurs écrire alternativement cette moyenne sous la forme:

où J⁺(x ) est le jacobien de la restriction Df (x )|E^{u x} : E^{x u}E^{u f(x )} de la dérivée de f au sous-espace instable en x ). On comparera cette dernière caractérisation au principe variationnel de la mécanique statistique d’équilibre. Remarquons que, avec la fonction Log J⁺(x ), on mesure le taux suivant lequel on perd de l’information sur la position de f ⁿ (x ) quand on ne connaît pas exactement la position de x ; la non-nullité de l’entropie métrique h size=1^猪(f ) est donc une manifestation de la dépendance sensible des conditions initiales , cause essentielle du passage d’une description géométrique des mouvements asymptotiques à une description statistique.

Remarquons enfin qu’en général 猪⁺ size=1^行i ne maximise pas l’entropie métrique de f ; il existe une mesure 猪⁰ size=1^行i réalisant le maximum (Parry, Bowen), c’est-à-dire telle que h size=1^猪0 size=1^行i (f ) = h (f ) (entropie topologique); cette mesure régit la répartition asymptotique des points périodiques de f dans 行ⁱ .

Nous ne rentrerons pas dans la démonstration du théorème de Ruelle, qui peut se faire en utilisant la dynamique symbolique des ensembles basiques, mais nous donnerons plutôt une idée de la mesure dans l’exemple, décrit au chapitre 6, du difféomorphisme f de S¹ 憐 D² dans lui-même ayant pour attracteur le solénoïde A. La mesure de Lebesgue (de Haar) d /2 神 sur le cercle |z | = 1 est une mesure de probabilité invariante ergodique pour l’application 嗀 du cercle sur lui-même définie par 嗀(z ) = z ² (c’est-à-dire 嗀( ) = 2 si z = e ⁱ size=1). Le théorème de Birkhoff fournit alors un sous-ensemble Y de mesure nulle dans le cercle tel que, si 見: S¹R est continue,

pour tout z n’appartenant pas à Y. Notons Z l’ensemble des points (z , w ) du tore solide S¹ 憐 D² dont la première coordonnée z appartient à Y. La mesure de Lebesgue de Z dans S¹ 憐 D² (qu’on peut dérouler en R 憐 D² 說 R³) est nulle et il n’est pas difficile de montrer que si 﨏: S¹ 憐 D²R est continue et si (z , w ) n’appartient pas à Z,

où la mesure d’équilibre 猪^+A est définie de la manière suivante:

1. On calcule la moyenne de 﨏 (pour la mesure de Lebesgue d /2 神) sur chacun des brins du solénoïde A faisant une fois le tour du tore solide (morceaux de variétés instables), ce qui fournit une fonction 﨏 sur l’ensemble de Cantor K 黎=0, 1^Z intersection d’un disque z = C ^te avec l’attracteur.

2. On moyenne 﨏 sur K muni de la mesure produit correspondant à p = 1/2 décrite plus haut à propos du shift. Ce résultat devient plausible lorsqu’on remarque que la distance de f ⁿ (z , w ¹) et f ⁿ (z , w ²) tend vers 0 lorsque n tend vers + 秊 (invariance du feuilletage stable formé des disques z = C ^te , sur lequel f est contractante).

Si on perturbe un peu le difféomorphisme f , on conserve un attracteur homéomorphe au solénoïde (stabilité structurelle) mais la mesure d’équilibre change; il faut maintenant affecter de poids les deux opérations de moyenne.

Dans cet exemple, l’absolue continuité par rapport à la mesure de Lebesgue de la mesure d’équilibre sur les variétés instables de l’attracteur A est évidente, et le comportement du couple (f |A, 猪^+A) ne se distingue pas d’un jeu de pile ou face au niveau probabiliste (en particulier, h size=1^猪+A(f ) = Log2 et coïncide donc dans ce cas très simple avec l’entropie topologique h (f )). Ce dernier fait est général: dès que la restriction de f à un attracteur hyperbolique 炙 est topologiquement mélangeante c’est-à-dire si pour tout couple (U, V) d’ouverts de 炙, il existe un entier N 閭 0 tel que, si n 礪 N, on ait f ⁿ (U) 惡 V 禮), le couple (f | 炙, 猪⁺ size=1^炙) est isomorphe au niveau probabiliste à un schéma de Bernoulli défini comme précédemment en remplaçant éventuellement0, 1^Z par0, 1, ..., q 漣 1^Z et en affectant, dans la définition de la mesure produit, chacun des entiers i = 0, ..., q 漣 1 d’une probabilité p ⁱ , avec p ⁱ = 1.

Autrement dit, les propriétés probabilistes de (f | 炙, 猪⁺ size=1^炙) sont celles d’un jeu de pile ou face biaisé généralisé, nec plus ultra du comportement aléatoire.

La théorie ergodique des systèmes dynamiques ne vérifiant pas l’axiome A n’existe pas encore, mais une vaste généralisation de la notion d’hyperbolicité a été développée à partir de la notion d’exposant caractéristique (ou exposant de Liapunov , car introduite par ce dernier dans sa théorie de la stabilité). Soit f un difféomorphisme d’une variété compacte M, munie d’une métrique riemannienne, et 猪 une mesure de probabilité sur M invariante par f . Le théorème ergodique non commutatif d’Oseledets fournit un sous-ensemble B de M de mesure 猪(B) = 1, invariant par f , et pour chaque point x de B une décomposition:

de l’espace tangent en x à M en somme directe de sous-espaces ayant les propriétés suivantes:

1. 葉v 捻 Wⁱ (x ), les limites:

existent, sont égales et indépendantes de v (Df ^{x n} est la dérivée en x de f ⁿ ). On les noteⁱ (x ) avec la convention que¹(x ) 麗²(x ) 麗 ... 麗^{s (x )}(x ). Lesⁱ (x ) sont les exposants caractéristiques de f au point x .

2. Lesⁱ et les Wⁱ dépendent «mesurablement» de x (mais pas continûment a priori).

3. Pour tout point x de B, s (f (x )) = s (x ), ⁱ (f (x )) =ⁱ (x ) et Wⁱ (f (x )) = Df ^x (Wⁱ (x )).

Il est naturel de noter:

et:

et de les appeler respectivement sous-espace stable et sous-espace instable en x .

Vaste généralisation du théorème de décomposition spectrale d’une matrice, ce théorème a son pendant non linéaire (existence de variétés stables pour f ) démontré par Piesin en 1977: pour tout point x de B tel que E^s (x ) ne soit pas réduit à0, l’ensemble:

est une sous-variété immergée dans M, admettant E^s (x ) pour espace tangent en x , contenue dans B, et égale en fait à:

La démonstration est une application sophistiquée de la méthode des «transformées de graphes» qui sert à établir les théorèmes d’existence de variétés stables dans le cas d’un point fixe ou d’un ensemble hyperbolique.

Lorsque la mesure invariante 猪 est définie à partir de la mesure de Lebesgue par une densité différentiable, les feuilletages stables et instables ont la propriété d’absolue continuité des mesures conditionnelles par rapport à la mesure de Lebesgue des feuilles; Piesin en déduit la formule de l’entropie (démontrée par ailleurs directement par Mañe):

qui, ainsi qu’il fut déjà remarqué, rend manifeste le lien entre la positivité de h size=1^猪(f ) et la dépendance sensible des conditions initiales évidemment liée à la présence d’exposantsⁱ (x ) positifs.

Enfin, dans le cas où aucun exposant n’est nul (hyperbolicité généralisée), le difféomorphisme f ressemble «presque partout» (au sens de la mesure invariante 猪) à un difféomorphisme d’Anosov non uniforme (les rapports de dilatation et de contraction varient de façon non contrôlée mais, en oubliant un sous-ensemble de petite mesure, on peut supposer les feuilletages stables et instables continus).

Il est alors possible de montrer un théorème de structure, cousin ergodique du théorème de décomposition spectrale de Smale: le sous-ensemble B des théorèmes d’Oseledets et de Piesin se décompose (à un sous-ensemble de mesure nulle près) en parties (boréliennes) disjointes Bⁱ invariantes par f , éventuellement en nombre infini, telles que (f |Bⁱ , 猪|Bⁱ ) soit ergodique pour tout i (composantes ergodiques, dont la dynamique peut se décrire à l’aide de schémas de Bernoulli (à un isomorphisme mesurable près)).

Il est vraisemblable que toute cette théorie s’applique aux équations de Navier-Stokes et donc au problème de la turbulence hydrodynamique ^{[cf. TURBULENCE]}: si on admet que ces équations définissent effectivement un système dynamique dans un certain ouvert 六 de l’espace des champs de vitesses de divergence nulle définis dans un domaine 阮 de R³ et nuls au bord de 阮, muni de la norme de Sobolev H¹ (ce qui semble être le cas dans les régimes faiblement turbulents), tout sous-ensemble attractant 炙 contenu dans 六 est compact, de dimension de Hausdorff finie (Mallet-Paret) et admet un plongement injectif dans un espace de dimension finie (Mañe). On peut alors définir des exposants caractéristiques qui forment une suite tendant vers 漣 秊. Le lecteur trouvera des exposés de synthèse sur cette question dans les travaux de Ruelle et se reportera aux belles expériences de Gollub et Swinney, et de Libchaber.

Au terme de cette excursion dans la théorie ergodique, signalons que le théorème d’Oseledets et la notion d’exposant caractéristique ont un sens pour des «applications fibrées» (skew products , ou produits gauches ) F = f | 憐 A: X 憐 RⁿX 憐 Rⁿ de la forme F(x , v) = (f (x ), A(x )v), où A est une application de X dans GL(n , R), groupe des isomorphismes linéaires de Rⁿ . Lorsque X est une variété, de fibré tangent X 憐 Rⁿ , et que A(x ) est la dérivée de f en x , on retrouve la situation précédente. Dans le cas général, on voit qu’il faut remplacer Df ^{x n} = Df (f ^{n -1}(x )) 獵 ... 獵 Df (x ) par A^{f n} (x ) = A(f ^{n -1}(x )) 獵 A(f ^{n -2}(x )) 獵 ... 獵 A(x ). Si, par exemple, f : face=F00190, 1, ..., n^Z索0, 1, ..., n^Z est l’homéomorphisme de décalage, on retrouve la théorie des produits aléatoires de matrices. On généralise de même la notion de structure hyperbolique (décomposition du fibré trivial X 憐 RⁿX en somme directe continue de deux sous-fibrés invariants, l’un uniformément contractant, l’autre uniformément dilatant). Un théorème récent d’Herman fournit, à partir de la théorie des fonctions pluri-sous-harmoniques, des minorations d’exposants dans ce cadre. Par exemple, si F: S¹ 憐 R²S¹ 憐 R² est défini par f (e ² size=1^神i size=1, v) = (e ² size=1^{神i (} size=1⁺ size=1^見), A( )v), où:

et:

on obtient une minoration du plus grand exposant par Log(/2 + 1/(2)).

Notons que les matrices A( ) ne sont hyperboliques que si:

elliptiques sinon, et qu’il s’agit donc d’une lutte entre comportement hyperbolique et comportement elliptique. Notons également que pour des raisons évidentes d’homotopie [ 料 A( ) est un générateur du groupe de Poincaré de SL(2, R)], F ne possède pas de structure hyperbolique bien qu’ayant au niveau probabiliste un comportement hyperbolique pour 礪 1.

Des difféomorphismes fibrés apparaissent naturellement dans de nombreux problèmes; en particulier, la structure du spectre de l’équation de Hill discrétisée est intimement liée à l’existence de structures hyperboliques ou d’exposants caractéristiques pour de tels difféomorphismes.

Que la tension elliptique-hyperbolique soit un phénomène subtil est attesté par l’étude au voisinage d’un point fixe elliptique d’un difféomorphisme local du plan préservant les aires: on a vu au chapitre 4 que le comportement elliptique se produit sur un ensemble de mesure positive (courbes invariantes données par le théorème de K.A.M.), mais on ne sait pas s’il en est de même du comportement hyperbolique; il existe bien sûr en général une infinité de points homoclines, ce qui assure la positivité de l’entropie topologique (existence de fers à cheval), mais tout cela pourrait ne concerner qu’un sous-ensemble du plan de mesure nulle auquel cas l’entropie métrique serait nulle. Pour le moment personne n’en sait rien.

9. Propriétés génériques

Que dire d’un système dynamique «général»? Voilà une question bien vague...

Certains théorèmes sont universellement valables et limitent en quelque sorte l’«espace» de tous les systèmes dynamiques (théorèmes d’existence et d’unicité locales de Cauchy-Kowalewska, théorème ergodique de Birkhoff, théorème de Piesin...) mais, de même que le lemme de Sard pour la théorie des fonctions différentiables, ces théorèmes contraignent peu la géométrie du système.

Peut-on alors répéter ce qui a si bien réussi pour les fonctions ^{[cf. SINGULARITÉS DES FONCTIONS DIFFÉRENTIABLES (la théorie mathématique et ses applications)]}, c’est-à-dire définir une classe de systèmes dynamiques assez particulière pour que la dynamique des éléments de cette classe soit géométriquement descriptible, et assez générale pour contenir «presque tous» les systèmes dynamiques (ou «presque tous» les systèmes d’un certain type: gradient, hamiltonien, avec symétries, etc.)? Et si les éléments de cette classe étaient structurellement stables, quel soulagement pour les physiciens, ingénieurs, et autres chercheurs de modèles! Plus de problème métaphysique sur l’imprécision des prédictions due à l’imprécision des coefficients numériques s’il y a toutes les chances pour que l’équation écrite soit stable et si de plus la dynamique d’une telle équation admet une description géométrique simple.

Ce rêve, un moment caressé par Smale (les succès obtenus par Peixoto en dimension deux y incitaient), s’est avéré doublement naïf, mais sa poursuite est à l’origine de nombreuses découvertes passionnantes. Naïve en effet l’idée que la stabilité structurelle d’un système implique la possibilité de prédictions précises: pour les systèmes qui, tels les flots géodésiques ergodiques de la fin du chapitre 7, présentent une «dépendance sensible des conditions initiales», seule une prédiction de type statistique est possible (chap. 8). Naïf également l’espoir que «presque tout» système dynamique soit structurellement stable, ou même descriptible géométriquement. Et puis, quel sens donner à «presque tout»? Déjà dans sa belle adresse au congrès international des mathématiciens de 1954, Kolmogorov soulevait ce problème, opposant la notion topologique (catégorie de Baire) à la notion métrique (mesure), la première bien adaptée à la mise en avant des situations hyperboliques, la seconde tournée plutôt vers l’elliptique mais souffrant de la non-existence d’une mesure naturelle sur les espaces fonctionnels.

Plus précisément, soit 﨑^r (M) l’espace, muni de la topologie C^r , des champs de vecteurs r fois continûment différentiables sur la variété compacte M: deux champs sont considérés comme proches s’ils le sont en chaque point ainsi que leurs dérivées jusqu’à l’ordre r (dans une carte locale, un champ est simplement une application de Rⁿ dans Rⁿ si n est la dimension de M).

La stabilité structurelle a été définie à la fin du chapitre 2; on a introduit une relation d’équivalence sur 﨑^r (M) (existence d’un homéorphisme de M envoyant orbite orientée du premier champ sur orbite orientée du second) et appelé structurellement stable tout élément contenu dans l’intérieur de sa classe d’équivalence (c’est-à-dire dont la classe d’équivalence est ouverte): autrement dit, un champ qu’aucune déformation suffisamment petite ne fasse sortir de sa classe d’équivalence. Des définitions analogues valent pour l’espace Diff^r (M) des difféomorphismes r fois continûment différentiables de M, la relation d’équivalence étant cette fois l’existence d’un homéomorphisme conjuguant les deux difféomorphismes considérés (cf. chap. 7, la discussion de la stabilité des difféomorphismes d’Anosov).

Un sous-ensemble 阮 de 﨑^r (M) (resp. Diff^r (M)) est dit résiduel (au sens de Baire) s’il est intersection d’une famille au plus dénombrable d’ouverts denses: rappelons qu’un sous-ensemble D d’un espace topologique X est dit dense si tout point de X peut être approché arbitrairement près par un point de D. L’intersection de deux sous-ensembles denses peut très bien être vide en général (par exemple Q et R 漣 Q dans R) mais pas celle de deux ouverts denses; de plus, dans les espaces 﨑^r (M) et Diff^r (M), un sous-ensemble résiduel est dense (propriété de Baire ). Une propriété (P) concernant les champs de vecteurs (resp. les difféomorphismes) est dite générique (au sens topologique) si l’ensemble des éléments de 﨑^r (M) (resp. Diff^r (M) qui la possèdent est résiduel. Remarquons que la réunion d’une famille dénombrable de propriétés génériques est encore générique.

Le théorème de Peixoto (qui fait suite en 1962 aux travaux de Andronov et Pontriaguine, de Baggis, Pliss) affirme que, pour tout r 閭 1, les champs de vecteurs structurellement stables sur une variété orientable M de dimension deux (surface) forment un ouvert dense 崙 de 﨑^r (M).

La stabilité structurelle est donc une propriété générique dans ce cas. De plus, les champs appartenant à 崙 sont caractérisés par une structure géométrique très simple (champs de Morse-Smale): nombre fini de points singuliers, tous hyperboliques, nombre fini d’orbites périodiques, toutes hyperboliques (c’est-à-dire telles que l’application de Poincaré correspondante ait un point fixe hyperbolique), aucune connexion homocline entre cols, ensemble 諸-limite et 見-limite d’une orbite périodique toujours égal à un point singulier ou une orbite périodique (revoir les définitions données à la fin du chap. 3, première excursion topologique). Le cas particulier des champs de vecteurs sans singularité sur le tore qui possèdent une courbe fermée rencontrant transversalement chaque orbite correspond (via l’application de premier retour) à un théorème analogue pour les difféomorphismes du cercle et l’étude du nombre de rotation esquissée au chapitre 5 montre que, même dans ce cas qui semble favorable, il existe au moins deux notions de généricité antagonistes: si (X size=1^猪) size=1^猪 size=1^{捻 [0, 1]} est une famille à un paramètre de tels champs de vecteurs sur le tore, l’ensemble 琉 des valeurs du paramètre pour lesquelles X size=1^猪 n’est pas structurellement stable a une mesure de Lebesgue positive (et pouvant même être très proche de 1) dès que les nombres de rotation des applications de premier retour de X⁰ et X¹ différent; cela n’empêche pas que, si la famille (X size=1^猪) est générique au sens topologique dans l’espace des familles de ce type, le complémentaire de 溜 soit un ouvert dense. Ce genre d’opposition entre «presque tout topologique» et «presque tout métrique» est familier aux arithméticiens et aux analystes et je ne saurais mieux faire que recommander au lecteur le merveilleux petit livre d’Oxtoby: Measure and Category .

En dimension supérieure à deux (un pour les difféomorphismes), un certain nombre de propriétés génériques (au sens topologique) ont été dégagées, mais elles sont bien loin de contraindre la géométrie du système au point de la rendre descriptible. Citons en particulier la propriété de Kupka-Smale (1963) (points singuliers et orbites périodiques tous hyperboliques et transversalité des variétés stables et instables des points singuliers et orbites périodiques dans le cas d’un flot, même définition avec les points fixes et les points périodiques pour les difféomorphismes), tout en rappelant les résultats de Newhouse, cités au chapitre 7 (fig. 71), qui montrent la non-généricité de la transversalité des variétés stables et instables d’ensembles invariants hyperboliques et empêchent donc de considérer comme marginale la complexité associée à la présence de contacts homoclines ou hétéroclines. Citons également, pour les difféomorphismes, la densité des points périodiques dans l’ensemble non errant malheureusement démontrée dans le seul cas où r = 1 (closing lemma de Pugh, 1967).

Une classe particulière de difféomorphismes ayant la propriété de Kupka-Smale est formée par les difféomorphismes axiome A n’ayant pas de cycle; nous avons déjà dit au chapitre 7 qu’ils étaient 行-stables; en particularisant un peu plus, on arrive aux difféomorphismes axiome A ayant la propriété de transversalité forte (cf. chap. 7) qui sont, eux, structurellement stables et dont la dynamique est raisonnablement descriptible.

Ni la 行-stabilité ni, a fortiori, la stabilité structurelle ne sont des propriétés génériques des difféomorphismes (resp. des flots) d’une variété de dimension supérieure ou égale à deux (resp. trois). Les premiers exemples de systèmes dynamiques non approximables par un système dynamique structurellement stable (resp. 行-stable) ont été donnés par Smale en 1966 (resp. Abraham et Smale en 1970). On en connaît beaucoup aujourd’hui, par exemple les flots possédant des «attracteurs de Lorenz géométriques» (cf. chap. 6), dont la topologie peut même varier continûment en fonction de paramètres.

Des affirmations (beaucoup!) plus faibles que la généricité ont pu être démontrées à propos des difféomorphismes structurellement stables, par exemple leur densité pour la topologie C⁰ (Shub) et leur existence dans toute classe d’isotopie (c’est-à-dire composante connexe par arcs) de Diff^r (M) (Smale). En particularisant de plus en plus, on obtient les difféomorphismes de Morse-Smale ( 行(f ) union finie de points périodiques hyperboliques dont les variétés stables et instables se coupent transversalement), génériques en dimension un, mais qui en plus grande dimension n’existent même pas dans toutes les classes d’isotopie (travaux de Shub, Sullivan, Franks).

Notons que, comme les champs de gradients qu’ils généralisent, les systèmes dynamiques de Morse-Smale vérifient des inégalités de Morse reliant leurs points singuliers et leurs orbites périodiques (resp. leurs points fixes et leurs points périodiques) à la topologie de la variété ambiante et sont donc un bon outil d’investigation de cette topologie. Quant à la recherche de «jolis» représentants d’une classe d’isotopie de Diff^r (M), elle a fait en particulier l’objet de beaux travaux de Thurston (difféomorphismes pseudo-Anosov sur les surfaces).

Des résultats de généricité existent également à l’intérieur de classes particulières de systèmes dynamiques, en particulier les champs de gradients (travaux de Palis et Smale) et les champs hamiltoniens (travaux de Takens et de Robinson). Rappelons cependant que le théorème de K.A.M. montre que l’ergodicité dans les surfaces d’énergie constante n’est pas une propriété générique des systèmes hamiltoniens dès que la classe de différentiabilité r est assez grande, et ce bien que l’ergodicité soit générique pour les homéomorphismes d’une variété compacte de dimension supérieure ou égale à deux qui préservent le volume (théorème d’Oxtoby-Ulam; le problème est ouvert pour les difféomorphismes une fois continûment différentiables).

Beaucoup moins de choses sont connues si l’on remplace la considération d’un système dynamique individuel par celle de familles paramétrées (cf. comparer à la théorie du déploiement versel décrite dans l’article SINGULARITÉS DES APPLICATIONS DIFFÉRENTIABLES). Sotomayor a pu caractériser une classe de familles à un paramètre génériques de champs de vecteurs sur les surfaces, mais l’exemple des champs sans singularité sur le tore montre que l’ensemble de bifurcation de la famille (ensemble des valeurs du paramètre pour lesquelles le champ n’est pas structurellement stable) est «génériquement» compliqué; cela est vrai même pour des familles «stables» (en un sens qu’on laisse au lecteur le soin de deviner), ainsi que l’on montré récemment Guttierez, Lins et de Melo en étudiant sur le tore certains champs de vecteurs avec singularité dont la construction remonte à Cherry.

Des phénomènes analogues ont été mis en évidence par Newhouse, Palis et Takens dans les familles à un paramètre de difféomorphismes qui «sortent» du sous-espace des difféomorphismes de Morse-Smale.

La description de familles à un paramètre génériques d’équations différentielles, de difféomorphismes, ou même d’endomorphismes, dans lesquelles la structure des attracteurs devient de plus en plus complexe, est actuellement l’objet d’une grande activité sous le nom recherche de routes possibles vers la turbulence : Il semble en effet que l’établissement de la turbulence hydrodynamique puisse être modelé par de telles familles en petite dimension. Les premiers modèles ont été ceux de Lorenz, puis de Ruelle-Takens (suite de bifurcations de Hopf menant à un attracteur torique de petite dimension et apparition d’un attracteur étrange, cf. chap. 6 et la discussion dans le chap. 5 du problème de la disparition des courbes invariantes des difféomorphismes du plan). Le phénomène de l’accumulation de dédoublements d’orbites périodiques, auquel on a fait allusion au chap. 6 à propos des endomorphismes de l’intervalle, semble également très pertinent (expériences de l’équipe de Libchaber). Une façon simple de le produire est de modifier la position du point de suspension d’un pendule avec une période moitié de la période propre du pendule, le paramètre étant l’amplitude de la modification.

Pour des familles à deux (ou plus de deux) paramètres, les choses se compliquent encore. Takens vient par exemple de montrer (1983) que la stabilité structurelle n’est pas une propriété générique des familles à plus de huit paramètres de champs de gradients sur une variété de dimension supérieure ou égale à trois, alors qu’elle l’est pour les familles à un paramètre (Palis et Takens) et que l’absence de récurrence dans de tels champs aurait pu faire croire qu’elle l’était quel que soit le nombre de paramètres. Dans le cas des difféomorphismes des surfaces, il faut noter le rôle de centre organisateur des difféomorphismes conservant les aires (hamiltoniens): les déformations d’un difféomorphisme qui possède un point fixe elliptique font apparaître des ensembles de bifurcation dont la structure est à deux niveaux: un niveau robuste du type «théorie des catastrophes» auquel se superpose un niveau fin très instable dont la complexité est analogue à la complexité hamiltonienne décrite au chap. 4.

Autrement dit, c’est toute la structuration des espaces 﨑^r (M) ou Diff^r (M) qui doit être décrite à plusieurs niveaux, au fur et à mesure que l’on prend en compte des phénomènes plus fins, et une théorie générale du déploiement versel semble aujourd’hui très utopique.

La signification de tout cela, et en particulier de la non-généricité de la stabilité structurelle, pour la physique n’est pas très claire: il semblerait naturel que d’une équation (ou d’une famille d’équations) décrivant un «phénomène» on exige une certaine stabilité, au moins dans la classe des équations à laquelle elle appartient naturellement; les coefficients numériques qui interviennent ne sont en effet connus qu’avec une certaine approximation et il serait raisonnable que le comportement qualitatif de la dynamique ne soit pas trop sensible aux approximations. Cependant, ce point de vue «évident» doit être tempéré pour plusieurs raisons. Tout d’abord, notre notion de stabilité fondée sur l’équivalence topologique des dynamiques est sans doute trop fine et non pertinente: les expérimentateurs «voient» tous les jours des phénomènes quasi périodiques qui, au moins dans les systèmes dissipatifs, sont très instables, bien que fréquents au sens de la mesure, mais comment distinguer une orbite dense sur un tore d’une longue orbite périodique (ou encore comment distinguer un rationnel d’un irrationnel)? Il serait important également d’introduire dans la théorie diverses échelles comme cela se fait depuis toujours en physique (Harthong a avancé que l’analyse non standard était un bon outil pour cela).

En effet, comme le remarque Thom, la plupart des bifurcations d’un système complexe n’ont sans doute aucun effet visible sur la phénoménologie qu’il sous-tend, du moins à une échelle donnée.

Enfin, une expérience fournit en général des suites de mesures et c’est sur de telles suites qu’il faut juger de la dynamique sous-jacente éventuelle. D’intéressants résultats de Takens ont paru récemment, qui valident et prolongent un certain nombre de pratiques utilisées en physique et permettent, dans les cas favorables, de prouver le caractère «déterministe» d’un processus et même de calculer approximativement la dimension de Hausdorff et l’entropie topologique d’un attracteur. À la base, se trouve une jolie démonstration du classique théorème de Whitney sur le plongement des variétés dans l’espace euclidien: si M est une variété compacte de dimension n , f : MR une fonction différentiable, X un champ de vecteurs sur M de flot ( 﨏^t )^t size=1^{捻 R}, 見 un nombre réel non nul, k un entier supérieur ou égal à 2n + 1, et si les données (X, f , 見) sont «génériques», l’application F de M dans R^{2n +1} définie par F(x ) = (f (x ), f ( 﨏 size=1^見(x )), ..., f ( 﨏^(k size=1^漣1) size=1^見(x ))) est un plongement (application injective dont la dérivée est en tout point de rang maximal). Pour représenter un attracteur dont on pense qu’il habite un espace de petite dimension, on peut donc essayer l’opération suivante: on mesure à des intervalles de temps 見 l’observable représentée par f , ce qui fournit une suite a ¹, a ², a ³, ..., a ^p , ... de nombres, et on construit dans R^k (en général k = 2 ou 3) l’ensemble des points de la forme (a ^{i +1}, ..., a ^{i +k} ), i = 0, 1, 2, ... Les résultats sont propres à étonner.

Que retenir des divers exemples évoqués, sinon un mélange d’admiration devant la beauté de certaines percées théoriques et de modestie devant le foisonnement des situations non maîtrisées.

Seule la classe des systèmes «hyperboliques» se prête pour le moment à une théorie unifiée conduisant à une bonne connaissance de la dynamique, et ce grâce à la propriété d’approximation des quasi-morphismes par des morphismes, indissolublement liée aux propriétés de stabilité; mais quelles difficultés pour montrer qu’une partie au moins d’un système réel se trouve dans cette classe!

Les problèmes «elliptiques» nécessitent en général des méthodes plus fines au niveau analytique; on comparera par exemple les théorèmes affirmant l’apparition par bifurcation de Hopf d’une courbe invariante d’un difféomorphisme et les théorèmes du type Kolmogorov-Arnold-Moser montrant, dans le cas d’un difféomorphisme du plan conservant les aires, l’existence de courbes invariantes de bon nombre de rotation au voisinage d’un point fixe elliptique: les premiers n’utilisent que l’existence d’un point fixe pour une contraction d’un espace de Banach, les seconds nécessitent dans leur généralité un théorème de fonctions implicites «raffiné» dans un espace de Fréchet.

Quant aux systèmes hamiltoniens «complètement intégrables» ^{[cf. MÉCANIQUE ANALYTIQUE]}, qui jouent dans l’univers hamiltonien le rôle tenu par les systèmes linéaires dans la théorie dissipative, leur classification n’est que très partielle et la construction sporadique de nouveaux exemples, par des méthodes en général algébriques, reste bien mystérieuse.

Bien entendu, dans les problèmes «réels», l’hyperbolique et l’elliptique se mélangent inextricablement et des méthodes adaptées doivent être péniblement élaborées: lorsqu’il s’agit d’équations différentielles en dimension trois, on arrive quelquefois à réduire la dimension en se ramenant (essentiellement) à un difféomorphisme en dimension deux par le choix d’une section de Poincaré (ou quelque chose d’approchant) puis à un endomorphisme, éventuellement discontinu, en dimension un si par bonheur ledit difféomorphisme laisse invariant un feuilletage, enfin à la dimension zéro par le biais de la dynamique symbolique d’un tel endomorphisme.

Ce schéma est à peu près suivi dans l’étude des attracteurs de Lorenz géométriques, mais la situation n’est pas toujours aussi idyllique: l’exemple de l’attracteur de Henon, dont on ne sait toujours pas si c’est un attracteur étrange ou l’adhérence d’une réunion d’orbites périodiques de périodes arbitrairement grandes, est significatif de la difficulté qu’il y a à appréhender géométriquement une situation donnée analytiquement.

Et pourtant, sans nier l’indéniable apport conceptuel de la théorie générale à des questions telles que stabilité, généricité, déterminisme, relation local-global, catastrophes..., ce sont aussi les équations «naturelles» qui nous intéressent, et pas seulement celles qui habitent un espace de dimension finie. Nous avions commencé en citant, avec Lefschetz, deux grands modèles d’équations différentielles: les équations de Van der Pol du côté dissipatif, celles du problème des trois corps du côté conservatif; nous finirons en rappelant les problèmes posés par des équations aux dérivées partielles telles que celles de Navier-Stokes (cas dissipatif) et Kortweg-de Vries (cas conservatif) (cf. équations aux DÉRIVÉS PARTIELLES - Équations aux dérivées partielles non linéaires). Ces équations peuvent raisonnablement être considérées comme des «équations différentielles» en dimension infinie, ce qui permet de poser à leur propos des problèmes de nature plus géométrique, tels que la dimension des attracteurs pour Navier-Stokes ou l’existence de solitons pour des perturbations conservatives de Kortweg-de Vries. Une littérature importante est consacrée à ces thèmes; pour s’en convaincre, il n’est que de remarquer le caractère familier des deux mots «turbulence» et «solitons» dans le vocabulaire scientifique d’aujourd’hui.

Certes, on ne peut toujours pas décider si, comme le demandait Feynman, l’odeur des violettes se trouve ou non dans l’équation de Schrödinger et la «turbulence développée» dans les équations de Navier-Stokes, mais certains aspects des équations sont venus à la lumière qui ont bouleversé nos intuitions; en particulier, la richesse des comportements possibles d’une équation différentielle en dimension supérieure à deux incite à croire que les équations aux dérivées partielles de la physique contiennent «beaucoup» de phénomènes.

En ce qui concerne les méthodes, on est loin de la belle unité que présentent certaines théories plus formalisées.

Certains s’en désoleront, d’autres penseront, et c’est l’opinion de l’auteur, que c’est dans ce foisonnement un peu hétéroclite que se manifeste le mieux la richesse et le «naturel» de ce domaine, et donc aussi son charme.

Encyclopédie Universelle. 2012.